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FOREWORD 

This volume comprises the Final Report prepared in 
complete fulfillment of the requirements of Contract 
NAS 7-434, Phase II. The Final Report covers the 
period April 1, 1967 to April 1, 1968. Contributors 
for the technical work reported are as follows: 

J. A. Findlay 
A. J. Orsino 
H, J. Sneck 

Significant contributions were made by J. A. Reilly, 
in computer programming, data reduction, literature 
review and editorial and publication assistance, and 
by H. T. Robison, in carrying out a successful 
experimental program and assisting with the literature 
review. 
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I. SUMMARY 

The final report consists of a single volume and presents the results 
and recommendations of the technical work performed under the National Aeronautics 
and Space Administration Contract NAS7-434, Phase II, from April 1, 1967 to 
April 1, 1968. 

The technical work is of a fundamental nature intended to establish a 
basic understanding of seals and sealing mechanisms. Analytical as well as 
experimental efforts are described in the area of Newtonian and non-Newtonian 
sealing concepts. Viscoelastic experiments using a seal-like configuration were 
continued. A hydrostatic gas bearing was used in measuring torque in these 
tests. The data obtained from these tests was used in the attempt to determine 
a simple correlative model for non-Newtonian phenomenon. Correlation included 
the shear-stress shear-rate correlation and the normal stress shear-rate 
correlation. 

Several series of Newtonian tests were run to study the effects of seal 
geometry on leakage and to consider the relationship of the different geometries 
to the phenomenon of "seal pumping". Corresponding analytical work was also 
carried out and through the use of digital computers, correlations were made and 
the derived theories were substantiated. The literature analysis and review 
was updated. 

All technical work was accomplished by personnel of the Research and 
Development Center, General Electric Company, Schenectady, New York. Names of 
individual contributions appear on the appendix containing the full report of 
their particular activity. Project Engineer is J. A. Findlay. Project Manager 
is G. R. Fox. The work was accomplished under the direction of Natipnal 
Aeronautics and Space Administration Technical Manager, Mr. R. S. Weiner, 

Liquid Propulsion Section, Jet Propulsion Section, Jet Propulsion Laboratory, 
California Institute of Technology, 4800 Oak Grove Drive, Passadena, California. 
National Aeronautics and Space Administration Project Manager is Mr. F. E. 
Compitello, Liquid Propulsion Technology, RPL, Office of Advanced Research and 
Technology, NASA Headquarters, Washington, D.C. 20546. 
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II. DESCRIPTION OF PROGRAM 

A, Introduction 

The Research and Development Center, General Electric Company, has been 
engaged in a fundamental study of dynamic and static seals for the National 
Aeronautics and Space Administration under Contract NAS7-434, Phase II. The 
principal work of the study is concerned with the technology of seals and seal- 
ing mechanisms on an analytical as well as an experimental basis. The scope of 
the program is extremely broad and requires detailed knowledge of several tech- 
nical disciplines. In view of these particular characteristics of the program, 
it has been necessary to select specific areas for detailed study in order to 
assure the development of meaningful information and data. It is expected that 
the results of these technical efforts would serve to form a basis for formula- 
tion of seal design criteria. A considerable amount of preliminary work, 
intended to provide a measure of order in the available technology of seals and 
sealing mechanisms was accomplished on contract NAS7-102. The work of this 
contract draws heavily from the findings and results of contract NAS7-102 and of 
Phase I of the current contract. The technical results of this work are contained 
in one volume. , 

B. Program Objective 

The objective of this program is to advance the technology of dynamic and 
static seals for liquid rocket engines. A study of the fundamentals of sealing 
processes is to be carried out, in order to provide the understanding necessary 
to formulate design criteria for seals in future liquid rocket engines. 

C. Description of Program Plan 

The detailed program plan covering the technical effort under this contract 
is outlined below. The task designations are evolved from the contract work 
statement and include suitable sub-headings to facilitate description of the 
effort involved. 

Program Plan 

Task 1 - Non-Newtonian Face Seal Experiments 

- Effect of Land Width 

- Torque Measurements 

- Leakage Speed Test 
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Task 2 


- Non-Newtonian Correlative Model 

- Review Existing Models 

- Formulation of a Tentative Model 

- Suggest Experiments 

- Evaluate Data and Model 

Task 3 - Face Seal Geometry Effects Study 

- Cavitation Experiments 

- Eccentric Seal Experiments 

- Inertia Experiments 

- Two- fluid Study 

Task 4 - Seal Design 

- Correlate Seal Data 

- Evaluate Seal Theories 

- Extend and Improve theories 

- Design Data 

Task 5 - Literature Review 

- Updating 

The planned analytical studies and experimental investigations cover areas 
which heretofore have not been explored in significant depth. These areas were 
considered to possess a high potential for improved understanding of seals and 
sealing mechanisms. During execution of the program, some variations in original 
plans were necessary. This was somewhat anticipated since it is not always 
possible to predict the course of fundamental technical investigations. As the 
work progressed, the knowledge acquired was used to reevaluate, and if necessary, 
to redirect the subsequent effort for each task. In all cases, the overall 
program objective was followed as closely as possible to assure useful results. 

A discussion of the results and accomplishments for each task will be 
presented in a subsequent section of this report. 
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III. DISCUSSION OF ACCOMPLISHMENTS 

1. General 

The technical accomplishments of this study are reported in a single 
volume. The technical work reported was initiated April 1, 1967 and continued 
to April 1, 1968. 

The results of the technical work were highly encouraging, especially the 
Newtonian face seal tasks. The phenomenon of "inward pumping" was found to be 
predictable by a theory based on radial eccentricity and misalignment. This 
theory is discussed in Appendix D and includes the effects of gas cavities and 
two different fluids on either side of the seal. These theories are confirmed 
by the experiments reported in Appendix C as is the inertia theory of H. J. 
Sneck's reported in Phase I of this contract. The Newtonian face seal theories 
are extended in Appendix E and F. 

The very difficult task of developing a non-Newtonian face seal correlation 
is discussed in Appendix A and B. 

The seals literature review is brought up to date in Appendix G. 

2. Discussion of Tasks 

A. Task 1 - Non-Newtonian Face Seal Experiments 

The non-Newtonian face seal experiments conducted during this contract 
phase were actually an extension of experimental work previously performed. The 
extension was due primarily to the marginal correlation obtained with the 
original jet capillary experiments. During the initial face seal experiments 
it was found that sealing effectiveness was several orders of magnitude below the 
value predicted from the capillary experiments. 

In the non-Newtonian face seal experiments covered by this report an 
effort was made to provide data for the development of a suitable correlative 
model. Such a correlation might provide a better understanding of the behavior 
of non-Newtonian fluids in face seals. Knowledge of torque level at the seal 
interface was considered a vital part of any effort to establish a suitable 
correlation. Under this task effort a hydrostatic torque device was constructed 
to facilitate evaluation of seal interface torque under dynamic conditions. 

The experimental results provided data at several fixed clearances over a 
broad speed range to 12000 rpm. Tests were performed with JP5 +3.2$ poly- 
isobutylene, JP5 with no polymer additive and SF96 silicone oil of two different 
viscosities. Leakage data as well as torque, clearance and speed were obtained 
for all tests. No evaluation of the data was performed under this task effort. 
The correlation of test results was accomplished as part of a separate task of 
this contract phase. The test data is presented in Appendix A. 
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Bo Task 2 - Non-Newtonian Correlative Model 

The objective of this investigation was to determine if a simple correlative 
model could be developed with the aid of the data obtained under Task 1. While 
there is considerable viscometric data available for seal-like configurations, 
none of this data was obtained in the presence of an applied pressure gradient 
and non-zero leakage rates. Most viscoelastic fluid flow theories have been 
developed for the much simpler case of zero-net leakage in configurations which 
do not allow for applying pressure gradients.. 

It was found from an examination of the torque data obtained in Task 1 
that over the entire range of shear rates tested the tangential shear stress 
was described quite accurately by a form of the invariant power-law. Although 
it has been often suggested that the invariant form of the power-law might be 
applicable to three dimensional flows, this is one of the few instances where 
such a correlation has actually been made. 

As anticipated, the sealing effects of the normal stresses were found to 
be relatively unimportant at low shear-rates. This is because, as the results 
of no net leakage tests show, the normal stresses are related to the tangential 
shear-rate. 

At high tangential shear-rates where this shear-rate predominates over the 
radial shear-rate due to leakage, the normal stresses can again be correlated as 
a function of the tangential shear-rate. Unfortunately, this correlation is not 
the same as the no leakage-rate correlation at low shear-rates. 

In the intermediate range of tangential shear-rates where the radial and 
tangential shear-rates are of the same order of magnitude, it was not possible 
to generate a general correlation which accounted for the normal stress effects 
on leakage sealing. 

Attempts to superimpose simultaneously the effects of degrading viscosity 
and improving normal stress effects failed probably due to the inherent non- 
linearity of viscoelastic phenomena itself. Failure to obtain a correlative 
theory on a global basis, i.e., by not considering the details of the phenomena 
is disappointing but not unexpected. The results obtained here will be very 
useful in the future efforts to develop such detailed theories. 

C„ Task 3 - Seal Geometry Effects 

The experiments in this task have been keyed towards the accurate measure- 
ment of face seal geometry. The basic problem in obtaining data to substantiate 
the theories of Newtonian face seal operation has been the uncontrolled variations 
in the surface geometry. Through the use of relatively large film thicknesses 
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and optically flat glass as the seal face, the effects of uncontrolled surface 
variations was minimized. Previous literature concerning seal geometry effects 
have not dealt with such controlled test conditions and, thus, the data obtained 
under this task appears to be the most useful for correlating theories. 

The experiments carried out in this task include those with concentric 
and misaligned seals, eccentric seals and those tests studying the effects of 
inertia and of fluid viscosity on seal leakage. The results of the cavitation 
experiments showed that the gas cavities had no effect on the leakage. The full 
film results verified the fact that the effect of seal eccentricity can be 
added directly to the hydrostatic leakage. Both of these results compared very 
well with the existing theories and these specific correlations are presented 
in Appendix C. 

Tests were also conducted under this task which are concerned with the two- 

fluid problem in face seal operation. This phenomenon considers the changes in 

seal hydrodynamics which occur when the fluids on either side of the seal face 
are of different viscosities. These tests were run using viscosity ratios, 

ranging from 3.5 to 0.002. The data taken was correlated using this ratio. 
Further tests were carried out in this task which consider the inertial effect 
on leakage. In this series of tests, as was with all tests, GE Time Sharing 
Computer programs were used to reduce the data. The results of both the two 
fluid and the inertia tests can be seen in Appendix C. 

D. Task 4 

1. Seal Design I , 

The main effort of this task has been the extension of basic face seal 

theories. The "short bearing" theory developed in Phase I of this contract 

has been used to predict gas cavity shapes and inception points. The analysis 
has been extended from the wavy film thickness to the tilted seal case used in 
the laboratory experiments. 

A full film solution has been derived for the eccentric seal case using 
the "short bearing" approximation. This was done to show that a combination of 
seal surface waviness or tilt coupled with an eccentricity between the two seal 
surfaces can result in a pumping of fluid against a hydrostatic pressure 
differential. 

The GE 635 face seal computer program was used to compute seal leakage 
and pumping due to eccentricity in the misaligned seal. This program has been 
improved to give faster convergence and to compute the short bearing approxi- 
mation as well as the two dimensional solution. 
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A further analysis of the two-fluid case was conducted. Many of the basic 
relationships were derived. A summary of the equations and approach to their 
solution is given in Appendix D. 

2. Seal Design II and III 

The analyses of face seal behavior are invariably based on the assumption 
of isothermal flow through the seal clearance. This assumption allows both 
the density and particularly the viscosity to be taken as constants. The 
analysis of thermal effects in face seals presented here formalizes the pro- 
cedures whereby the designer may estimate when such simplifying assumptions are 
valid. 

Reversed flow effects in seals are particularly interesting when the 
externally applied pressure in a face seal exceeds the internal pressure. Under 
this condition it is possible to establish counter flows within the clearance 
space so that there is a net exchange of internal and external fluids. The 
rate at which these fluids are exchanged between the inner and outer boundaries 
is an important design consideration. Methods for predicting these exchange 
rates are developed in the reversed flow analyses. 

E. Task 5 - Literature Review 

The literature review was patterned after the similar study under Phase I 
of this contract. The review of all available reference sources was updated as 
far as possible with most sources reviewed up to January 1968. 

Abstracts were prepared for all the material considered of interest in the 
area of seals and sealing mechanisms for both dynamic and static seals. More 
than one hundred abstracts from the more than eleven major source references were 
prepared. In some cases the abstract was taken from the authors own prepared 
abstract. In other cases the abstract was prepared from a complete reading of 
the text. 
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APPENDIX A 


NON-NEWTONIAN FACE SEAL EXPERIMENT 
BY 

A. J. Orsino 

I. INTRODUCTION 

The experimental study by Bernd, Reference 1, predicted increased sealing 
effectiveness in a face type dynamic seal for fluids containing small amounts 
of certain polymer additive. As a result of the work reported in Reference 1, 
an experimental program using a seal like configuration was undertaken in an 
effort to obtain a more detailed understanding of non-Newtonian fluid behavior 
in dynamic face seals. The results of this effort are reported in Reference 2 
and show only a minimal degree of correlation with the results of the jet 
capillary experiments by Bernd, The fundamental reasons for the poor cor- 
relation is probably due to the significant differences in geometry between 
the capillary and the dynamic face seal. Transit time for the fluid in the 
capillary test is extremely short and perhaps of a lower value than the fluid 
relaxation time. This is not the case for the seal configuration. The 
comparatively tortuous path for the fluid in the seal, with its correspondingly 
longer transit time, may give rise to increased possibility for rupture of 
intramolecular bands. Rupture of fluid molecular bands has a strong influence 
on fluid properties. 

In order to better understand the behavior of non-Newtonian fluids in thin 
films, and thereby establish a basic understanding of face seals operating with 
such fluids, a new series of experiments was conducted. The objective of this 
new testing was the acquisition of experimental data for use in evolving a 
suitable correlative model. The model would provide a basis for making 
meaningful predictions regarding face seal performance. Testing during this 
contract phase was an extension of previous work using the identical test 
equipment plus some modification to allow for torque measurements. It was 
expected that torque data would provide an independent means for development 
of a suitable correlative model. Knowledge of the torque present at the 
fluid interface allows detailed evaluation of viscosity and its variation with 
shear rate. 
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This Appendix describes the design of the torque device, the test pro- 
cedure used and the results of the experimental data obtained for the several 
fluids tested, A description of the data evaluation and a discussion of the 
correlative model evolved are contained elsewhere in this report under the 
results of the Task 2 effort* 

II. HYDROSTATIC TORQUE MEASURING DEVICE 

The single most important feature of the test device used for the non- 
Newtonian seal experiments is the hydrostatic bearing arrangement designed to 
provide a high stiffness characteristic and a low parasitic torque. Previous 
designs considered a set of preloaded angular contact ball bearings with an 
oil mist lubrication system* However, calculations for a twenty millimeter 
duplex angular contact ball bearing with sufficient preload to assure a 
1°10^ Ib/in stiffness, revealed a bearing internal torque at least one order 
of magnitude greater than the anticipated interface fluid torque. An added 
limitation on such a bearing would be the torque "ripple” due to variations 
in bearing surface finish as well as surface waviness. To circumvent the 
limitations of the rolling element bearing arrangement a hydrostatic fluid 
film bearing was designed. Pressurized gas was used as the bearing working 
fluid in order to preclude the possibility of contamination of the seal test 
fluid* The risk of pneumatic instability as well as lack up can be very 
troublesome in such bearings and, therefore, some caution is required to 
assure satisfactory operation. The arrangement selected for detail design 
featured a double acting thrust bearing combined with a journal bearing having 
a double row of admission holes. The design of the hydrostatic torque device is 
very conventional and is based on the criteria contained in Reference 3 and 4. 

Figure 1 presents the overall seal assembly test device and shows the 
bearing arrangement in cross-section. The upper thrust bearing is supplied by 
36 holes of 0.016" diameter equally spaced and supplied from a common manifold. 
The lower bearing is likewise supplied by 36 holes of 0.016 diameter equally 
spaced. Alternate gas supply ports for the lower thrust bearing are connected 
to the journal bearing which is fed by two rows of 18 holes each. The journal 
bearing holes are also 0 o 016 n diameter. By this arrangement it is possible to 
supply the journal bearing and the lower thrust bearing from one common 
manifold* The two manifolds were supplied from separate pressure sources. 
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Torque measurements were made by conventional strain gauge technique 
employing a bridge circuit connected to an Ellis amplifier. Linearity of the 
strain gauge and torque arm arrangement is checked by dead weights. During 
operation a fixed amplifier gain is maintained to further assure system 
linearity. Speed was monitored on a Hewlett Packard counter connected to a 
magnetic pickup which received its signal from three pins mounted on the 
spindle shaft. 

Seal interface clearance was adjusted by two sets of opposed screws which 
allow accurate dependable setting in the low: microinch range. Determination of 
actual clearance is made by three "Electro-Jet" Sheffield magnetic reluctance 
probes mounted in close proximity to the interface. The "Electro-Jet"Hus 3 
direct contact instrument capable of providing accurate measurement in the low 
microinch range. Provision was made for use of "Wayne Kerr" capacitance probes, 
however it was not possible to use them at high speed due to the difficulty of 
maintaining a full fluid film under the probe tip. A partial film gives rise 
to variations in capacitance and consequently erratic and uncertain data. The 
Sheffield probes are equipped with a retracting mechanism to prevent damage 
to the runner surface. 

As for previous tests of this type a graduated tube was used as the 
pressurized fluid reservoir connected to the center of the stationary head by a 
flexible hose. An overall view of the non-Newtonian Test Assembly is shown 
in Figure 2, The Sheffield meter located next to the graduated reservoir is 
shown connected to the distance probes. Other related instrumentation can be 
clearly identified. Figure 3 is a close-up view of the hydrostatic torque 
measuring device. Flexible connections to the manifolds can be seen as well 
as the connecting plastic hose from the reservoir to the seal center position. 
The torque lever and strain gauge arrangement is clearly visible. Component 
parts for the hydrostatic test device including the stationary seal mating 
face are shown in Figure 4. 

III 0 TORQUE DEVICE PERFOjRHANCE TEST 

To assure adequate performance of the torque measuring device, a series 
of calibration tests were performed to determine the stiffness characteristics 
of the thrust bearings. This was accomplished by applying dead weight loads 
directly to the bearing and subsequently measuring the deflection in the 
direction of loading. Figure 5 is a plot of the load versus displacement for 
the upper thrust bearing at several selected pressure levels. Manifold 
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pressures were set rather arbitrarily but were based on equal clearance for 
each, thrust bearing. From Figure 5 it can be seen that the required stiffness 
of l.?10 6 Ib/in can be easily obtained with a pressure level of approximately 
100 PSIG. For modest increases in pressure a stiffness as high as 2°10^ lb/in 
can be obtained to match the stiffness of the support structure. Data in 
Figure 5 was taken on the inverted unit for convenience of loading with dead 
weights. 

Total axial clearance in the thrust bearing is 0.00074 inches and is 
sufficient to provide the necessary face clearance during all conditions of test 
loading. For example, at a seal interface pressure of 40 PSIG a total axial 
load of approximately 95 pounds can be developed. With a supply manifold 
pressure of 100 PSIG the displacement will be limited to approximately 100 
microincb.es. At the higher supply pressure of 140 PSIG the displacement is 
limited to less than 50 microinches . The actual axial movement is not critical 
providing sufficient clearance is maintained under load in the assembly. 

Clearance settings were made with pressure at the seal interface so that total 
axial movement of the hydrostatic torque device was properly taken into 
account. As part of the torque device performance testing the system breakaway 
torque was determined with an ARPO hand type torque meter. A breakaway torque 
of 0.22 oz inches was found and it was judged that steady state drag torque 
did not exceed approximately 0.04 oz inches. 

IV. DESCRIPTION OF TEST PROCEDURE 

As an initial step, a preliminary interface clearance was set without 
fluid in the reservoir. The final setting was accomplished with the reservoir 
at the proper pressure for test and with zero speed. It was necessary to 
establish the correct fluid level after the final clearance setting. The 
Electro-Jet gauges were retracted for all speed tests. 

With the test fluid under pressure the desired speed was established and 
data was recorded. Torque,, speed, temperature and flow data were taken. The 
flow was determined by timing the fluid level change for a fixed fluid volume 
change. Immediately after each test the speed was reduced to zero and with the 
reservoir under pressure the clearance was determined with the Electro-Jet gauges 
The running clearance then became the average of the readings before and after 
test. This correction accounted for any heating which developed during test 
within the drive spindle and was of only second order of significance. Minimum 
thermal expansion was present in the seal mating parts since they were made 
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from invar. No change in interface clearance was present due to speed alone. 

This had been determined independently with the capacitance probes operating 
without interface fluid. 

V. LOW SHEAR RATE VISCOSITY TEST 

Viscosity data, employing a conventional Cannon and Fenske glass capillary 
viscometer, was obtained for samples of JP5 containing 3.2$ Vistanex polymer 
additive. The results indicate a viscosity of 4471 centistokes (20,650 Saybolt Secon 
Seconds) at 77°F. This value of viscosity is approximately 2200 times greater 
than the viscosity of the. base JP5 fluid. In effect this data provides the 
zero shear rate viscosity value and agrees quite well with data for the same 
fluid taken on the Weissenberg Rheogoniometer . 

Additional viscosity data was obtained from the rheogoniometer located at 
the General Electric Research and Development Center. The data for JP5 with 
3.2 $ Vistanex of molecular weight 2 x 10^ and 10^ are shown in Figure 6. The 
data covers only the low shear rate range. Also shown is the normal stress 
for the Vistanex 200 mixture. This data was obtained in an effort to provide 
additional information useful in establishing the correlative model of Task 2. * 

VI. DISCUSSION OF TEST DATA 

Tests were conducted at nominal clearances of 0.001, 0.002, 0.003 and 
0.004 inches with JP5 plus 3.2 $ polyisobutylene. The major portion of the 
testing was conducted with a 1/2 inch wide interface land as used in prior 
testing. Tests with a 1/4 inch wide land were also completed. Data for the 
JP5 plus 3.2 $ polyisobutylene is contained in Tables 1, 2, 3 and 4 for each of 
the nominal clearances and at several pressures. The data is plotted in 
Figures 7, 8, 9 and 10 and shows the same trend as previously seen with other 
viscoelastic fluids. This is characterized by a change in rate of torque 
increase as speed is increased. In each case the change occurs between 1000 
and 1500 RPM. This slope variation is due to the influence of shear rate on 
viscosity. The specific relationship will be analyzed in greater depth as part 
of Task 2. The data covers testing to 12000 RPM. In the lower speed range, 
below 100 RPM, the torque readings contain some scatter due to the low level 
of interface torque. This difficulty was recognized at an early date and 
provision was made to obtain low shear rate data from a separate source. The 
low shear data was obtained on a rheogoniometer available at the General 
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Electric Research and Development Center. The data at the lower shear is shown 
in Figure 6. This low shear rate data will be used as an additional input to 
the correlation study of Task 2. It should be noted that the rheogoniometer 
is limited to shear rates of 10 reciprocal seconds. 

As a further comparison Table 5 presents data for plain JP5 at a reservoir 
pressure of 5 PSIG. Due to the extremely low viscosity it was possible to test 
at a clearance of 0.001 inches only. In Figure 11 the torque verses speed from 
Table 5 is plotted. The straight line relationship characteristic of such a 
fluid is evident. 

The data from the tests with the narrow land is shown in Table 6 and Table 
7 for the 0.002 inch clearance and the 0.003 inch clearance respectively. The 
distinct change in slope is seen for each clearance in Figures 12 and 13 for 
the same clearances. 

As an additional source of information for correlations of the viscosity 
data, several tests were performed with silicone fluids. Silicone SF96 of two 
different viscosity values (50 and 100) were tested up to 1000 RPM. Fairly 
good correlation was achieved between measured torque values and torque 
values computed from flow rates. Test data for the silicone fluids evaluated 
is shown in Table 8. 
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■Newtonian Seal Test Device (see parts list. Table 9) 
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Figure 2. Non-Newtonian Test Assembly 
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Figure 3. Hydrostatic Torque Device 




Figure 4. Hydrostatic Torque Device (component parts) 
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TABLE 1 


JP5 + 3.2$ VISTANEX 200 
NOMINAL CLEARANCE 0.001 IN. 


(STANDARD 


INITIAL TEST 


Q 

h 

RPM 

in^/sec 

in. 

0 

.00218 

.00114 

200 

.0048 

.00114 

300 

.00451 

.00114 

700 

.00188 

.00113 

1000 

.00204 

.00113 

2000 

.0019 

.00108 

4000 

.0025 . 

.00103 

6000 

.00277 

.00097 

8000 

.00371 

.00086 

10000 

.00389 

.00077 

12000 

.00399 

.00070 


LAND WIDTH) 

RESERVOIR PRESSURE 20 PSI 
AMB. TEMPERATURE 77°F 


Torque 

oz-in 

1.65 

2.55 

4.5 
4.8 
6.45 
8.1 

9.6 
10.8 
11.85 
13.05 


REPEAT TEST 


0 

.0015 

.00111 

-- 

200 

.00383 

.00113 

1.42 

300 

.0040 

.00112 

1.8 

700 

.00191 

.00111 

3.83 

1000 

.00164 

.00110 

5.34 

2000 

.00156 

.00109 

6.83 

4000 

.00183 

.00105 

7.8 

6000 

.00222 

.00100 

9.16 

8000 

.00297 

.00095 

10.65 

10000 

.00356 

.00087 

11.86 

12000 

.00361 

.00078 

13.07 
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JP5 + 3.2 $ VISTANEX 200 

NOMINAL CLEARANCE 0.002 IN 

(STANDARD LAND WIDTH) 
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6000 .00383 .00390 4.27 

8000 .00562 .00385 3.45 


TABLE 5 


JP5 

NOMINAL CLEARANCE O.OOl IN 
(STANDARD LAND WIDTH) 

RESERVOIR PRESSURE - 5 PSI 
AMB. TEMPERATURE 80 F 



Q 

h 

Torque 

R PM 

itr/sec 

In. 

Oz-In 

0 

.0289 

.00113 

-- 

200 

.0294 

.00107 

.045 

300 

.0278 

.00107 

.0526 

700 

.0279 

.00105 

.12 

1000 

.0266 

.00104 

.233 

2000 

.0277 

.00104 

.51 

0 

.0259 

.00109 

— 

1000 

.0271 

.00106 

. 308 

2000 

.0232 

.00105 

.413 

4000 

.0258 

.00102 

.99 

6000 

.0294 

.00095 

1147 

8000 

.0367 

.00089 

1.89 

10000 

.0402 

.00081 

2.4 

12000 

.0281 

.00070 

3.00 
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TABLE 7 

JP5 WITH 3.2 $ VISTANEX 200 

NOMINAL CLEARANCE 0.003 IN 

(STANDARD LAND WIDTH) 

RESERVOIR PRESSURE - 10 PS I 
AMB. TEMPERATURE 75 °F 


RPM 

9 

in 3 /sec 

h 

in. 

Torque 

Oz-In 

0 

.0563 

.00306 


200 

.0728 

.00306 

.675 

400 

.0789 

.00306 

.90 

600 

.0823 

.00306 

1.27 

800 

.0823 

.00306 

1,43 

1000 

.0823 

.00306 

1.65 

2000 

.0104 

.00303 

3.3 

4000 

.00992 

.00303 

4.28 

6000 

.0118 

.00298 

4.87 

8000 

- - 

.00292 

5.4 

10000 

.0203 

.00290 

6.15 

1000 

.0842 

.00296 

1.65 

1500 

.0315 

.00294 

2.4 

2000 

.0115 

.00292 

3.15 
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NOMINAL CLEARANCE 0.002 IN 
(STANDARD LAND WIDTH) 
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TABLE 9 

PARTS LIST 
(Drwg. No. 422D466) 

NON NEWTONIAN SEAL TEST DEVICE 


Part No. 

Description 

1 

Mounting Plate 

2 

Journal Bearing 

3 

Journal 

4 

Manifold Ring 

5 

Thrust Bearing 

6 

Thrust Bearing Manifold 

7 

Spaur Ring 

8 

Torque Arm Support 

9 

Seal 

10 

Seal 

11 

Seal 

12 

Seal 

13 

Sealbolt 

14 

Connector 

15 

Screw 
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APPENDIX B 

NON-NEWTONIAN CORRELATIVE MODEL 
BY 

H„ J„ SNECK 

I. INTRODUCTION 

Viscometric tests of Newtonian fluids containing viscoelastic additives 
has shown that it is possible to generate a radial pressure gradient between 
two rotating disks such that the pressure at the center of the disks is larger 
than at the perimeter. If such a pressure gradient can also be created in a 
face seal, this pressure gradient could have a beneficial effect cm the outward 
leakage rate. 

Early experiments (Ref. 1) indicated that such a beneficial effect can be 
obtained. Further tests have been performed, this time including torque 
measurements, with the view to obtaining if possible a simple correlative model 
which describes the seal performance characteristics in terms of measurable 
fluid properties. The results of this attempt are reported on the following 
pages . 
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II 0 CALIBRATION TESTS USING A NEWTONIAN FLUID 

An initial check-out of the test-rig performance was made using several 
known Newtonian fluids. The fluids used were JP-5 and two silicone fluids, 
SF96(50) and SF96(100) o 

Because JP-5 has such a low viscosity it was not feasible to make tests 
with an applied pressure differential above 5 psi. The viscosity calculated 
from torque measurements was found to be constant at 1„5(10“7) reyns over the 
shear-rate range 10 ^ to 10^ sec“^ 0 This viscosity was then used to compute the 
leakage rate from the equation 


Q = 


ttC~ 


6p,l 


M 


(pi-p 2 + |q p(R^) 2 


n i 



The agreement between measured and predicted leakage-rate was within 10$, which 
is an acceptable correlation considering the effects which small clearance 
measurement errors can introduce,, (The nominal clearance at 5 psi was 0,001 
inches <,) 

The two silicone fluids were also tested at pressures of 10, 30, and 40 psi 
using a nominal gap of 0.002 inches. Torque measurements yielded viscosities 
of 29 centipoise for SF96(50) and 69 centipoise for SF96(100)« As with the 
JP-5 the calculated leakage was somewhat larger than the measured values, the 
maximum error in this case being of the order of 15$ for shear rates ranging 
from 5 x 10^ to 5 x 10^ sec“^ 0 

The measured viscosities are 60 to 70$ of the manufacturers indicated 
viscosity and about 80 to 90$ of the values required to compute the measured 
leakage-rates. The Newtonian fluid check-out of the apparatus indicates that 
reasonably accurate viscosities can be obtained from the torque measurements, 
and what is more important, these viscosities are constant throughout the 
shear-rate range and consistant with regard to leakage-rates. 


III. ANALYTICAL FORMULATION OF THE PROBLEM 

The equations of motion for the fluid are 



dr 


(P-Prr) + 


^rr-^00 , ^ T zr 

— — — -b — — 

r dz 


(1) 
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( 2 ) 


0 — — (p-Pzz) 


0 = 


dT z 9 

3z 


(3) 


Where the P rr , Pqq , and P zz are the deviatoric components of the stress. The 
last equation indicates that the shear stress is constant across the clearance. 
For laminar flow this suggests the approximation 


z 

c 


V = m> I 

Integrating the first equation across the clearance yields 


zr 




V- P rr 


dz 


z2dz + f ( r .6) 

where P = p-P rr . Assuming that the viscosity isotropic allows the shear- 
stresses to be related by 


(4) 

(5) 

( 6 ) 


T “ T a 
zr z9 


du 




dz/dz I 


(7) 


where t z q is independent of z. 

Integrating the radial momentum equation once, again with respect to z gives 


rw 

T z6 7" 


|zf(r,6) + H + 



p -P 
r Q9 rr 


dfdz - 


o o 


12 


+ g(r,9) 


( 8 ) 


If it is assumed that P90-P rr is independent of z, then 

2 I 


in 

o o 


p ea- p rr 


d^dz. 


2 

£1 6Q rr 


(9) 


Applying the boundary conditions that u = o at z = o, c allows the arbitrary 
functions f(r,9) and g(r,0) to be evaluated. As a result 


ru) 

iC 


— p _p " 

dP . 80 rr 

z^-cz) 

£ IM 2 

4 3\ 

Z -Cz 

dr + r 

. . 

1 z 

r Ic ) ! 

12 / 


( 10 ) 
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The volumetric leakage rate is calculated from 

Q = 

o 

When this flow-rate is evaluated it becomes 

Q T z0 |"dP . ^99“^rr c^ . 1 ( 3 \ p |raJ\2 5 

= • [a? + 7 j 12 + 12 ITol r Irl c 


/ 2nurdr 


(ID 


( 12 ) 


Integrating this equation with respect to r using the boundary conditions 

P = @ r = Ri (13) 

P = P 3 @ r = R 3 


yields 




dr 



dr 

(14) 


The left hand side of this 
rR2 

<7 V* 

dr 


62 _ 

3 

TTC 


zr 


J 

*1 



fdv ' 

r 

[dz , 


equation may also be written in the form 


(15) 


IV. NO-ROTATION LEAKAGE CORRELATIONS 

For the purpose of determining the basic behavior of the fluid the seal 
may be considered as a slightly divergent rectangular slit. When there is no 
rotation the slit may be considered analogous to a circular pipe in that the 
flow is caused only by a pressure differential. Rabinowitsch and Mooney (1, 2) 
have devised a general expression relating pressure drop, leakage-rate and 
geometry which is independent of fluid properties. A log-log plot of the 
parameters in this expression is shown in Figure 1 for the no-rotation test 
results. for 3.2$ Vistanex 200 in JP-5. These results are for pressure 
differences between 5 and 30 psi and clearances of 1, 2, 3 and 4 mils. It 
may be concluded that the fluid follows a power-law relation from the fact 
that the results generally lie along a straight line in this plot. 

The quantity 


is related to the radial shear stress T zr for parallel disks by the 



The data shown in Figure 1 is correlated using the expression 


C 1V P 2 

V R 1 


= 2.86 (10“ 3 ) 



If it is assumed that 



the last two equations yield the power law expression 


T = 1.222 ( 10 “ 3 ) 
zr v 



\c r 


0.6 


3u 

dz 


where the viscosity can be taken to be 


p = 1.222 (10' 3 ) 


_2_r0.6 
c r| 


(reyns) 


V. RESULTS OF TORQUE TESTS 

The torque test data can be reduced assuming that 


t 2 0 "'H- = M- 


rw 


z 9 


in the absence of actual velocity profile data. If p, is constant, 
torque is given by 



expression 


( 16 ) 


( 17 ) 


( 18 ) 


( 19 ) 


( 20 ) 


( 21 ) 

then the 


( 22 ) 
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The viscosity calculated from the data using this equation is shown in 
Figure 2 as a function of 

* £ . 

c 

There is considerable scatter in the results up to about 


Rjd 

__Z_- 5(i 0 5)- 
c 

after which there is a good correlation for all applied pressure differences 
and clearances. In this range it is obvious that the viscosity follows a 
power-law relationship of the form 

in 


ru) 

“-“or 


(23) 


If the torque is now calculated on the basis of the power-law viscosity 
relation, the equation for the torque becomes 


T = 2 


f p. oun TT \ 1 < 

u\ n 

■ n+4 n+4* 

R 2 " R 1 

1 c H 

cl 

n+4 


(24) 


A first estimate of n can be obtained from Figure 2 and then the test data 
processed through the revised form of the torque equation to determine the 
constant p- 0 . Carrying out this calculation yields the correlation 

,-0.60 

-I 

c 


.-0.60 

p. = 2. 7(10“ 3 ) (t-J (reyns) 


(25) 


for 

— > 5(10 4 ) 
c 

VI. INVARIANT POWER LAW CORRELATION 

The scatter in the low shear rate viscosity is undoubtedly due to the fact 
that the radial and tangential shear rates are of the same order of magnitude 
in this range. At higher rotational speeds, the tangential shear rate 
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predominates so that the simple correlation given above is quite accurate. The 

invariant form of the power-law is . 

n- 1 

„ I T" 

2 . 2 


T ij = ^o 


+ 


dv 

dz 


where the viscosity is given by 


M- * V> t 


!i | 2 + (K ) 2 




( 26 ) 


(27) 


From the two previous correlations it is clear that n-1 = -0.6. Since 

dv J ru) | 

dz 'c ' 

the quantity p, o is equal to 2.7 (10"3). The quantity Q/c^r used in equation 
(20) is not equal to du/dz but only a measure of it. Bearing this in mind, the 
discrepancy in coefficients can be accommodated by writing the viscosity as 



-0,30 


( 28 ) 


If k is chosen equal to 3.75, then this expression also agrees with the no- 
rotation correlation for (U = o. 

The final expression for the viscosity is then 


H- 


2.7 


h’.' (“I 


0.30 

(reyns) 


(29) 


The total torque can now be calculated from the expression 



(30) 


The predictions of this equation are in fairly good agreement with the experi- 
mental results over the entire operating range of the seal. 
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VII. NO -LEAKAGE PRESSURE RISE 

Tests were also performed with no initially applied internal pressure. It 
was found that the internal pressure increased with time due to inward pumping 
of the fluid. This pressure eventually reached an equilibrium value when the 
net leakage-rate vanished. This effect is obviously a manifestation of the 

•A // 

so-called Weissenberg effect. 

Figure 3 shows the result of these tests where AP tot t ^ le measured 
pressure corrected for centrifugal effects. The pressure difference is plotted 
in the figure against the parameter 




because the data indicates a size effect. 

It should be recalled that this pressure difference is related to the 
normal stresses through the expression 







dr 


(31) 


When the right-hand-side of this equation is evaluated for a laminar Newtonian 
flow it is found to be dependent on the shear-rate, as well as the clearance. 

This normal stress term for Newtonian flow is extremely small and ordinarily 
neglected. The important fact is that even in such a simple case this term is 
not just a function of the shear rate. 

It should also be noted that it was assumed earlier that PQ8“Prr was 
independent of the coordinate z, and following that assumption three integrations 
with respect to z where performed on P0Q-P r r* The effect of performing these 
integrations is to introduce' the clearance into the integrated result. The 
way in which the normal stresses depend on the local shear rate determines 
how the clearance enters the integrated expression given above. 

If it is assumed that P^-Pqq is proportional to shear stress square as 
some experimenters suggest, then 


P 68 P rr 


dr 


d'. |a» 

c 


0.8 


(32) 


46 


Since 

r 0 ) 
c 

is only an approximation for the tangential shear rate, a small error in it 
could account for the experimentally observed ~\fc~ factor. 

Correlating the result of Figure 3 yields 



dr = 1.582 (IQ- 3 ) 



p . 

where the constant absorbs the effects of the radius R, . 


( 33 ) 


VIII. LEAKAGE AT LOW RPM WITH AN APPLIED PRESSURE DIFFERENTIAL 

Having examined the performance of the fluid under the special conditions 
of no-leakage and no-rotation, we are now in a position to compare these 
results with the actual performance of the fluid under more general operating 
conditions. 

Since such an excellent correlation for the non-rotating performance was 
obtained, it seems logical to first investigate the perturbing effects of 
modest rotational speeds. Because the rotational speeds are small, the 
centrifugal pumping effect will be negligible and any changes in leakage rate 
can be attributed to changes in fluid properties. 

Of the quantities so far investigated, probably the term containing 
P rr -Pe0 is the one in greatest doubt and least understood from a fundamental 
point of view. On the other hand the term containing 

/ru) 

T z e/r 

which is equal to the viscosity, is by comparison a much better understood 
quantity and, therefore, on a much firmer foundation. This suggests that the 
best way to examine the experimental data for general operating conditions 
is to compute the quantity 


f?. 


-P„ + ir (XU 2 I R 2 p? I - I - 


P r P 2 + 20 p* 


2 *1 


3 J r 

TTC 


dr 


(34) 
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which is equal to the apparent value of 


J F ee- p r r 


R, 


dr 


( 35 ) 


and compare the result with the values of the latter quantity obtained for the 
special case of Q = o. In the first expression above (33), is given by the 
invariant power law derived earlier ^equation (28) j . 

The results of such computations clearly indicate that when Q/rc 2 > »u/c 
the first expression (33) is very nearly equal to zero (i.e., within probable 
experimental error), so that the Weiss enb erg effect or contribution of 


R 


z 



p ■ -p 
96 rr 

r 


dr 


is negligible or overwhelmed by the viscous effects. Typically, this occurs 
at the higher pressure differentials (Pq-P 2 ^20 to 30 psi) which, of course, 
correspond to high leakage rates. Actually, for tests made, m) /c was always 
less than lO^sec"* when rotf/c < Q/rc 2 . According to the no-leakage test results, 
negligible Weissenberg effects would be expected under these circumstances. 

The only disquieting aspect of the calculated results was the fact that 
the small residual was always negative indicating that p- 0 may be a little too 
large by a factor of about 1.2. 


IX. LEAKAGE AT HIGH RPM WITH AN APPLIED PRESSURE DIFFERENTIAL 

When the data is run through the equation above for high ru)/c where rw/c 
is about ten times Q/cr 2 , the apparent value of 


ft*. 



R. 


exhibits several interesting features. 

The first is the rapidity with which the viscous term becomes insignificant. 
At total shear rates of the order of lO^sec - -*-, where the total shear rate is 
given by 
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the pressure and inertial terms entirely overwhelm the viscous term so that the 
apparent value of 







is positive and practically independent of Q. This condition is almost the 
same as the no-leakage case studied earlier and warrants comparison with it. 

Such a comparison reveals that the apparent values are always higher than those 
measured for the no-leakage case, sometimes almost twice as much for the same 
shear rate without leakage. The no-leakage form of the equation 

J] * 99 / 17 -) dr = 1 * 582 ao“ 3 ) -f VT 

R • 

is, therefore, a conservative estimate of its value for seal operation when 
Q * o. 

X. LEAKAGE at intermediate rp m with an applied pressure differen tial 

This regime of operation is the most interesting and also the most 
difficult, since it cannot be compared directly with either of the limiting 
cases previously examined. 

The first important feature of note in this regime is the relative weakness 
of the viscous effects. While more influential than in the limiting case of 
extremely high shear rate, the viscous term is still quite small compared to 
the pressure and inertia effects. The apparent value of the normal stress term 
amounts to a variation around the applied pressure differential with slight 
corrections for inertial effects. The result is that in the total shear rate 
range 2 (10^) to (10^) the quantity ** s proportional to 



for a given clearance and pressure difference. 
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In this range the leakage contribution to the total shear-rate is quite 
small so that to a high degree of accuracy 


P -P 
*99 rr 


a 



which is clearly quite different from the no-leakage results. 

Because the pressure term dominates the viscous term in equation 33 the 
normal stress difference is also proportional to the applied pressure difference 
■ p l" p 2 so t ^ iat f° r a given clearance 



There seems to be no simple expression which will correlate the effects of 
clearance. As noted, the proportion above holds only for a given nominal 
clearance. When the clearance is changed the normal stress difference changes, 
but not in a regular or consistent manner. In fact, there seems to be an 
optimum clearance (about 2 mils) at which the viscoelastic effects are the 
strongest. This is at variance with the no-leakage results where there was a 
definite nonotonic trend with clearance changes. 


XI. CONCLUSION 

It is clear from these findings, as well as those for very high shear-rates, 
that the advent of radial leakage has a significant effect on the performance 
of the viscoelastic fluids on face seals at high shear-rates (of the order 10^ 
or greater). Since these effects can only be inferred indirectly, rather 
than measured directly, it is natural that one would like to assign observed 
unexplainable peculiarities to poor experimental values of the measured 
quantities. 

In particular, it would be nice if the peculiar behavior at moderate rpm 
could be laid at the doorstep of poorly measured viscosity. This is not a 
plausible explanation since it is in this range of shear-rates where the most 
consistent viscosities were obtained. 

Even more surprising is the fact that the best overall correlation (and the 
most believable) was obtained at low shear-rates where the viscosity correlation 
is its poorest. At least in this region the normal stress difference is 
relatively unimportant and, therefore, does not introduce added complications. 
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Considering all the evidence, it is obvious that the normal stress 

difference P n -P is dependent on more than just the tangential shear-rate as 
uU rr 

might be supposed from the examination of viscometer test data where the net 
leakage-rate is zero. As pointed out by several investigators, zero net leakage 
flow fields cannot be the simple flows which are usually pictured, but must 
contain rather complicated secondary flows involving radial and axial velocity 
components. 

From the results obtained here, it is obvious that radial velocity com- 
ponents play a significant role in the behavior of the fluid. The major clue 
in this regard is the normal stress difference dependence on the applied pres- 
sure difference. This pressure difference is probably unimportant in itself^ 
its importance lies in the fact that it is symptomatic of higher leakage-rates 
and high radial shear-rates. 

The key to understanding high shear performance is in determining how 
normal stress difference depends on the radial shear rate. The analysis used 
above avoids a detailed investigation of the velocity field and deals only 
with the gross features of the seal performance, i.e., leakage, applied pres- 
sure difference, etc. Obviously, this is not sufficient when the normal stress 
difference is of the order of the applied pressure difference. Its dependence 
on clearance is evidence that more detailed knowledge is required. 

To carry out a detailed analysis is by no means a simple task. It requires 
first a model constitutive equation, and second, a computer program for solving 
the equations of motion and predicting the overall performance characteristics 
since these equations will not be simple. This is a formidable task, but 
unfortunately there seems to be no alternative. 
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FIGURE 2 

SHEAR-STRESS, SHEAR-RATE CORRELATION 
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MAXIMUM SHEAR RATE Y a , SEC 





NOMENCLATURE 


c 

n 


P 

P 

P 


ee 


= seal clearance 
= invariant power- law index 
= isotropic pressure 
= p-P 

r rr 

= radial component of the deviatoric normal stress 
= tangential component of the deviatoric normal stress 


APfot = 3? l" P 2 + 20 P* ^ R 2 -R P 

Q - volumetric leakage-rate 
r = radial coordinate 

R^ = seal inner radius 

R 2 = seal outer radius 

T = total seal torque 

u = radial component of the velocity 
v = tangential component of the velocity 
z = axial coordinate 

= total shear-rate 

= tangential component of the shear-rate 

pi = viscosity 

|jl o = consistency index 

p = density 

T = tangential component of shear stress 

T = radial component of shear stress 
zr 

u> = angular velocity of seal rotation 
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APPENDIX C 

SEAL GEOMETRY EFFECTS 
BY 

J. A. FINDLAY 


I. INTRODUCTION 

The purpose of this study is to provide an experimental verification of 
the Newtonian face seal theories developed under this contract. These theories 
include the inward pumping theory, the effects of inertia and turbulence on 
dynamic face seals and cavitation effects in the fluid passages which are 
typical of these thin film seals. 

The basic problem in deriving an adequate theory for face seal behavior is 
the lack of accurate data. In general, this is a difficult experimental area 
because at the close clearances required for face seal operation, small, uncon- 
trolled variations in surface geometry play a major role. Thus, the key to 
successful face seal experiments is the controlled, accurate measurement of the 
face seal geometry. These experiments have been designed with this factor as 
an outstanding guide. This allows accurate measurement of face seal performance 
characteristics, such as, leakage, normal force and power loss, as a function of 
known face seal geometry. Thus, the objective of obtaining data adequate for 
precise correlation and understanding of the physics of face seal behavior may 
be attained. Thus, the situation of having what appears to be comprehensive 
seal data, when actually the seal interface geometry is unknown, will thereby 
be avoided. 

The experiments of this task substantiate the theories within the limits 
of the assumptions made. The actual correlation of the experiments with the 
theories follows in this Appendix. 

II. TEST EQUIPMENT AND PROCEDURE 

The rotating and stationary faces of the test seal are shown in Figures 
1, 2, and 3. The metal face of Figure 1 is B6H3 steel which was heat treated 
before grinding was completed. After grinding, the measured runout was * 0.000025 
in. and the flatness was within ~ 0.000005 in. both of which are within 
specifications. Figure 2 shows one of the transparent seal faces mounted on the 
bridge plate which is used to position it. Figure 3(a) shows another transparent 
seal face. This one has a narrower land than the one shown in Figure 2. Both 
of these transparent seal faces are 7740 pyrex glass optical polished plates from 
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Corning Glass. These glass faces are mounted on lucite, a transparent, heat 
resistent acrylic compound. The metal stationary seal face shown in Figure 3(b) 
has the same dimensions as the transparent seal face shown in Figure 2. The 
pressure taps can be seen in the metal face. 

A cross-section of the face seal experiment is shown in Figure 4* In this 
figure, the center feed pipe and film thickness probes are shown. The Sheffield 
displacement sensors were used to set the initial film thicknesses. Three 
sensors were used and from these three readings any combination of average film 
thickness, h^ 9 and tilt eccentricity, e^, could be set. It was found, however, 
that as the supply pressure to the seal was increased the film thickness opened 
up. A static calibration was made using the Sheffield's which showed that the 
film thickness increased uniformly with pressure. Therefore, if only one sensor 
could be read the change in h^ and e could be determined by assuming the same 
increase in film thickness at the other two sensor locations. This was very 
important since the Sheffield's could only be used for static measurements and 
during the dynamic tests the Wayne-Kerr capacitance probe had to be used. 

Because Wa yne-Kerr probes were useful only when a film of oil could be 
maintained between them and the rotating seal face, only one of the probes could 
be used. If three were used they would wipe the oil off the rotating, face and 
no readings of any value could be obtained. 

The Wayne-Kerr probe was calibrated using the test fluid, 350 centistoke 
silicone oil, and the measurements taken with it during the dynamic tests were 
used to correct h^ and e t , This was accomplished using a data reduction program 

Two other corrections to the raw data were made using the data reduction 
program. They were: (1) a pressure correction due to the hydrostatic pressure 
head of oil in the graduated supply reservoir (see Figures 5 and 6) and (2) a 
viscosity correction due to temperature changes measured by a thermocouple in 
the leakage oil. The leakage rate was also calculated using the data reduction 
program using as input the initial and final oil levels in the graduated reservoir 
and the time interval. 

This data reduction program is referred to as No. 1 and its listing is 
given at the end of this Appendix. Data reduction program No. 2 is used for 
calculating the pumping component of leakage in the eccentric seal tests. Its 
listing is also given. Data reduction program No. 3 is similar to No. 1, but 
has the inertia component of leakage added to it. 

These three data reduction programs take the raw data from the tests and 
reduce it down to a form suitable for comparison with the theoretical curves 
shown in Figures 10 and 11 of Appendix D. 
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Photographs were taken for all the test conditions where gas cavities 
appeared. Figures 9, 10, 11 and 12 are samples of some of these figures. Fig- 
ures 9. and 10 are for the concentric seal showing cavities of different sizes. 

v> 

Figures 11 and 12 are for the eccentric seal. Both of these latter figures 
show air being pumped into the center of the seal against a moderate pressure 
differential (about 1 psi). In Figure 11 a small spiral of bubbles may be seen. 
In Figure 12 much more air is being pumped in and it appears as the blurred 
black regions in the center of the seal. 

Figure 5 shows a photograph of the transparent test seal set up. The glass 
face seal is clearly shown as are the Sheffield displacement probes and the 
Wayne-Kerr capacitance probe. The tube connected to the center of the seal is, 
of course, the oil supply line. A more complete view of the oil supply line, 
graduated reservoir and supply pressure system may be seen in Figure 6, This 
figure also shows the movie camera used to take motion pictures of the gas 
cavities and the inward pumping of air. 

III. CAVITY EXPERIMENTS ( CONCENTRIC «, MISALIGNED SEAL) 

A. INTRODUCTION 

The first seal theory to be studied experimentally was that of gas cavities 
in the fluid film between the seal faces. The short bearing theory developed 
during Phase I of this contract indicated that these gas cavities, under 
certain conditions would influence the hydrodynamic forces but not the leakage. 
This was true at least for the concentric seal. That is, the leakage would be 
the same as for the hydrostatic seal. However, the numerical solutions obtained 
on the GE 635 Computer Program gave leakage rates that increased with A and, 
therefore, cavity size. It was suspected at the time that this was the result 
of inadequately describing the boundary shape of the cavity in the GE 635 
Program (Reference 1) . 

B. TEST PROCEDURES 

The experimental program schedule was set up to settle the difference 
between the two methods of analyses described above. Tests were performed on 
the transparent test seal described in Part II. 

In these experiments the oil was fed to the inside of the seal under 
pressures from 1 to 40 psig and the outside diameter of the seal was exposed to 
air. The stationary seal member was held at a fixed tilt to the rotating 
member so that a "converging-diverging" film shape was, obtained between the seal 
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faces with a film thickness variation of about two-to-one. For these conditions, 
air was "sucked" into the seal interface in the diverging film region forming 
a gas cavity at ambient pressure. A range of pressures and speeds were used so 
that relative magnitude of the hydrodynamic effects could be varied. Therefore, 
leakage measurements were obtained for conditions ranging from the hydrostatic 
up to A = 1000. A summary of the concentric seal experiments is given in Table 
I. It will be noted that only a few tests were made with the wide land seal. 

This seal was cracked during one of the tests and the remainder of the experi- 
ments were carried out using the narrow land seal. The wide land seal is shown 
in Figure 2 and the narrow land seal in Figure 3(a). 

C. RESULTS 

1. Flow Rates 

a. Wide Land Seal 

The data taken witfr the side land seal are given in Table II and shown cor- 
related in Figures 14 and 15. The correlation is not as good as hoped for but 
this is attributed to the fact that no correction was made for the hydrostatic 
pressure head in the supply line and the difficulty experienced in obtaining 
close control on the film thickness measurements. 

The maximum error in this correlation of the flow coefficient was about 
100$. This could be due to an error of only 30$ in the film thickness. Most 
of the data correlated within about 25$. This would represent a film thickness 
error of only 8$. 

The important thing to note is that the rotational speed had no measurable 
effect on the leakage. Since the gas cavities, which were photographed, in- 
creased in size with speed it can be concluded that the presence and size of 
the gas cavities had no effect on the seal leakage. 

Observing a set of points, in Figure 14 and 15, for a constant speed it can 
be noted that the pressure level did effect the accuracy of the correlation. 
Since the pressure level influences the film thickness, this further indicates 
that the accuracy of the film thicknesses may be the weak link in this 
correlation. 

b. Narrow Land Seal 

The correlation of data with theory is much better in these experiments. 
Again the scatter of the data is attributed mainly to the accuracy of the film 
thickness measurements. Also the conclusion that the cavity does not affect 
the leakage, is again substantiated. 
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The data for these experiments are summarized in Table III and Figures 16, 
17, 18, 19 and 20. 

2. Cavity Shapes and Inception 

One of the objectives of these experiments was to see if the observed 
cavity shapes could be predicted, and if the conditions under which the cavity 
first appeared could be predicted. These correlations are not yet complete, 
but some general observations can be made at this time. 

First of all, for a given geometry the theory says that the cavity will 
first appear for a given, fixed value of A. That is, as the pressure across 
the seal is increased so must the rotational speed be increased. This observa- 
tion has been confirmed. What is more difficult to confirm is the exact value 
of A for which the cavity first appears. This is due to the fact that the 
cavity appears gradually and an observation cannot be made until a finite 
cavity appears. 

Secondly, the cavity shapes are generally the same as predicted in the 
Phase I study for the wavy seal. However, cavity shapes have not yet been 
calculated for the misaligned seal although the two cases are similar in 
geometry. The theory assumes that there are many small ribbons of oil flowing 
through the cavity. However, with the relatively large film thicknesses used 
in these experiments there are only a few large ribbons of oil. This will 
affect the correlation of the theory and experiments but not drastically. 

A summary of cavity shapes and corresponding seal conditions is given in 
Figure 8 and Table XI. Host of data taken for the concentric seal leakage tests 
had photographs. 

Many of these, while showing the gas cavities, are not clear enough for 
publication. However, some measurements were taken from these photographs. As 
noted before, the value of A necessary for a cavity to first form (inception) 
was very difficult to measure. A correlation between 9 a (cavity nose) and A 
has not been pursued to any extent. However, a few measurements were made of 
J^crit vs A. A correlation of these measurements against theory is shown in 
Figure 34. So far this correlation, also, is rather crude. The observations 
made do indicate that the cavity shapes are given to a fair approximation by 
the theory. This correlation could be extended in the future. 

The theoretical curve shown in Figure 34 was obtained from the Phase I 
cavity theory for the wavy seal and is, therefore, an approximation to this 
case. The data points were taken from experiments on a misaligned or tilted 
seal. 
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IV. SEAL PIMPING EXPERIMENTS (ECCENTRIC. MISALIGNED SEAL)., 

A. INTRODUCTION 

In Appendix D, an analysis is developed for the wavy, radially eccentric 

i 

face seal showing the existence of a "pumping" leakage. This pumping leakage 
occurs because of the combined effects of a once per rev film thickness 
variation, due either to waviness or misalignment, and a radial displacement of 
the seal center from the center of rotation. 

The assumption made in this analysis was that a complete fluid film exists 
between the seal surfaces. Since in the misaligned seal experiments gas cavities 
usually occurred, the pumping experiments were broken down into three groups: 

(1) inward pumping, with cavities, (2) inward pumping, no cavities, and (3) out- 

/ 

ward pumping, no cavities. A fourth category, "full film with different fluids 
on either side of the seal," is discussed in Part V. 

B. TEST PROCEDURES 

In the group (1) experiments the seal radial eccentricity was in such a 
direction, in relation to the seal misalignment, so that the hydrodynamic pump- 
ing action was inward, or- opposed to the hydrostatic leakage. Also in these 
experiments the outer circumference of the seal was exposed to air. This re- 
sulted in the air being drawn into the seal film where the surfaces diverged, 
thus forming gas cavities. The gas cavities reduce the generation of hydro- 
dynamic pressures in the seal film and thus reduce the pumping action that 
occurs because of the seal eccentricity. For the high A tests some air was 
pumped through the seal to the inside against the pressure. 

The inward pumping experiments with no cavities (group 2) were accomplished 
by building up an oil dam on the rotating seal face outside of the stationary 
seal face; and supplying oil between the dam and the seal. Thus, the outside 
circumference of the seal was flooded with oil, and no gas cavities were formed 
in the seal film during these experiments. The pumping action in the seal was 
exactly as predicted and in a direction to oppose the hydrostatic leakage. 

Under high A conditions a net inward leakage was obtained against the hydrostatic 
pressure head of the oil in the graduated reservoir. The experimental set-up 
for this experiment is shown in Figure 13. This figure shows the rubber tubing 
and bulb used to supply the oil to the outside of the seal. This oil was dyed 
blue and as it was pumped into the center of the seal it filled the normal 
supply line and flowed up into the oil reservoir. This dark colored oil can be 
seen in Figure 13 also. 
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The group (3) experiments were performed with the radial eccentricity set 
in a position 180° from that in the group (1) and (2) experiments. This changed 
the cos a term in the pumping equation from -1 to +1. In this case, the pump- 
ing action was outward, in the same direction as the hydrostatic leakage. No 
gas cavities formed during these experiments and the data correlated very well 
with the theoretical predictions. The leakage in these experiments was greater 
than the hydrostatic leakage and increased linearly with A. 

C. RESULTS 

From the preceeding discussion it should be clear that the inward and out- 
ward experiments with no cavities should correlate with the theory developed in 
Appendix D„ This correlation is shown in Figures 21, 22 and 23. Figure 21 
shows the q pump/e t data from the outward pumping experiments and its comparison 
with the identical curve for e r = 0.091. As can be seen, the correlation is 
very good. 

Figure 22 shows the leakage coefficient data from the inward pumping, with- 
out cavities, experiment. Again the comparison with the theoretical curve is 
very good. 

Both of the preceeding sets of data are plotted in Figure 23 using the 
invariant parameter q pump/Ae^t cos a. For the seal tested, = 0.83, the 
theoretical value of this parameter is 2.63 as obtained from Appendix D. The 
correlation on this basis does not look as good but all the errors are magnified 
in this type of plot. In addition, it is difficult to obtain accurate pumping 
flows since they must be determined by subtracting the measured leakage from the 
theoretical concentric seal leakage. Since the concentric seal correlations 
showed some scatter, it is obvious that the accuracy of the pumping flow mea- 
surements must be affected even more since a difference in flows is involved. 

Figures 24 through 29 show the inward pumping data for seals containing gas 
cavities. The first conclusion that can be drawn from these curves is that the 
presence of cavities decreases the inward pumping effect. 

In Figures 24 and 25 the same data is plotted in two different ways. From 
these curves, it can be seen that at the low values of A, where no cavities exist 
in the fluid film, the data compares quite closely to the noncavity seal pumping 
theory. However, at the higher values of A, where gas cavities do form, the 
effect of seal pumping decreases. 

Figures 26 through 29 also show the decrease in seal pumping with increased 
A. There is more scatter in this data, however. It is felt that the reason for 
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this is that the data in these figures is for smaller film thicknesses, and 
the accuracy of the measurements and, therefore, the leakage correlation is not 
as good as for the larger film thickness tests. Even the static leakage cor- 
relation for h^ = 0.001" was not sufficiently accurate to allow meaningful values 
of the pumping coefficient to be determined. However, in Figure 27, where some 
of this data is plotted, still shows the general trend of decreasing leakage 
with increasing A due to the seal eccentricity, i.e., inward pumping. 

Figure 29 illustrates the data for e fc = 0.20. At this low value of tilt 
eccentricity the pumping effect is small and, therefore, the accuracy in mea- 
suring it is relatively poor. Thus the data points in Figure 29 show quite a 
bit of scatter. 

The low pressure, inward pumping tests were repeated with a slightly 
different feed system. This was done because some of the low pressure inertia 
tests showed a great deal of scatter. In looking into the problem it was 
found that air was trapped in the supply tube and thus giving a lower hydro- 
static pressure at the inside of the seal than assumed. i 

The feed system was modified so that the reservoir stand pipe was placed 
directly over the feed hole in the center of the stationary face seal. This 
reservoir was then connected to the feed hole by a piece of transparent plastic 
tubing. All air was thereby eliminated from the feed system and the scaling 
pressure could be accurately determined by the height of the oil in the reservoir. 

This set up was used to: (1) repeat the low pressure, inward pumping 
experiments with no cavities, (2) repeat the low pressure inward pumping 
experiments with cavities, (3) perform the no cavity inward pumping tests with 
two different fluids on either side of the face seal, and (4) perform the 
inertia experiments. 

Figure 30 shows the results of the low pressure, inward pumping experiments, 
with and without cavities. The experiments without cavities correlates very 
well with the theory developed in Appendix D. Both of these experiments were 
carried out with silicone oils of two different viscosities, i.e., 350 and 100 
centistokes. The correlations between the two fluids is excellent. 

It will again be noted that the presence of gas cavities in the film re- 
duces the inward pumping action of the eccentric seal. At low values of A the 
two curves are identical. However, as A increases and a cavity starts to form, 
in the oil-air experiments, the two curves diverge and the inward pumping 
component of leakage reaches a limiting value. This is demonstrated by the fact 
that the net leakage reaches a steady value as A increases. This can also be 
seen in Figure 32 for values of A up to 4100. 
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The analytical relationships for the "inward pumping with cavity" case has 
been set up in Appendix D, but no solutions have as yet been obtained. Part of 
this has been set up on the computer, but the final programming was not 
accomplished in the time allowed. It is recommended that this work be continued. 

V. TOO FLUID a PUMPING EXPERIMENTS 

A* INTRODUCTION 

Up to this point, the very important phenomenon of seal pumping has been 
investigated for the single fluid cases. Another very important application is 
where the face seal separates two different fluids. 

If a hydrodynamic film exists between the two seal surfaces, then the two 
fluids will meet in the seal interface and the seal hydrodynamics will be 
modified. 

This problem has been set up in Appendix D, but again, time did not permit 

the completion of this solution. However, experiments have been run and the 

data correlated with the viscosity ratio. As suspected, if the viscosity of 
the fluid on the outside of the seal is lower than the viscosity of the fluid on 
the inside of the seal, then the pumping action is reduced. 

B. TEST PROCEDURE 

The experimental apparatus was set up with stand pipe reservoir directly 
over the feed hole in the center of the stationary face seal, as explained in 

IV. Silicone oil of one viscosity was supplied to the center of the seal under 

hydrostatic pressure by means of this stand pipe. A fluid dam was mounted on 
the rotating seal face, the same as for the inward pumping - no cavity experi- 
ments. This outside reservoir, at ambient pressure, was flooded with a silicone 
oil of a different viscosity. 

Two things modified the procedure somewhat. The first was that at high 
speeds (high A) cavities started to form because not enough oil could be supplied 
to the outer reservoir. Therefore, the maximum value of A for these experiments 
was limited. The second was that the net leakage of the fluids through the 
seal would "contaminate" the other fluid. That is, at low A the inside fluid 
leaked to the outside reservoir and mixed with the other oil there. Thus, the 
viscosity of the oil on the outside of the seal tended to be variable. At high 
A the outer fluid was pumped into the center of the seal, thus mixing with the 
fluid there and modifying its viscosity. 
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Therefore, periodically the experiments were stopped and the inner and 
outer reservoirs flushed before continuing the measurements. 

G. RESULTS 

The measurements taken in these experiments are correlated in Figures 31 
and 32. Viscosity ratios, p. 0 /p«£» ranging from 3.5 to 0.002 were used. Although 
for the single fluid case, i. - 1, the data showed very little scatter, and 

lay on a straight line; the two fluid data scattered and it is not certain that 
the correlations should be straight lines or not. Theoretical solutions have not 
yet been obtained for this case. Because of the fluid mixing, as discussed 
above, the two fluid data scattered and repeating the same conditions would 
yield slightly different results depending on whether or not the reservoirs had 
just been flushed or whether the net leakage had been inward or outward. Since 
it takes a finite length of time to obtain a single data point, and fluid mixing 
starts immediately, it was not possible to eliminate the scatter even when the 
reservoirs were flushed frequently. 

VI. INERTIA EXPERIMENTS 

A. INTRODUCTION AND TEST PROCEDURE 

In the Newtonian study performed under Phase I by H. J. Sneck, laminar 
flow with inertia effects was considered. The results of this analysis showed 
that the leakage solution was unchanged except for the addition of a centrifugal 
inertia term. Thus, the concentric, misaligned seal leakage coefficient is 



The experiments performed up to this point have been at low speed and with 
high viscosity oils. Thus, large hydrodynamic pressures could be generated while 
the inertia terms were minimized. In order to obtain significant inertia effects, 
it was necessary to run the seal experiments at high speeds. To do this, a low 
viscosity oil (5-10 centistokes) had to be used in order to reduce the friction 
and, therefore, the torque required and heat generated. A compromise between 
oil viscosity and average film thickness had to be arrived at so that the leakage 
rates would not be too great nor the film thickness error too large. Values of 
10 centistokes and 0.003 inches were used. Both aligned tests and a tilt 
eccentricity, ® t , of 0.56 were used. 
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B. RESULTS 

The experimental measurements are correlated with the theoretical values 
in Figure 33 „ The leakage equation was rearranged as follows 


-q lnR i /2n(l + a e^> = 1 
The inertia term 


3_ 

20 





assumes that all the inertia effects take place within the seal interface. 
However, it was quite obvious that because of the rotating plate, rather than a 
rotating ring, the inertia effects were felt throughout the seal. 

The final correlations showed that if an inertia term of 
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was used, the correlation was quite good, This can be seen in Figure 33. 


VIo CONCLUSIONS 

Although some scatter was found in certain of the data, there was sufficient 
consistency in most of the data to confirm the misaligned seal leakage theory, 
with and without inertia effects, and the eccentric seal pumping theory. 

Since the eccentric, misaligned seal theories for pumping with cavities 
and with two fluids have not been completed, no verification of these theories 
have been made. However, sufficient data has been obtained to evaluate the 
effect and give a basis for evaluating the theory. 

These theories should be completed and incorporated in the digital computer 
SEAL PROGRAM. 
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Figure 1. Rotating Seal Face 
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Figure 2. Transparent Stationary Seal Face Mounted on the 
Bridge Plate. (Sheffield displacement sensors 
are shown. ) I. D. = 2. 25 in. - O. D. = 3. 50 in. 
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(a) 

Transparent Stationary 
Seal Face 
I. D. = 2. 25 in. 

O. D. = 2. 75 in. 


(b) 

Steel Stationary 
Seal Face 
I. D. = 2. 25 in. 
O. D. = 3. 50 in. 


Figure 3 
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Figure 4. Newtonian Face Seal Experiment 
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Figure 5. Narrow Land, Transparent Face Seal Setup 




Figure 6. Face Seal Experiment Showing Movie Camera 
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(c) (d) 

Figure 7. Gas Cavities in a Transparent Face Seal. (The dark areas are 
the gas cavities and the light areas oil. Rotation is in the 
clockwise direction with the minimum film thickness, of the 
misaligned seal faces, coinciding with the full film region) 


74 










(c) (d) 

Figure 8. Gas Cavities in the Narrow Land Transparent Face Seal, H- = 0. 83. 

(The dark areas are the gas cavities and the light areas oil. Rotation 
is in the clockwise direction with the minimum film thickness, of the 
misaligned seal faces, at the left of the photos. ) See Table XI. 
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Figure 8 (Continued) 
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(o) (p) 

Figure 8 (Continued) 
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Figure 10. Concentric, Transparent Face Seal Showing Gas Cavity. 
(Smaller cavity, lower speed) 
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Figure 11. Eccentric, Transparent Face Seal Showing Gas Cavity and Inward 
Pumping of Air. (Note spiral of bubbles in center) 
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Figure 12. 


Eccentric, Transparent Face Seal Showing Gas Cavity. 
(Note many bubbles in center due to higher speed) 
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Figure 13. Face Seal Experiment Setup for Flooded Seal 
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FIGURE 1.4 
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FIGURE 15 

MISALIGNED SEAL LEAKAGE 
H = 1 +R€ t COS0 
R| =0.665 
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FIGURE 17 

MISALIGNED SEAL LEAKAGE 
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FIGURE 19 

MISALIGNED SEAL LEAKAGE 
H s I + R € t COS 0 
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LEAKAGE COEFFICIENT 


FIGURE 20 

MISALIGNED SEAL LEAKAGE 
H =1 + R6 t COS0 




pump 


FIGURE 21 
OUTWARO PUMPING 
DUE TO 

RADIAL ECCENTRICITY 
(NO CAVITY) 
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FIGURE 22 

LEAKAGE COEFFICIENT 
SHOWING THE EFFECTS 
OF INWARD PUMPING 
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FIGURE 23 



CORRELATION OF INWARO 
AND OUTWARD PUMPING 
EFFECTS DUE TO 
SEAL ECCENTRICITY 
(NO CAVITY) 
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THEORY 2.63 
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FIGURE 25 

PUMPING COEFFICIENT 
FOR 

ECCENTRIC FACE SEAL 
WITH CAVITIES 
R 1=0.83, € r =0.091 
hj =0.004", € t =0.50 
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FIGURE 26 

PUMPING COEFFICIENT 
FOR 

ECCENTRIC FACE SEAL 
WITH CAVITIES 
Rj = 0.83, 6 r =0.091 
hj =0.003", €f =0.50 
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FIGURE 27 

LEAKAGE COEFFICIENT 
FOR 

ECCENTRIC FACE SEAL 
WITH CAVITIES 
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FIGURE 28 

PUMPING COEFFICIENT 
FOR 

ECCENTRIC FACE SEAL 
WITH CAVITIES 
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FIGURE 29 

PUMPING COEFFICIENT 
FOR 

ECCENTRIC FACE SEAL 
WITH CAVITIES 
Rj 8 0.83, € r 8 0.091 
hj 8 0.003", €f 8 0.20 
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FIGURE 31 

INWARD PUMPING EXPERIMENTS 
TWO FLUIDS 
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FIGURE 32 

INWARD PUMPING EXPERIMENTS 
TWO FLUIDS 
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FIGURE 33 

INERTIA EXPERIMENTS 
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FIGURE 34 

CORRELATION OF CAVITY SHAPES 
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TABLE I 


SUMMARY OF CONCENTRIC SEAL TESTS 


Table No. 

h i9 in 

e t 

Remarks 

Figures 

II - 1 

.004 

.50 

Wide land seal 

14 

II - 2 

I! 

.60 

it ii it 

15 

II - 3 

II 

.70, .80 

, ii ii ii 

15 

III - 1 

II 

.50 

Narrow land seal 

16 

III - 2 

11 

.20 

11 It It 

17 

III - 3 

II 

.80 

II II II 

18 

III - 4 

• 003 

.50 

11 11 11 

19 

III - 5 

• 003 

.50 

II II II 

20 


I 

I 
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TABLE II 


Wide Land Seal Data - Concentric 
d Q = 3.469 inch; dj = 2. 312 inch; Rj = 0. 665 inch 

1. This set of data had no dynamic film thickness measurements or 
viscosity temperature correction (Figure 14 ) „ 


p i, ; 

N 

h i 

Q 


psi 

rpm 

iHi 

in 3 / sec 

q 

5 

0 

0.00417 0.581 

0.0176 

28.9 

10 

0 

.00430 .560 

. 0351 

26.6 

15 

0 

.00436 .555 

.0540 

25.7 

20 

0 

.00448 .535 

. 0766 

25.2 

30 

0 

.00468 .510 

.1225 

23. 5 

40 

0 

.00492 .485 

. 1905 

23.8 

5 

200 


0.0181 

29.7 

10 

200 


.0361 

27.4 

15 

200 

(Assumed to be the 

.0619 

29.4 

20 

200 

same as above) 

.0746 

24. 6 

30 

200 


. 1190 

22.8 

40 

200 


.1856 

23.2 

5 

400 


0.0189 

31.2 

10 

400 


. 0373 

28. 4 

15 

400 

(Assumed to be the 

. 0534 

25.4 

20 

400 

same as above) 

.0723 

23. 8 

30 

400 


. 1238 

23. 6 

40 

400 


. 1810 

22, 6 

5 

600 


0. 0156 

25. 6 

10 

600 


. 033 1 

25.2 

15 

600 

(Assumed to be the 

. 0495 

23.6 

20 

600 

same as above) 

. 0702 

23.2 

30 

600 


. 1145 

22.0 

40 

600 


. 1679 

21. 0 

5 

600 


0.0146 

24. 0 

10 

800 


. 0273 

20. 8 

15 

800 

(Assumed to be the 

. 0427 

20.4 

20 

800 

same as above) 

. 0672 

22. 2 

30 

800 


.0970 

18. 6 

40 

800 


. 1881 

23. 5 
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Table II (continued) 


2. 

This set of data includes a temperature correction on viscosity 



and dynamic film thickness connections (Figure 15 ). 



Pi 

N 



Q 


psi 

rpm 

in. 

et 

in 3 / sec 

q 

5 

400 

0.004107 

0.595 

0.0193 

32.8 

10 

400 

.00415 

. 590 , 

. 0378 

29.6 

15 

400 

. 00425 

.575 

.0576 

29.0 

20 

400 

.004375 

.558 

.0794 

25.9 

30 

400 

.00455 

.538 

. 129 

23.8 

5 

600 

0.00416 

0.587 

0.0230 

35.4 

10 

600 

.004147 

.590 

.0418 

30. 1 

15 

600 

. 00416 

. 587 

. 0605 

29. 7 

20 

600 

.004188 

.585 

.0794 

27.7 

30 

600 

.004336 

.565 

. 130 

26.9 

40 

600 

. 004472 

.545 

.206 

29. 2 

5 

800 

0. 003812 

0.640 

0.0227 

42.0 

10 

800 

.003853 

.635 

.0408 

39.2 

15 

800 

. 003853 

.635 

.0518 

30.4 

20 

800 

.003893 

.627 

.0734 

31. 1 

30 

800 

. 004027 

.605 

.1130 

28.6 

40 

800 

.004188 

. 585 

. 1630 

28.3 

3. 

hj = 0. 004; e t = 0.70; 

0.80 (Figure 15 ) 

• 



5 

400 

0.004153 

0.760 

0.0247 

43.2 

10 

400 

. 004153 

,760 

.0471 

38.9 

15 

400 

.004181 

.756 

.0640 

33. 6 

20 

400 

. 004234 

.746 

.0878 

33.3 

30 

400 

. 004368 

.723 

. 1280 

29.6 

40 

400 

.004516 

.700 

. 2050 

32. 5 

5 

600 

0. 003929 

0.804 

0.0224 

43. 5 

10 

600 

. 003929 

.804 

. 04 14 

39. 1 

15 

600 

.003993 

. 791 

. 0624 

33. 3 

20 

600 

. 004033 

.785 

. 0827 

32. 6 

30 

600 

. 004181 

. 756 

. 1320 

33. 8 

40 

600 

.004315 

.732 

. 1890 

33. 0 

5 

800 

0.004261 

0.741 

0. 0238 

32.4 

10 

800 

.004153 

. 760 

. 0517 

39.2 

15 

800 

. 004153 

.760 

.0750 

34.6 

20 

800 

. 004181 

.756 

. 0974 

33.4 

30 

800 

.004261 

.741 

. 1490 

32.4 

40 

800 

. 004422 

. 715 

. 2150 

31. 2 
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TABLE III* 


Narrow Land Seal Data - Concentric 
d Q = 2. 750 inch; dj = 2/281 inch; - 0.83 inch 


1. = 0.004 inch, e^. = 0. 50 (Figure 16 ). 


Pi 

N 

h. 

l 


Q 


psi 

rpm 

in. 

e t 

in 3 /sec 

q 

2 

0 

0.00400 

0.492 

0.0138 

46.8 

4 

0 

.00403 

.490 

.0206 

39. 5 

6 

0 

. 00406 

.484 

.0301 

39. 9 

8 

0 

, 00401 

.480 

.0424 

42.2 

10 

0 

.00413 

.478 

.0513 

40. 6 

15 

0 

.00421 

.469 

.0798 

40. 9 

20 

0 

.00428 

.459 

. 1117 

41.2 

30 

0 

. 00443 

.445 

.1750 

39.4 

2 

200 

0.00397 

0.497 

0.0107 

37.4 

4 

200 

. 00402 

.490 

. 0197 

38. 5 

6 

200 

.00407 

.484 

.0311 

40.8 

8 

200 

.00409 

.482 

.0410 

41,2 

10 

200 

.00410 

.480 

.0521 

41.9 

15 

200 

.00418 

.471 

.0787 

40.7 

2 

400 

0.00402 

0.490 

0.0073 

23.4 

4 

400 

. 00403 

.489 

.0195 

36. 1 

6 ■ 

400 

.00403 

.489 

.0301 

39.3 

8 

400 

.00405 

.487 

. 0380 

37.9 

io ; 

400 

. 00407 

.484 

.0515 

41. 1 

15 

400 

.004 15 

.474 

.0756 

38.5 

2 : 

600 

0. 00395 

0.499 

0.0090 

29.6 

4 

600 

. 00394 

. 500 

,0193 

37.4 

6 

6Q0 

.00394 

. 500 

.0254 

34.4 

8 

600 

.00397 

.497 

.0365 

37. 5 

10 

600 

.00398 

. 495 

.0461 

38.2 

15 

600 

. 00405 

.487 

.0713 

37. 6 

2 

800 

0.00378 

0.522 

0. 0059 

21. 7 

4 

800 

.00378 

.522 

.0151 

31.6 

6 

800 

.00378 

.522 

. 0230 

33. 1 

8 

800 

.00378 

. 522 

.0344 

38.2 

10 

800 

. 00378 

.522 

.0396 

35.4 

15 

800 

.00381 

.518 

. 0586 

34.4 

*NOTE: 

These data were reduced using data reduction program No. 

1. 
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Table III (continued) 


2. hj = 0.004 inch, e^. = 0.20 (Figure 17) 


Pi 

N 

h i 


Q 


psi 

rpm 

in. 


in 3 / sec 

q 

2 

0 

0. 00405 

0.206 

0.0108 

35.2 

4 

0 

.00409 

.204 

.0163 

30.2 

6 

0 

. 00413 

.204 

. 0266 

33.9 

8 

0 

.00414 

. 199 

.0344 

33.5 

10 

0 

.00417 

.200 

. 0430 

33.4 

15 

0 

.00426 

. 194 

. 0677 

33.7 

20 

0 

.00432 

. 191 

. 0938 

34.0 

30 

0 

.00455 

. 161 

. 1549 

32,5 

2 

200 

0.00433 

0. 169 

0,0084 

23. 0 

4 

200 

. 00437 

.168 

.0150 

23. 1 

6 

200 

.00440 

. 167 

. 0276 

29.5 

8 

200 

. 00443 

. 166 

. 0345 

28.0 

10 

200 

.00444 

. 165 

. 0451 

29.6 

15 

200 

. 00453 

. 162 

. 0690 

29. 0 

2 

400 

0. 00433 

0. 169 

0.0102 

29. 0 

4 

400 

. 00435 

. 169 

.0176 

28.3 

6 

400 

. 00437 

.168 

. 0281 

31.0 

8 

400 

.00440 

. 167 

. 0356 

29.2 

10 

400 

. 00443 

. 166 

.0445 

28.8 

15 

400 

. 00445 

. 165 

. 0661 

28.6 

2 

600 

0.00427 

0. 172 

0. 0069 

20.0 

4 

600 

. 00428 

.171 

.0159 

26. 1 

6 

600 

.00427 

. 172 

. 0261 

30.3 

8 

600 

. 00429 

.171 

. 0305 

26.9 

10 

600 

. 00433 

. 169 

. 0389 

27.2 

15 

600 

.00439 

.167 

. 0625 

28.7 

2 

800 

0.00421 

0. 174 

0. 0075 

23 . 1 

4 

800 

.00416 

. 176 

.0126 

22.4 

6 

800 

. 00416 

. 176 

.0217 

26.8 

8 

800 

. 00416 

.176 

. 0260 

24. 1 

10 

800 

. 00415 

. 177 

. 0395 

29. 5 

15 

800 

. 00421 

. 174 

. 0562 

27. 3 


111 


Table III (continued) 


3. Tij = 0. 004 inch, = 0. 80 (Figure 18 ) 


Pi 

N 

hi 


Q 


IPsi 

rpm 

in. 

®t 

in 3 /sec 

q 

2 

0 

0. 00380 

0.915 

0.0163 

67.8 

4 

0 

. 00385 

.902 

. 0298 

69.3 

6 

0 

.00386 

.901 

. 0419 

68,0 

8 

0 

. 00389 

.892 

. 0394 

48.2 

10 

0 

.00392 

.887 

. 0645 

63. 1 

15 

0 

. 00399 

.870 

. 1029 

64.9 

2 

200 

0. 00401 

0.796 

0. 0126 

43.9 

4 

200 

. 00403 

. 791 

. 0263 

52.4 

6 

200 

. 00406 

. 786 

. 0401 

55.2 

8 

200 

. 00409 

. 780 

. 0497 

52. 0 

10 

200 

. 00410 

.778 

. 0621 

52.4 

15 

200 

. 00417 

. 765 

. 1009 

55.3 

2 

400 

0, 00401 

0.796 

0. 0082 

28.6 

4 

400 

. 00403 

.791 

. 0239 

47. 5 

6 

400 

.00403 

. 791 

. 0405 

56.8 

8 

400 

. 00406 

. 786 

.0510 

53.6 

10 

400 

. 00409 

.780 

. 0601 

50.6 

15 

400 

. 004 14 

. 770 

. 0879 

47. 6 

2 

600 

0.00393 

0.812 

0. 0107 

41.6 

4 

600 

. 00395 

.807 

. 0249 

52. 5 

6 

600 

. 00395 

.807 

. 0339 

50. 0 

8 

600 

. 00398 

.802 

-- 

— 

10 

600 

. 00398 

.802 

. 0629 

55. 1 

15 

600 

.00401 

.796 

.0902 

52.7 

2 

800 

0. 00382 

0.835 

■ 


4 

800 

.00387 

. 824 

. 0221 

48. 2 

6 

800 

. 00385 

.829 

. 0323 

50. 7 

8 

800 

.00385 

.829 

. 0447 

52. 5 

10 

800 

. 00385 

.829 

. 0541 

51.3 

15 

800 

. 00387 

.824 

. 0885 

55. 6 
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Table III (continued) 


4. hj = 0. 003 inch, = 0. 50 (Figure 19 ) 


Pi 

N . 

hi 


Q 


psi 

rpm 

in. 

e t 

in 3 / sec 

^ , 

5 

0 

0.00311 

0. 515 

0. 0134 

47.4 

10 

0 

. 00319 

. 502 

. 0282 

49. 6 

15 

0 

. 00327 

.490 

.0415 

46. 3 

20 

0 

. 00334 

.482 

. 0629 

50. 2 

25 

0 

. 00343 

.467 

. 0830 

49.4 

30 

0 

. 00351 

.461 

. 1031 

47. 7 

5 

200 

0, 00310 

0. 515 

0. 0106 

o 

I> 

CO 

10 

200 

. 00318 

. 502 

. 0276 

47.8 

15 

200 

.00326 

.490 

.0431 

4TA 

20 

200 

. 00335 

.478 

. 0544 

42.2 

25 

200 

.00340 

.471 

. 0804 

48. 0 

30 

200 

* 00345 

.463 

. 1025 

48.9 

5 

400 

0.00310 

0. 515 

0.0121 

43.3 

10 

400 

.* 00318 

. 502 

. 0217 

38.2 

15 

400 

00324 

.494 

. 0419 

46. 2 

20 

400 

.00329 

.486 

. 0605 

48. 2 

25 

400 

. 00335 

.478 

. 0793 

47.9 

30 

400 

. 00343 

.467 

. 0972 

46.3 

5 

600 

0. 00302 

0. 529 

0. 0089 

32.7 

10 

600 

* 00310 

. 515 

. 0212 

38*4 

15 

600 

. 00318 

. 502 

. 0379 

42. 6 

20 

600 

. 00324 

.494 

. 0528 

42.4 

25 

600 

. 00326 

.490 

-. 0732 

46. 3 

30 

600 

. 00332 

.482 

. 0902 

44 . 6 

5 

800 

0. 00294 

0. 544 

0.0155 

59.8 

10 

800 

. 00300 

. 534 

. 0202 

38.8 

15 

800 

. 00302 

. 529 

.0342 

43. 5 

20 

800 

. 00308 

.520 

* 0478 

42.3 

25 

800 

. 00310 

.515 

. 0622 

43.2 

30 

800 

. 00317 

.504 

. 0826 

43. 7 
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Table III < continued) 


5. hj = 0. 003 inch, e t = 0.50 (Figure 20 ) 


Pi 

N 

h i 


Q 


Mi 

5 

rpm 

in. 


in 3 / sec 

q , 

200 

0. 00309 

0.518 

0.0151 

53.5 

10 

200 

.00317 

.504 

. 0303 

53. 1 

15 

200 

.00324 

.494 

. 0450 

50.8 

20 

200 

. 00332 

,482 

. 0662 

52,7 

25 

200 

. 00337 

.474 

. 0879 

53.8 

30 

200 ' 

.00345 

.463 

. 1035 

49.4 

5 

400 

0.00308 

0. 520 

0.0151 

54. 2 

10 

400 

. 00316 

.507 

. 0352 

62. 5 

15 

400 

. 00321 

,498 

.0489 

56. 5 

20 

400 

.00326 

.490 

.0642 

52.5 

25 

400 

. 00332 

.482 

, 0860 

52.9 

30 

400 

. 00337 

.474 

. 1050 

51. 6 

5 

600 

0. 00300 

0.534 

0. 0159 

61.6 

10 

600 

. 00310 

. 515 

.0313 

58.5 

15 

600 

. 00313 

. 511 

. 0469 

58.4 

20 

600 

. 00318 

.502 

, 0586 

51. 6 

25 

600 

. 00324 

.494 

.0816 

54. 0 

30 

600 

. 00326 

.490 

.1020 

55.2 

5 

800 

0. 00289 

0. 554 

0. 0128 

54. 6 

10 

800 

. 00294 

. 544 

. 0286 

61. 0 

15 

800 

. 00297 

.539 

. 0330 

46. 1 

20 

800 

.00300 

. 534 

. 0557 

56.9 

25 

800 

. 00302 

.529 

.0651 

52. 1 

30 

800 

. 00309 

.518 

. 0848 

51. 5 
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A. Listing 


TABLE IV 

DATA REDUCTION PROGRAM NO. 1 


3 LET K =0 

4 FOR Rr 200 TO 400 STEP 200 

5 READ S 

8 FOR I = i TO S 

10 READ PI ,Y1,Y2,T,F,V 

11 LET Hl:%5 

12 LET H2= 4.7 
20 LET Xl = -2 
25 LET X2= 1 
27 LET D= 2.75 

30 LET Hl = H 1 * 1 E-3 
55 LET H2= H2*l E-3 
37 LET Br (H2-H1 )/<X2-Xl ) 

40 LET As H 1 -B*X l 
50 LET K: K+l 
60 LET W0= W* 1 E-6*268 
70 LET H= A+WO 

80 LET P2= <<<Yl+Y2)/2)-52)/66.7 

90 ( LET Ql = (Y2-Y1 )*.938/T 

100 LET P= PI-P2 

110 LET M= EXP (6.68 -.0 104* F) 

115 LET E= B*D/(2*H) 

120 LET 0= (1 . 7 E- 6* O* Q 1 )/(P*Ht3) 

125 LET L= 1.675E-7*R*M/(P*Ht2) 

130 PRINT “RUN" ;K ,"Q1 =“?Q1 ,"HBAR=“;H 
140 PRINT 

150 PRINT "E SUB T="? E, "FLOW=“;Q ,“LAMBDA=“;L 
160 PRINT 

210 NEXT I 

211 PRINT 

212 PRINT 
220 NEXT R 

900 DATA 6 

901 DATA 5,1 .8*2.4,52,81 ,.4 

902 DATA 10,2.9,3.4, 19.2,81 ,.5 

903 DATA 15, 4. 1,5, 19.4,81, .8 

904 DATA 20,6.5,7.5,15.1,81,1 

905 DATA 25,8.5,9.5,11,81,1.3 

906 DATA 30,10.5,12,12.8,81,1.5 


The input data at statement 901 is 

P = 5 psig F = 81°F 

= 1.8 W = .4 

Y 2 ‘ 2 - 4 
T = 52 sec 
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Table IV (continued) 

B. Program Description 

1. Input: PI = supply pressure read from gage, psig 

Y1 = initial reading from graduated supply reservoir 
Y2 = final reading from graduated supply reservoir 
T = time between Y1 and Y2, sec, 

F = oil outlet temperature, °F 
W = Wayne-Kerr proximity probe reading 
HI = initial (static) setting of Sheffield probe No. 1 
H2 - initial (static) setting of Sheffield probe No. 2 
D = Outside diameter of the seal, ip 

2. Output: 

3 

Ql, Q = leakage rate, in /sec 

HBAR, h = average film thickness, in. 

E Sub T, = tilt eccentricity 

FLOW, q = flow coefficient 

3. Calculations 

a. Leakage Rate: Ql = (Y2 - Yl)„938/T, in 3 /sec 

The constant .938 is obtained from the calibration of the graduated 
reservoir. 


b. Average Film Thickness, lu:H = A + WO 


C; 


A = HI - B XI = initial h^ setting 
- 6 

WO = W(26S x 10 ) = Wayne Kerr film thickness correction, in 

B - (H2 - H1)/(X2 - XI) = misalignment slope 
XI = location of HI probe from the center of the seal > in 
X2 = location of H2 probe from the center of the seal, in 


Tilt Eccentricity, e : E SUB T = 

L Zn 


This equation is derived as follows : 

First, the film thickness for the misaligned seal is given by 

h = h . (1 + R e. cos 9). 

1 t 


Therefore (h ) = h. (1 + e ) and (h ) . = h.(l - e . ) . Thus solving 

v o' max 1 N t o min 1 t' 

for the tilt eccentricity w§ have 
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Table IV (continued) 

e. = f(h ) v - (h ) 1 / 2*h . 

t I o max o min l i 

* misalignment slope x outside diameter / 2 
= B D/2A. 

D. Average Supply Pressure: P = PI - P2 

PI = supply pressure read from gage, psig 

P2 = (Yl+Y2)/2 - 52 /66.7 

= hydrostatic pressure head correction 

E. Viscosity: M = EXP (6.68 - .0104F) 

The constants 6.68 and .0104 were obtained from a graph of viscosity vs 
temperature for 350 centistoke silicone oil. 

F. Leakage Coefficient, q: FLOW = 1.7 x 10 ^M(Ql) / P(H)^ 
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TABLE V 


SUMMARY OF ECCENTRIC SEAL TESTS 


Radial Eccentricity 0.125" 

e = .091 
r 

Narrow Land Seal' 




R. = 

.83 




h. 



TABLE 

FIGURE 

in 

ft 

REMARKS 

VI 

21, 23 

.003 

.50 

Outward Pumping(cos ot = -1) 
No Cavity 

VI 

. 22, 23 

.003 

.50 

Inward Pumping(cos ot = 1) 
No Cavity (Oil dam used 
around outside of seal) 

VII 

24, 25 

.004 

.50 

Inward Pumping (cos Oi = 1) 
With Gas Cavities 

II 

26 

.003 

.50 

ii 

II 

27 

.001 

.50 

ii 

II 

28 

.003 

.80 

ii 

II 

29 

.003 

.20 

ii 
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ECCENTRIC, MISALIGNED FACE SEAL DATA 
NARROW LAND SEAL R. = .83 


O r-4 <J“ CO 

o ^ r> on h 


in o vo I s o 

NvTlO O 


os i-h <t m o 
m o os p^. 


cm co co oo oo cm oo co oo cs n m n cn 


oo o cm 


o N CM <}» O CN 00 CM 0 > O H CO 

• •••• •••»• * t • i * 

t-l I s - <t OS »-< CO CM v 0 i—l o CO 4 H CO CO 

co cm H h w Ml- co vo os vo in 


ON <f O 
m vo vo 

H N CO 


cm in in m vo h on in © oo in oo oo <t 


00 O CO O v 0 

^ ^ in m sf 


os n** <f: cm in 

on ts so vp in 


lO ON CM N CM 
CM 0Q I s **. VO v0 


00 H O 

CO ON s}- 


sec 

r*4 

m 

m 

Os 

in 

r-4 

CO 

tH 

T— l 

CO 


vO 

vO 

00 

Os 



co r-v 

CM 

vQ 

vO 

f— 1 

CO 

o 

VO 

CM 


00 

in 

CO 

O 

n 

p— 1 

co <1- 

vO 


CM 

<!■ 

in 


00 

CO 

Mt 

m 


Os 

o 

d 

o 

O o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

O 

O 

o 


• 

• • 

• 

• 

• 

• 

• 

• 

• ' 

• 

• 

• 

• 

• 



I s * CO On O On 

oo I s in m co 

•<t -d* -4- >d- 


H CO CO Mn 
on oo vo m 
<f <3- <t >d- 


sf H CO CO 
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This value was used to calculate 


TABLE' TOif V;. 
ECCENTRIC SEAL BATA?; 

Radial Eccentricity Osi25 M 
e r “ .091 

Narrow '..Land Seal 
R = .83 

1. h i = .004, e t = .50, (Figures 24,25) 


p i 

psi 

N 

r£m 

A 

h 

in 

ft 

Q 

in"*/sec 

_2_ 

q p»mp 

q 

^pump 

^ € r € t C0S 01 

2 

200 

273 

.00402 

.490 

.00357 

12.3 

63.4 

2.56 

4 

it 

157 

.00405 

.487 

.01386 

27.3 

32.8 

2.39 

6 

ti 

109 

.00407 

.484 

.02345 

31.9 

23.2 

2.34 

8 

ti 

84 

.00409 

.482 

.03180 

33.1 

20.7 

2.71 

10 

« 

67 

.00413 

.477 

.04560 

37.5 

11.2 

1.84 

15 

ii 

44 

.00421 

.468 

.07090 

37.6 

10.5 

2.64 

2 

400 

572 

.00399 

.494 

.00397 

14.4 

58.9 

1.13 

4 

ti 

333 

.00397 

.497 

.01400 

29.8 

27.9 

.96 

6 

ii 

225 

.00402 

.490 

.01934 

27.4 

32.5 

1.59 

8 

ii 

171 

.00405 

.487 

.02658 

28.5 

30.4 

1.95 

10 

it 

136 

.00407 

.484 

.03689 

31.2 

24. 7 

1.99 

15 

it 

90 

.00413 

.478 

.06171 

34.2 

18.3 

2.23 

2 

600 

826 

.00397 

.497 

.00695 

24.5 

38.6 

.51 

4 

it 

480 

.00397 

.497 

.00994 

20.3 

46.9 

1.07 

6 

ii 

338 

.00397 

.497 

.01563 

22.5 

42.4 

1.38 

8 

ii 

259 

.00397 

.497 

.02365 

26.1 

35.3 

1.50 

IQ. 

ii 

208 

.00399 

.494 

.03234 

28.5 

30.5 

1.61 

15 

ii 

132 

.00407 

.484 

.05302 

29.1 

29.0 

2.41 

2 

800 

1100 

.00394 

.500 

.00408 

14.5 

58.6 

.59 

4 

ii 

632 

.00394 

.500 

.00690 

14.0 

59.4 

1.03 

6 

ti 

450 

.00391 

.504 

.01078 

15.7 

55.9 

1.37 

8 

ii 

347 

.00389 

.507 

.01876 

21.3 

44.8 

1.36 

10 

ii 

278 

.00386 

.511 

.02255 

20.6 

46.1 

1.82 

15 

ii 

182 

.00394 

.500 

.04188 

24.5 

38.5 

2.32 
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2. = .003, e t = .50, (Figure 26) 


p i 

psi 

5 

N 

rpm 

200 

A 

h. 

l 

in 

ft 

Q 

in~V sec 

_2_ 

4 pump 
e t 

4 pump 

hee,. cos a 
r t 

225 

.00313 

.483 

.00651 

23.8 

40.1 

1.96 

10 

ii 

115 

.00321 

.471 

.01907 

34.5 

17.2 

1.64 

15 

it 

76 

.00327 

.463 

.03272 

38.5 

8.3 

1.20 

20 

n 

55 

.00335 

.452 

.04873 

40.5 

3.1 

.63 

25 

ii 

41 

.00348 

.434 

.06425 

38.3 

6.9 

1.85 

30 

ii 

34 

.00348 

.434 

.08156 

40.7 

1.3 

.42 

■ 5 

400 

467 

.00311 

.487 

.00381 

14.5 

59.0 

1.39 

10 

ii 

236 

.00319 

.474 

.01590 

29.8 

27.3 

1.27 

15 

it 

149 

.00324 

. 467 

.02759 

32.2 

22.2 

1.64 

20 

ii 

108 

.00332 

.455 

.03955 

32.5 

21.1 

2.15 

25 

ii 

84 

.00338 

.448 

.06212 

39.2 

5.9 

.77 

30 

it 

68 

.00343 

.441 

.07346 

37.0 

10.3 

1.67 

5 

600 

725 

.00303 

.500 

.00458 

18.5 

50.4 

.76 

10 

ti 

353 

.00311 

.487 

.01173 

22.6 

42.5 

1.32 

15 

it 

229 

.00319 

.474 

.02132 

26.0 

35.5 

1.70 

20 

ii 

171 

.00321 

.471 

.03518 

31.7 

23.2 

1.49 

25 

ii 

134 

.00327 

. 463 

.04786 

33.1 

20.0 

1.64 

30 

it 

105 

.00332 

.455 

.06700 

35.8 

13.7 

1.43 

5 

800 

999 

.00289 

.523 

.00507 

22.2 

42.8 

.47 

10 

ii 

512 

.00295 

.513 

.01048 

23.1 

41.2 

.89 

15 

. ii 

342 

.00297 

.509 

.01975 

28.9 

30.0 

.96 

20* 

ii 

241 

.00305 

.495 

.02931 

29.4 

28.7 

1.31 

25 

ii 

193 

.00300 

.504 

.03643 

29.8 

28.1 

1.60 

30 

it 

151 

.00311 

.487 

.05003 

30.8 

25.6 

1.86 

1 

200 

763 

.00303 

.500 

.00155 

19.8 

47.8 

.69 

2 

n 

485 

.00303 

.500 

.00295 

24.0 

39.6 

.90 

3 

ii 

355 

.00303 

.500 

. .00426 

25.4 

36. 7 

1.14 

4 

ii 

286 

.00300 

.504 

.00856 

41.5 

4.8 

.19 

5 

it 

236 

.00300 

.504 

.00782 

31.2 

25.2 

1.17 

2 

400 

980 

.00300 

.504 

.00159 

13.2 

61.0 

.69 

4 

it 

562 

.00300 

.504 

.00630 

29.9 

27.7 

• 54 

6 

n 

402 

.00297 

.509 

.00763 

26.2 

35.2 

.97 

.7 

ii 

363 

.00292 

.518 

.00902 

28.5 

30.8 

.94 

4 

600 

821 

.00305 

.495 

.00631 

28.7 

30.0 

.40 

6 

ii 

555 

.00311 

.487 

.00716 

21.6 

44.4 

.88 

8 

it 

423 

.00312 

.485 

.01062 

24.4 

38.8 

1.00 

10 

it 

347 

.00311 

.487 

.01657 

31.3 

24.5 

.78 
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p i 

£Si 

N 

rpm 

A 

h. 

X 

in 

ft 

Q 

in^/sec 

28.9 

q 

^pump 

% 

q 

^pump 

Ae e: cos a 
r t 

4 

800 

1165 

.00297 

.509 

.00581 

29.9 

.28 

6 

ii 

829 

.00295 

.513 

.00763 

27.2 

33.2 

.44 

8 


642 

.00293 

.516 

.00938 

26.1 

35.4 

.61 

10 


522 

.00293 

.516 

.01366 

30.8 

26.2 

.55 

12 

ii 

443 

.00292 

.518 

.01477 

28.4 

30.9 

.77 

14 

ii 

379 

.00293 

.516 

.01753 

28.8 

30.2 

.88 


3. h. 
i 

= .001, 

s = • 

50, (Figure 27 ) 





15 

125 

322 

.00121 

.413 

.00197 

42.4 



20 

11 

210 

.00131 

.383 

.00309 

40.4 



25 

II 

150 

.00139 

.361 

.00545 

48.0 



30 

II 

110 

.00148 

.338 

.00690 

41.8 



40 

II 

66 

.00167 

.300 

.01914 

61.1 



15 

200 

564 

.00116 

.432 

.00211 

52.2 



20 

n 

374 

.00124 

.404 

.00256 

39.2 

2.3 

.07 

25 

ii 

255 

.00135 

.371 

.00364 

35.0 

10.6 

.46 

30 

it 

198 

.00140 

.357 

.00552 

39.6 



40 

it 

112 

.00161 

.310 

.01538 

54.2 



15 

300 

932 

.00110 

.453 

.00105 

30.1 

26.3 

.31 

20 

it 

643 

.00116 

. 432 

.00167 

31.4 

22.8 

.39 

25 

i» 

483 

.00120 

.417 

.00243 

33. 2 

17.9 

.41 

30 

ii 

370 

.00125 

.399 

.00350 

35.1 

12.3 

.37 

40 

n 

207 

.00145 

.344 

.00924 

44.7 



15 

400 

1563 

.00100 

.501 

.0090 

35.6 

16.3 

.11 

20 

ti 

1156 

.00101 

.495 

.00144 

41.7 

3.7 

.035 

25 

ii 

907 

.00102 

.488 

.00167 

37.6 

11.7 

.14 

30 

ii 

686 

.00108 

.464 

.00218 

35.2 

15.4 

.22 

40 

ii 

409 

.00121 

.413 

.00372 

31.9 

20.7 

.55 

10 

o; 

0 

.00122 

.417 

.00153 

48.5 

0 

0 

20 

ti 

11 

.00141 

.365 

.00521 

56.0 

1! 

ft 

30 

ii 

11 

.00159 

.324 

.00932 

47.1 

II 

II 

40 

it 

11 

.00177 

.277 

.01924 

52.6 

It 

11 

50 

ii 

It 

.00197 

.254 

.03442 

55.4 

II 

II 
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4 . h i = .003, e t = .80, (Figure 28 ) 


p i 

psi 

2 

N 

mi 

o' 

A 

h i 

in 

ft 

Q 

.3, 

rn / sec 

a_ 

q 

pump 

^ € r 6 t COS 01 

0 

.00300 

.798 

.01078 

89.1 

0 

4 

It 

ii 

.00305 

.791 

.01363 

62.4 

ii 

6 

Y! 

ii 

.00305 

.787 

.01954 

62.7 

n 

8 



.00309 

.776 

.02565 

61.1 

ti 

10 

Yl 

ii 

.00313 

.771 

.03722 

•69.9 

ti 

15 



.00322 

.751 

.05154 

60.5 

it 

20 

tt 


.00330 

.733 

.07248 

60.0 

n 

2 

200 

527 

.00291 

.818 

.00249 

22.9 

.99 

3 

if 

388 

.00290 

.821 

.00439 

29.8 

1. 10 

4 

1! 

304 

.00291 

.818 

.00586 

31.0 

1.24 

5 

It 

243 

.00293 

.814 

.00821 

34.7 

1.49 

6 

tl 

207 

.00293 

.814 

.01017 

36.6 

1.63 

7 


181 

.00293 

.814 

.01288 

40.3 

1.60 

2 

400 

609 

.00290 

.822 

.00178 

9.5 

1.16 

4 

tt 

610 

.00290 

.822 

.00445 

23.7 

.85 

6 

n 

430 

.00289 

.826 

.00695 

26.2 

1.11 

8 

it 

327 

.00290 

.822 

.01242 

35.6 

1.09 

10 

fi 

261 

.00290 

.822 

.01746 

39.8 

1.15 

12 

it 

216 

.00293 

.814 

.01934 

36.1 

1.60 

14 

it 

184 

.00294 

.811 

.02490 

39.6 

1.60 

5 

200 

212 

.00307 

.775 

.01082 

37.8 

1.42 

10 

it 

111 

.00310 

.769 

.02443 

44.5 

1.80 

15 

it 

72 

.00318 

.749 

.04352 

50.2 

1.49 

20 

ii 

53 

.00323 

.737 

.06212 

51.8 

1.33 

25 

fi 

41 

.00331 

.719 

.08527 

53.2 

. 87 

30 

ii 

33 

.00337 

.708 

.10992 

54.8 

- 

15 

400 

144 

.00318 

.749 

.03703 

42.6 

1.53 

20 

ti 

107 

.00321 

.743 

.05814 

49.6 

1.30 

25 

it 

86 

.00323 

.737 

.07689 

51. 6 

.70 

30 

ii 

68 

.00332 

.719 

.09571 

50.0 

1.24 

35 

ii 

57 

.00337 

.708 

.12129 

52.0 

. 75 

40 

ii 

48 

.00342 

.696 

.15129 

54.3 

- 
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20, (Figure 29 ) 


5. h. = .003, €. - . 
l t 


F i, 

psi 

N 

rpm 

A 

1 

in 

e t 

• 3 / 
in /sec 

_a_; 

Ae e^cos a 
r t 

5 

0 

0 

.00307 

.195 

.01150 

41.1 

0 

10 

II 


.00317 

.188 

.02233 

39.0 


15 

11 

It 

.00325 

.180 

.03518 

38.7 


20 

11 

11 

.00333 

.176 

.05133 

39.9 


25 

11 

It 

.00342 

.173 

.06764 

39.2 


30 

11 


.00351 

.170 

.08716 

39.3 


35 

11 

If 

.00357 

.168 

.10619 

39.2 


40 

11 


.00367 

.163 

.13132 

39.1 


5 

100 

104 

.00313 

.191 

.00788 

26.7 

4.9 

10 

11 

53 

.00321 

.186 

.02132 

35.8 


15 

II 

34 

.00329 

.181 

.03350 

35.7 

- 

20 

11 

25 

.00334 

.178 

.04841 

37.3 

- 

25' 

11 

19 

.00344 

.173 

.06700 

38.3 

- 

30 

11 

15 

.00353 

.169 

.08527 

37.6 

- 

35 

11 

13 

.00358 

.166 

.10196 

37.0 


40 

II 

11 

.00368 

. 162 

.12849 

37.9 

- 

5 

100 

104 

.00313 

.191 

.00885 

29.9 

3.03 

10 

ft 

53 

.00321 

.186 

.01927 

32.3 

3.36 

15 

It 

34 

.00329 

.181 

.03212 

34.1 

1.95 

20 

II 

25 

.00334 

.178 

.04425 

34.0 

2.94 

25 

It 

19 

.00344 

.173 

.06194 

35.4 

-■ 

30 

11 

15 

.00353 

.167 

.07949 

35.1 

- 

35 

11 

13 

.00358 

.166 

.09670 

35.1 

- 

40 

11 

11 

.00368 

.162 

.12342 

36.4 

•• - ■ 

5 

200 

212 

.00310 

.192 

.00929 

32.2 

.87 

10 

1! 

110 

.00315 

.189 

.01857 

32.8 

1.40 

20 

11 

51 

.00332 

.180 

.04467 

35.2 

. 24 

30 

11 

32 

.00348 

.171 

.07688 

35.5 

- 

40 

11 

22 

.00361 

.165 

.12563 

39.1 

- 

50 

11 

16 

.00377 

.158 

.17370 

38.1 

- 

60 

It 

13 

.00393 

.152 

.23450 

37.9 

. - 

5 

200 

212 

.00310 

.192 

.00634 

22.0 

3.66 

10 

ii 

110 

.00315 

.189 

.01770 

31. 2 

2.18 

20 

it 

51 

.00332 

.180 

.04425 

34.9 

.42 

30 

it 

32 

.00348 

• .171 

.07817 

36.1 

- 

40 

ti 

22 

.00361 

.165 

,12026 

37.-4 

- 

50 

ii 

16 

.00377 

.158 

.16456 

36.1 


60 

ii 

13 

.00393 

.152 

,23848 

38.6 

- 
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TABLE VIII 


DATA REDUCTION PROGRAM No. 2 


A. Listing 

10 PRINT "RUN", "Q PUMP ", "LAMBDA" 

15 PRINT 

20 FOR I - 17 TO 23 
30 READ E,Q,L 
40 LET Q 1 - 25 + 37.5*E 
50 LET P = (Q1-Q)/(E*.091*L) 

60 PRINT I,P,L 
70 NEXT I 
900 DATA .191 

B. Program Description 

This program takes the output from data reduction program No. 1, q,A > £ t » 

£ r , cos o(, and calculates the pumping variable, qpump/ A^r^t cos o( * The 
pumping coefficient, qp U mp> f° un d by subtracting the theoretical q for 
£ r = 0 from the measure (input) q. The theoretical value of q for G r = 0 is 
calculated in the program from an approximate equation for q(£ t ) for the misaligned 
seal. 
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TABLE IX 


ECCENTRIC MISALIGNED FACE SEAL DATA 
NARROW LAND SEAL, Ri = 0.83 
RADIAL ECCENTRICITY 0.125", e r = 0.091 
COS a = 1, INWARD PUMPING 


1. No cavities (oil/oil), 350/350 centistokes 


Oil 

Level 

Y1-Y2 

N 

rpm 

A 

hi 

in 

it 

,Q 

in -/sec 

£ 

58.5 

- 53.5 

0 

0 

0.00290 

0.5690 

0.00109 

45.99 

52 

- 47 

0 

0 

0.00290 

0.5690 

0.00103 

44.97 

38 

- 38 

22 

298 

0.00290 

0.5690 

0.0 

0.0 

39 

- 39 

40 

537 

0.00290 

0.5690 

-0.00096 

- 45.33 

42 

- 45 

66 

863 

0.00290 

0.5690 

-0.00196 

- 89.86 

41 

- 40 

15 

200 

0.00290 

0.5690 

0.00029 

13.48 

42 

- 47 

47 

615 

0.00290 

0.5690 

-0.00161 

- 53.15 

44 

- 43.7 

21 

273 

0.00290 

0.5690 

0.00010 

4.72 

44 

- 45 

35 

453 

0.00290 

0.5690 

-0.00064 

- 28.91 

46 

- 49 

58 

737 

0.00290 

0.5690 

-0.00167 

- 74.17 

52 

- 55 

82 

1003 

0.00290 

0.5690 

-0.00252 

-107.71 

57 

- 61 

108 

1310 

0.00286 

0.5723 

-0.00330 

-141.91 


2. No cavities 

(oil/oil), 350/200 centistokes 



68 - 67 

28 

316 

0.00292 

0.5629 

0.00063 

24.59 

58.5 - 57 

12 

144 

0.00291 

0.5648 

0.00132 

55.29 

56.5 - 56 

29 

352 

0.00291 

0.5648 

0.00037 

15.84 

56 - 56.5 

49 

595 

0.00291 

0.5648 

-0.00031 

- 12.92 

57.5 - 60.5 

115 

1373 

0.00291 

0.5648 

-0.00296 

-123.04 

62 - 63 

60 

701 

0.00291 

0.5648 

-0.00081 

- 32.96 

64.5 - 66.5 

87 

999 

0.00291 

0.5648 

-0.00090 

- 36.09 

68 - 70.5 

127 

1428 

0.00291 

0.5648 

-0.00272 

-106.58 

73 - 77 

185 

2052 

0.00289 

0.5732 

-0.00532 

-207.30 

58 - 57 

18 

203 

0.00292 

0.5627 

0.00083 

32.55 

56.5 - 55.5 

19 

216 

0.00292 

0.5627 

0.00073 

29.01 

55 - 54 

20 

229 

0.00292 

0.5627 

0.00076 

30.22 

51 - 50 

15 

176 

0.00292 

0.5627 

0.00079 

32.19 

49 - 48 

12 

143 

0.00292 

0.5627 

0.00081 

33.62 

48 - 48 

33 

393 

0.00292 

0.5627 

0.0 

0.0 

48.5 - 49.5 

57 

675 

0.00292 

0.5627 

-0.00091 

- 37.51 

50.5 - 52.5 

80 

940 

0.00291 

0.5646 

-0.00188 

- 77.21 
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3, No cavities 

(oil/oil) , 350/100 centistokes 



Oil 







Level 

N 


hi 


Q 


Y1-Y2 

r£tn 

A 

in 

it 

in^/ sec 

£ 

58 - 57 

18 

214 

0.00290 

0.5667 

0.00042 

17.46 

57.5 - 58.5 

80 

950 

0.00290 

0.5667 

-0.00072 

- 29.93 

59 - 59 

47 

555 

0.00290 

0.5667 

0.0 

0.0 

61 - 64 

116 

1340 

0.00290 

0.5667 

-0.00253 

-102.35 

68.5.- 69.5 

78 

882 

0.00288 

0.5705 

-0.00043 

- 17.29 

66.5 - 65 

20 

230 

0.00288 

0.5705 

0.00095 

38.64 

64.5 - 64 

47 

546 

0.00288 

0.5705 

0.00028 

11.39 

64 - 65 

88 

1021 

0.00288 

0.5705 

-0.00091 

- 37.28 

63 - 63.5 

77 

894 

0.00289 

0.5686 

-0.00032 

- 13.03 

66 - 69 

126 

1427 

0.00289 

0.5686 

-0.00234 

- 93.02 

73.5 - 78.5 

166 

1810 

0.00287 

0.5704 

-0.00407 

-156.68 

*61 - 61 

52 

608 

0.00290 

0.5619 

0 

0 

62 - 64 

109 

1259 

0.00290 

0.5619 

-0.0018 

- 72.89 

**65 - 65.5 

71 

822 

0.00288 

0.5680 

-0.00026 

- 10.68 

59.5 - 61 

91 

1091 

0.00287 

0.5677 

-0.00095 

- 40.13 

53 - 54 

71 

887 

0.00287 

0.5677 

-0.00077 

- 33.89 

55 - 58 

139 

1704 

0.00287 

0.5677 

-0.00288 

-124.85 

**39 - 38 

12 

164 

0.00285 

0.5741 

0.00073 

35.62 

37 - 36 

11 

153 

0.00285 

0.5741 

0.00060 

29.68 

35.5 - 34.5 

11 

156 

0.00284 

0.5760 

0.00068 

34.43 

**76.5 - 76 

31 

341 

0.00284 

0.5760 

0.00020 

8.01 

75.5 - 75 

32 

354 

0.00284 

0.5760 

0.00023 

9.01 

75.5 - 74 

33 

367 

0.00284 

0.5760 

0.00044 

17.31 

74 - 73 

32 

357 

0.00284 

0.5760 

0.00046 

18.43 

**65.5 - 64.5 

25 

293 

0.00284 

0.5760 

0.00057 

24.04 

64 - 63 

26 

307 

0.00284 

0.5760 

0.00050 

20.93 

63 - 62 

26 

311 

0.00283 

0.5781 

0.00036 

15.59 


4. No 

cavities 

(oil/oil) , 350/20 centistokes 



**42.5 - 

41.5 

21 

279 

0.00294 

0.5706 

0.00055 

25.09 

41 - 

40.5 

56 

751 

0.00294 

0.5706 

0.00027 

12.43 

40 - 

40 

95 

1301 

0.00092 

0.5720 

0 

0 

**55 - 

55 

125 

1568 

0.00290 

0.5762 

0 

0 

55 - 

56 

168 

2101 

0.00290 

0.5762 

-0.00053 

- 23.00 

56.5 - 

57.5 

205 

2540 

0.00290 

0.5762 

-0.00085 

- 36.83 

59 - 

61 

310 

3834 

0.00288 

0.5825 

-0.00266 

-115.84 
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5. No cavities (oil/oil), 350/10 centistokes 


Oil 







Level 

N 


hi 


Q 


Y1 

-Y2 

rpm 

A 

in 

ft 

in-Vsec 

£ 

70.5 

- 69.5 

23 

255 

0.00297 

0.5600 

0.00098 

37.31 

65 o 5 

- 64.5 

95 

1092 

0.00296 

0.5643 

0.00075 

29.58 

64 

- 63 

168 

1949 

0.00296 

0.5643 

0.00055 

22.00 

63 

- 63.2 

246 

2860 

0.00296 

0.5643 

-0.00015 

- 5.88 

61.6 

- 61.8 

249 

2919 

0.00296 

0.5643 

-0.00015 

- 6.22 

62 

- 62.2 

303 

3544 

0.00296 

0.5643 

-0.00015 

- 6.14 

62.5 

- 63 

355 

4319 

0.00289 

0.5758 

-0.00052 

-22.08 

59.5 

- 58.5 

101 

1257 

0.00289 

0.5758 

-0.00070 

30.61 

58 

- 57 

149 

1871 

0.00289 

0.5758 

-0.00081 

35.87 

57 

- 57 

291 

3664 

0.00289 

0.5758 

0 

0 

57.5 

- 58.5 

434 

5508 

0.00287 

0.5774 

-0.00077 

-34.41 


6 . No cavities (oil/oil), 350/5 centistokes 


**47 

- 

46 

37 

492 

0.00290 

0.5712 

0.00068 

31.71 

45.5 

- 

45 

70 

938 

0.00290 

0.5712 

0.00057 

21.49 

44 

- 

43.5 

127 

1719 

0.00290 

0.5712 

0.00043 

20.23 

**39 

- 

38 

60 

836 

0.00290 

0.5712 

0.00081 

39.36 

37.5 

- 

37 

144 

2023 

0.00291 

0.5692 

0.00046 

22.44 

37 

- 

36.5 

238 

3355 

0.00291 

0.5692 

0.00047 

23.19 

36.5. 


37 

278 

3919 

0.00291 

0.5692 

-0.00032 

-15.81 

33 

- 

33 

284 

4111 

0.00291 

0.5692 

0.0 

0 

33 

- 

33.1 

337 

5002 

0.00287 

0.5749 

-0.00007 

- 3.44 


7. No cavities (oil/oil) , 350/2 centistokes 


46 - 

44 

40 

490 

0.00291 

0.5645 

0.00118 

50.54 

43 - 

42 

102 

1268 

0.00291 

0.5645 

0.00104 

45.27 

41 - 

40- 

144 

1812 

0.00291 

0.5645 

0.00097 

43.33 

39 - 

38 

163 

2077 

0.00291 

0.5645 

0.00102 

45.29 

37 - 

36 

228 

2942 

0.00291 

0.5645 

0.00094 

44.45 

35.5 - 

35 

316 

4110 

0.00291 

0.5645 

0.00053 

24.09 

34.5 - 

34 

344 

4577 

0.00288 

0.5707 

0.00049 

23.10 
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8. No cavities (oil/oil), 350/0.65 centistokes 


Oil 

Level 

Y1-Y2 

N 

IE! 2 

A 

hi 

in 

it 

,Q 

in J /sec 

R 

56.5 - 55.5 

105 

1321 

0.00290 

0.5738 

0.00086 

37.78 

53 - 51.5 

141 

1816 

0.00290 

0.5738 

0.00174 

78.21 

50.5.- 48.5 

182 

2385 

0.00290 

0.5738 

0.00189 

86.64 

47 - 45.5 

225 

3011 

0.00290 

0.5738 

0.00172 

80.53 

**42.5.- 40.5 

142 

1929 

0.00293 

0.5700 

0.00180 

84.68 

39 - 37 

302 

4202 

0.00293 

0.5700 

0.00164 

79.38 

36 - 34.5 

395 

5602 

0.00293 

0.5700 

0.00160 

78.93 

33.5 - 32 

430 

6082 

0.00296 

0.5643 

0.00135 

65.68 


9. Cavities 

(oil/air) 3 

350 centistokes/air 




99.5 - 98.5 

21 

204 

0.00292 

0.5676 

0.00064 

21.72 

98 - 97 

51 

501 

0.00292 

0.5676 

0.00040 

13.65 

95 - 94 

65 

648 

0.00292 

0.5676 

0.00029 

9.95 

94 - 93 

126 

1272 

0.00291 

0.5694 

0.00027 

9.52 

91 - 90 

150 

1537 

0.00291 

0.5694 

0.00026 

9.42 

90 - 89 

173 

1786 

0.00291 

0.5694 

0.00032 

11.52 

87 - 86 

201 

2107 

0.00291 

0.5694 

0.00023 

8 » 49 

86 - 85 

240 

2594 

0.00287 

0.5773 

0.00027 

10.43 

81 - 80 

267 

3004 

0.00285 

0.5789 

0.00022 

8.81 

77 - 76 

307 

3554 

0.00284 

0.5834 

0.00024 

9.73 

74 - 73 

343 

4104 

0.00282 

0.5899 

0.00031 

13.24 

77 - 72 

0 

0 

0.00291 

0.5692 

0.00151 

58.52 

71 - 66 

0 

0 

0.00292 

0.5696 

0.00149 

59.04 

65 - 60 

0 

0 

0.00293 

0.5700 

0.00138 

56.06 

58 - 53 

0 

0 

0.00293 

0.5703 

0.00138 

57.88 

95 - 90 

0 

0 

0.00292 

0.5698 

0.00161 

55.93 

88 - 83 

0 

0 

0.00292 

0.5698 

0.00153 

54.89 

82 - 77 

0 

0 

0.0Q292 

0.5698 

0.00152 

56.11 

96 - 94 

0 

0 

0.00293 

0.5700 

0.00161 

54.74 

93 - 83 

0 

0 

0.00293 

0.5700 

0.00155 

54. 75 

81 - 79 

0 

0 

0.00293 

0.5700 

0.00153 

56.25 

75 - 65 

0 

0 

0.00293 

0.5700 

0.00141 

54.71 

63 - 61 

0 

0 

0.00293 

0.5702 

0.00134 

54.43 
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10. No cavities (oil/oil), 100/100 centistokes 


Oil 







Level 

N 


hi 


Q 


Yl 

-Y2 

rpm 

A 

in 

f_t 

in^/sec 

£ 

73 

- 63 

0 

0 

0.00291 

0.5694 

0.00362 

42.63 

60 

- 50 

0 

0 

0.00291 

0.5694 

0.00344 

43.72 

48 

- 38 

0 

0 

0.00291 

0.5694 

0.00321 

44.09 

36 

- 26 

0 

0 

0.00291 

0.5694 

0.00301 

44.99 

72.5 

- 67.5 

10 

34 

0.00289 

0.5733 

0.00321 

38.21 

65 

- 62 

21 

74 

0.00289 

0.5733 

0.00255 

31.49 

53.5 

- 51.5 

43 

162 

0.00289 

0.5733 

0.00153 

20.19 

51 

- 51 

80 

303 

0.00289 

0.5733 

0 . 

0 . 

51.5 

- 53.5 

116 

436 

0.00289 

0.5733 

-0.00194 

-25.58 

54.5 

- 58 

160 

587 

0.00289 

0.5733 

-0.00368 

-47.48 


11 . 

No cavities 

(oil/oil) , 100/350 centistokes 



52 


50 

15 

56 

0.00290 

0.5691 

0.00184 

23.87 

49 


48 

35 

132 

0.00290 

0.5691 

0.00184 

8.48 

48 


47 

45 

170 

0.00290 

0.5691 

0.00049 

6.46 

47.5 

_ 

51 

122 

457 

0.00290 

0.5691 

0.00343 

-44.86 

53 

- 

54.5 

65 

237 

0.00290 

0.5691 

- 0.00122 

-15.55 

56 

- - 

59 

131 

471 

0.00289 

0.5757 

-0.00335 

-42.88 

**55 

- 

54 

22 

80 

0.00289 

0.5757 

0.00123 

15.75 

54.5 


57 

88 

319 

0.00289 

0.5757 

-0.00250 

-31.81 

58.5 

_ 

61 

100 

354 

0.00289 

0.5757 

-0.00298 

-37.06 

62.5 

- 

65 

162 

566 

0.00288 

0.5752 

-0.00339 

-41.83 


12. No cavities 

(oil/oil), 100/200 centistokes 



87.5 - 83.5 

30 

92 

0.00291 

0.5692 

0.00306 

32.91 

80 - 76 

32 

103 

0.00291 

0.5692 

0.00281 

31.39 

74 - 74 

74 

242 

0.00291 

0.5692 

0 . 

0 . 

74 - 71 

35 

116 

0.00291 

0.5692 

0.00244 

28.08 

68.5 - 65 

19 

65 

0.00291 

0.5692 

0.00312 

37.19 

**55 - 55 

64 

237 

0.00289 

0.5733 

0 . 

0 . 

55.5 - 56.5 

70 

258 

0.00289 

0.5733 

-0.00065 

- 8.45 

57.5 - 59 

90 

326 

0.00289 

0.5733 

-0.00148 

- 18.85 

**56.5 - 61.5 

99 

355 

0.00288 

0.5753 

-0.00431 

- 54.40 

64-68 

149 

520 

0.00288 

0.5753 

-0.00457 

- 56.20 

71 - 74 

196 

659 

0.00288 

0.5753 

-0.00421 

- 49.88 
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13. Cavities (oil/air), 100 centistokes/air 


Oil 

Level 

Y1-Y2 

N 

rpm 

A 

hi 

in 

it 

Q 

in^/sec 

3 . 

62.5 - 

59.5 

35 

123 

0.00289 

0.5665 

0.00223 

27.33 

57 - 

54 

14 

51 

0.00289 

0.5665 

0.00288 

36.53 

52 - 

50 

51 

190 

0.00289 

0.5665 

0.00164 

21.45 

49 - 

47 

74 

281 

0.00289 

0.5665 

0.00127 

16.87 

46 - 

44.5 

77 

297 

0.00289 

0.5665 

0.00121 

16.32 

43.5 - 

42.5 

133 

522 

0.00289 

0.5665 

0.00094 

12.84 

37 - 

36 

129 

533 

0.00288 

0.5683 

0.00081 

11.86 

35.5 - 

34.5 

164 

687 

0.00288 

0.5705 

0.00069 

10.24 

67 - 

66 

284 

995 

0.00285 

0.5789 

0.00060 

7.52 

51.5 - 

50.5 

318 

1223 

0.00285 

0.5789 

0.00066 

9.03 

50 - 

49 

372 

1461 

0.00283 

0.5832 

0. 00051 

7.13 
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TABLE X 


HIGH SPEED INERTIA TESTS 
CONCENTRIC MISALIGNED FACE SEAL DATA 
NARROW LAND SEAL, Ri = 0.83 


1. 10 centistoke oil, e t = 0.0 


Oil 








Level 

N 


hi 

o 2 

,Q 



Y1-Y2 

rpm 

A 

in 

e t 

in' 3 /sec 

«L 

Pcorr 

95 - 75 

0 

0 

0.00300 

9.25E-6 

0.0346 

1.010 

0.0066 

65 - 55 

0 

0 

0.00300 

9.25E-6 

0.0346 

1.010 

0.0066 

45 - 25 

0 

0 

0.00300 

9.25E-6 

0.0258 

1.016 

0.0080 

80 - 60 

250 

74 

0.00300 

2.34E-6 

0.0388 

1.165 

0.1169 

90 - 70 

500 

142 

0.00298 

0 . 

0.0435 

1.171 

0.1702 

100 - 80 

750 

203 

0.00297 

0 . 

0.0486 

1.134 

0.2167 

100 - 80 

1000 

274 

0.00295 

0 . 

0.0537 

1.124 

0.3037 

100 - 80 

1250 

348 

0.00293 

0 . 

0.0601 

1.110 

0.4104 

100 - 80 

1500 

417 

0.00292 

2.47E-6 

0.0664 

1.057 

0.4895 

100 - 80 

1650 

464 

0.00290 

0 . 

0.0721 

1.064 

0.5934 

Realign Spindle 







100 - 80 

1000 

266 

0.00299 

9.38E-6 

0.0561 

1.124 

0.3034 

100 - 80 

2000 

535 

0.00297 

5.94E-7 

0.0935 

1.031 

0.7938 

75 - 55 

2500 

778 

0.00295 

5.43E-6 

0.1098 

0.9594 

1.0617 

90 - 50 

3000 

924 

0.00291 

5.58E-6 

0.1403 

0.9511 

1.5139 

90 - 50 

4000 

1256 

0.00287 

5.74E-6 

0.2020 

0.8745 

2.3388 

80 - 40 

5000 

1816 

0.00275 

2.78E-6 

0.2295 

0.7640 

2.8972 


2. 1C 

i centistoke 

oil, et 

= 0.56 





100 - 

80 

0 

0 

0.00288 

0.3065 

0.0468 

1.4627 

0.0514 

95 - 

75 

250 

75 

0.00287 

0.3087 

0.0468 

1.4797 

0.0715 

100 - 

80 

500 

146 

0.00286 

0.3112 

0.0515 

1.4905 

0.1147 

100 - 

80 

750 

215 

0.00286 

0.3112 

0.0561 

1.4361 

0.1464 

100 - 

70 

1000 

291 

0.00286 

0.3112 

0.0631 

1.4321 

0.2295 

100 - 

70 

1500 

440 

0.00285 

0.3132 

0.0815 

1.3830 

0.4281 

80 - 

40 

2000 

675 

0.00284 

0.3180 

0.0990 

1.3019 

0.6503 

40 - 

10 

2500 

1082 

0.00280 

0.3275 

0.1222 

1.2795 

0.9893 

80 - 

40 

3000 

1021 

0.00280 

0.3275 

0.1578 

1.1727 

1.2002 

* Reservoir v 

first 

flooded with 350 < 

centistoke 

oil, siphoned 

out and 

then fill 


with 100 centistoke oil. 


**Siphon out and refill reservoir. 
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TABLE XI 


SEAL CONDITIONS 
FOR THE 

CAVITY SHAPE EXPERIMENTS 
(SEE FIGURES 8(a) - (s>) 
R t = 0.83, h = 0.003 in 


1. e = 0.56, e = 0.091 
t r 


Oil 


Level 

Y1-Y2 

N 

rpm 

A 

hi 

in 

it 

,Q 

irrVsec 

£ 

Figure 

No*.., 

79 - 

78 

39 

432 

0.00293 

0.5633 

0.00052 

20.02 

8(a) 

77.5 - 

76.5 

101 

1128 

0.00293 

0.5633 

0.00053 

20.45 

8(b) 

76 - 

75 

160 

1801 

0.00293 

0.5633 

0.00051 

19.71 

8(c) 

74.5 - 

73.5 

190 

2197 

0.00291 

0.5692 

0.00048 

19.25 

8(d) 

69 - 

68 

234 

2793 

0.00291 

0.5692 

0.00046 

19.10 

8(e) 

67.5 - 

66.5 

298 

3647 

0.00288 

0.5730 

0.00044 

18.75 

8(f) 

85.5 - 

84.5 

173 

1866 

0.00292 

0.5651 

0.00026 

9.79 

8(g) 

80 - 

79 

34 

378 

0.00292 

0.5651 

0.00047 

18.09 

8(h) 

78.5 - 

77.5 

78 

874 

0.00292 

0.5651 

0.00031 

11.99 

8(i) 

77 - 

76, 

126 

1423 

0.00292 

0.5651 

0.00029 

11.25 

8(j) 

75.5 - 

74.5 

210 

2409 

0.00291 

0.5623 

0.00025 

9.93 

8(k) 

70.5 - 

69.5 

281 

3325 

0.00291 

0.5621 

0.00028 

11.76 

8(1) 



600 


-- 




8(m) 



600 


— 

-- 

-- 

-- 

8(n) 

58 - 

57 

385 

5134 

0.00284 

0.5739 

0.00026 

12.18 

8(o) 

53 - 

52 

503 

7057 

0.00282 

0.5801 

0.00036 

18.44 

8(P) 


2. e fc = 0.56, 

Oil 

Level 

Y1-Y2 

e = 
r 

N 

rpm 

0.0 

A 

hi 

in 

it 

Q 

in 3 /sec 

£ 

Figure 

No* 

62 - 60 

31 

376 

0.00294 

0.5658 

0.00143 

60.09 

No 

59 - 57 

121 

1497 

0.00293 

0.5656 

0.00138 

59.24 

Cavity 

No 

53-51 

154 

1980 

0.00293 

0.5656 

0.00134 

59.75 

Cavity 

No 

49 - 47 

203 

2680 

0.00293 

0.5656 

0.00136 

62.34 

Cavity 

No 

46 - 44 

241 

3246 

0.00293 

0.5656 

0.00144 

67.25 

Cavity 

No 

43 - 41 

312 

4360 

0.00291 

0.5670 

0.00134 

65.45 

Cavity 

8(q) 

37-35 

380 

5541 

0.00291 

0.5670 

0.00122 

62.03 

8(r) 

34 - 32 

505 

7706 

0.00288 

0.5728 

0.00120 

60.30 

8(s) 
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TABLE XII 

DATA REDUCTION PROGRAM NO. 3 


A, Listing 

1 PRINT " ORDER OF OUTPUT” 

2 PRINT 

3 LET K* 0 

4 PRINT ” RUN NO”# "FLOW”# ”H BAR”# ”E SUB T»2' 

5 READ S# SO# ED 

6 PRINT 

7 PRINT ” 9”*”LAMBDA 1”#”0PUMP”#"9PUMP/L*£R*E" 

8 FOR I* I TO S 

10 READ R#Y1#Y2#T#H1#H2#F 
13 LET P1»0. 

15 LET W* 0. 

20 LET Xla-2 

25 LET X2* 1 

27 LET D* 2.75 

30 LET Hla HI* IE- 3 

35 LET H2* H2*lE-3 

37 LET B* CH2-H1 )/(X2-Xl) 

40 LET A* HI -B*X1 
50 LET K* K+1 
60 LET WOa W* IE- 6*268 
70 LET H* A* WO 

80 LET P2a <<<Yl*Y2)/2>+107.5>*.00351 

85 LET PO* 2.68E-7*R»2 

90 LET 91« (Y1-Y2)*. 02525/ T 

100 LET Pa P1+P2+P3 

110 LET Ma EXP<2.9876-.0089*F) 

US LET E* B*D/C2*H) 

120 LET 9* 1.76E-6*M*91/CP*H» 3)*. 0296 
125 LET LI « 1 . 67E-7*R*M/<CP-P3)*H» 2) 

130 LET 90* l • +SO*Et 2 
135 LET 99» 9/90 
140 LET R1«0. 

150 PRINT 

160 PRINT K#91#H#E*2 
165 PRINT 9#L1#99#R1 

210 NEXT I 

211 PRINT 

212 PRINT 

900 DATA 9# 1. 18#0 

B. Program Description 

This program is very similar to data reduction program No. 1 but it has the 
inertia component of leakage added to it. This can be seen as P3 in the pro- 
gram listing. This program was used to reduce the data to that shown in Table X. 
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NOMENCLATURE 


A - seal area, in^ 

2 

A^ - flow area, in 

h - film thickness, in 

- average film thickness at inside radius, in 

e - displacement between seal center and shaft center, in 

y - coordinate across the film, in 

p - coordinate measured from center; of. shaft, in 

3 

p - density, lb/in 

r - coordinate measured from center of seal, in 

0 - angular coordinate, radian 

P - a + 9 , radian 

a - phase angle between waviness and eccentricity, radian 

y - angle measured to p , radian 

i|r - tilt phase angle, radian 

t - time, seconds 

U - fluid velocity, in/sec 

tw - rotational speed, radian/sec 

3 

Q - flow rate, in /sec 

p - pressure, psi 

T - shear stress, psi 

W - seal load capacity, lb^ 

SM - shaft friction moment, in- lb^ 

0^ - specific heat, Btu/lb°F 

k - thermal conductivity, Btu/sec in °F 

T - temperature, °F 

$ - dissipation function, sec ^ 
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M> - dynamic viscosity, psi-sec 
Dimensionless Terms 


R 

H 

P 

A 

e 

e„ 


U' 

I 

X 


radius ratio: = r/r 


film thickness ratio = h/E, 


pressure ratio - p/p 


seal number 


( !o ) 2 

P o V 


waviness parameter 
radial eccentricity 
tilt eccentricity 
load coefficient = - 


W 


P r 
r o o 


friction coefficient 


SM 


r3 2 #, 
h.p r /6 
i*o o 


flow coefficient 


12 UQ 

r3 

h i p o 


pumping velocity variable 


U 

r-<« 

o 


radius variable * 


r -r. 
o i 


b “ h 

film thickness variable *> o i 


n - number of circumferential waves 

E^, - numerical grouping 

- numerical grouping 
Subscripts 

0 - at the outside radius of the seal 

1 - at the inside radius of the seal 

C - within the cavitated film 
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APPENDIX D 
SEAL DESIGN I 
BY 

J. A. FINDLAY 


I. INTRODUCTION 

The work on this task has been concerned, primarily with the extension of 
the basic theories of seal cavitation to the misaligned eccentric seal and thus 
to include inward pumping. Expressions have been derived so that either the 
"wavy" seal or the tilted seal may be considered. The theories derived in this 
task have been used in making supporting calculations for the cavitation and 
pumping experiments of Task III. This was accomplished using a GE 635 computer 
program. The iterative method used in this program was improved so that a 
numerical solution to the short bearing approximation is obtained in one sweep. 
Subsequent iterations include the two-dimensional effects. Convergence is 
attained in fewer iterations than the previous method. Work has also been done 
towards modifying this computer program so that cavity solutions can be obtained 
for the eccentric seal. 

Further analytical work has been cbne on the two-fluid problem. This con- 
siders the case when the fluids on either side of the face seal are of different 
viscosities. The effect of the viscosity of both the inside and outside fluids 
on the leakage rates has been studied. 


II. INWARD PUMPING IN MECHANICAL FACE SEALS 

A. INTRODUCTION 

Hydrodynamic effects have been shown to exist in mechanical face seals. 

Seal forces and leak rates are primary considerations in the design and applica- 
tion of face seals. Therefore, it is important to understand the basic princi- 
ples of face seal operation. (See Figure 1.) 

An anomolous behavior observed in face seal experiments by Denny [l] , 
Nauf^and others, is that of "inward pumping". This is a condition where fluid 
or air leaks through the seal in the opposite direction to the hydrostatic 
pressure drop across the seal. The term "inward pumping" was attached to this 

^-Numbers in parentheses refer to similarly numbered references in bibliography 
at end of Appendix. 
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phenomena because in the seal experiments where it was observed the inside of 
the seal was pressurized. Therefore, any leakage in the outward direction was 
attributed to the hydrostatic pressure differential while observed leakage 
toward the center of the seal and against the pressure was attributed to a 
pumping action in the seal. 

So far no satisfactory explanation of this phenomenon has been set forth. 
Nau|sj in a recent paper based on experiments, using a transparent seal apparatus 
confirmed his earlier postulate [ 4 ] that lubricant cavitation was the cause for 
"inward pumping". He also concluded that this pumping action occurred only in 
the inward direction. 

Lyman and Saibel £sj, using an approximate solution to the face seal hydro- 
dynamic problem with gas cavities, were able to calculate an "inward pumping" 
rate for the case of no pressure difference across the seal. 

Findlay jbj, £ 7 ] and showed analytically that cavity solutions using 
the Swif t-Floberg ^14j , ^15| boundary conditions, i.e., conservation of flow at 

the cavity boundaries, did not result in any "inward pumping". Subsequent 
experimental results reported here and in [ 9 ] confirm this conclusion and 
substantiate the short bearing theory put forth in [^j . 

A theory is developed here which shows that a combination of radial eccen- 
tricity and misalignment, or one circumferential film thickness wave will result 
in a "pumping" leakage which may be either inward or outward. Experimental 
results reported here and in ^ 9 ] confirm this theory, which was developed for a 
full film, i.e., no gas cavities. Experiments and theory for the eccentric seal 
with gas cavities show that rather than causing "inward pumping" the cavities 
reduce the pumping effect. 

B. DERIVATION OF THE SHORT BEARING PUMPING SOLUTION 
1. General Equations for the Wavy Seal 

Reynolds* equation for the radially eccentric face seal is 


h Hi + I !e n 3 H = A <fl +te [ HU e 


where, from Fig. 2, the relative seal surface velocities are given by the follow- 
ing expressions 


= e r cos £ 


U(J = R + e r sin 


Using the short bearing approximation, where the dp/90 term is neglected, 
and assuming that the film thickness is a function of 0 only, as in the wavy seal 
case, equation (1) reduces to 



Integration of this expression gives the following general solution, 


A / ' dH ^1 

r '^\- +Re r 3l “ S ] d0 + ^3 /nE+ C 2 


where and are constants of integration which are determined from the 
boundary conditions on the pressure. 

For the case where gas cavities may be present in the fluid film j^6 J , there 
are three separate film regions and three sets of boundary conditions where 
equation (5) must be applied. These regions, shown in Fig. 3, are: 

1. The full film region where the film extends across the face of the 
seal. 

2. The outer film region where the film extends from the. outer edge 
of the cavity to the outer edge of the seal. 

3. The inner film region where the film extends from the inner edge 
of the seal to the inner edge of the cavity. 

There is also the cavity region where the pressure is constant and 
equations (4) and (5) do not apply. 

The general form of the boundary conditions is 

P = P @ R = R and P = P @R = R (6) 

m m n n 

where the specific values for each region are given by the following table: 



p 

p 

R 

R 

R 


m 

n 

m 

n 

c 

Full film region 

P. 

i 

P 

o 

R. 

X 

1 

- 

Inner film region 

P. 

X 

P 

c 

R. 

i 

R 1 

R i 

Outer film region 

P 

c 

p 

o 

R 2 

1 

R 2 
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The pressure and radial pressure gradient in any region have the following 
general form 

^n(R/R ) A HH 1 

(P “ V “/n(R /rT') (P ni Pn) + \ (R2 “ 

x m n 4H 


/■>«V 

+ 4 (R - V e r sin B + /[l(VE[i) 


(Rf - Rf) 
' m n' 


+ 4 (R m - R n ) e r sin P 


dP 


_ = (? m ~ V _A_ dH 

dR " R^nR m /R n 4h 3 d0 


2R 


(R? - R 2 ) 
n m 


R/n(VV 


+ 4 


(R - R ) " 
n m 

e sin P 1 
r I 

R/n(R /R ) 
m n 


I 


(7) 


( 8 ) 


2. Full Film - Single Fluid Case 

If it is assumed that the same fluid is present on both sides of the seal 
and that the fluid film is complete, i.e„, no gas cavities, then, equation (8) 
becomes 


dP _ (P i " V + A dH 


dR R/n R. 

l 


4H 


3 dQ 


2R + 


(1 - R^) -1 

R/n R. 
x 


+ 4 e sin P 
r 


1 + 


(1 - R.) 


R/n R i 


The leakage coefficient for the radially eccentric seal is defined as 
^TT 

t 3 dP' 


:/ 


A H U' - H' 

r dR 


R dP 


(9)' 


( 10 ) 


Therefore, integrating at R = R = 1 we have 
-2n ° 

o (P. - P ) 
{-A He cos p - H / _• ° 

] r \ 


A dH 
4 d9 


f. , 1 - R ll , dH 

1 2 + / n R. “ A dP 


e sm 
r 


(con 1 t) 


141 


dH . . (1 - B i> 


ai “f 


(P. " P ) 
1 o 


/ n R ± 


d9 


:TT 


:tt 


JnR 


H d0 - A e 


(1 - R.) 


X O' 


r /n R 


x o- 


dH . Q . 0 
de SLn p de 


A 

4 


r* 9 ~ 

rl TT 



_2n 

1 ‘ R ? 

/ 1 A TJ\ 



[at 

J 

J 111 

OX 

de 

- A e 

V 

de [ H sin 13 d0 J 


If a wavy seal surface is assumed the film thickness is given by 
H = 1 + e cos n 9 

For this case the leakage coefficient becomes 

r2n 

(p - P ) I 

q =- j? a R ° / (He cos n e ) 3 de 


X O- 


- A e. 


(1 

7 


AhL- f 

n R. I 

1 rs-J 


2n 


(- e n sin n 9) sin (cr + 0) d0 

(1 - *,) 


_ (P. - P ) 

= — i — ° a 2rr (1 + lo5 e^) + A e e tt ~n — - - - -- cos oi 

/ n R. r n R, 


( 12 ) 


(13) 


for n = 1„ It will be noted that for n > 1 the hydrostatic leakage is unchanged 
but the pumping leakage becomes zero* 

The first term in equation (13) is the seal leakage due to the hydrostatic 
pressure differential* The second term is the leakage due to the radial 
eccentricity of the seal coupled with the circumferential variation in the film 
thickness* This second leakage term may either be inward or outward depending 
upon of, the phase angle between the surface waviness and the eccentricity 
vector (Fig. 2), or upon the direction of rotation, to, which appears in A. 

The shaft-seal radial eccentricity that will make the net leakage zero is 


e 


= -e . 
r 


r 

o 


(p i " P o ) h i (1 + 1 * 5 ® 2) 


3 u u) (r - r .) e cos a 
r o i 


(14) 


142 


(15) 

(16) 


For the minimum value of e in equation (14) we find that 

cos a = 1; o; = 0° 

(1 + 1.5 e 2 )/e = 2.45 e= 0.817 

Note that 0 < e < 1 and at € = 1, (1 + 1.5 e^)/e = 2.50. 

Substituting (15) and (16) into (14) and using the test seal conditions, 
i.e., r Q = 1.75 in, r ± = 1.125 in, \s> =4.95 x 10~ 5 psi/sec (350 centistokes) , 
the result is 

e = 2.52 x 10 5 (p i - p Q ) h? / N (rpm) (17) 

In order to estimate the order-of-magnitude of e let the following conditions 
be assumed, 

(P t “ P Q ) = 5 psi 

h j[ = 10" 3 in 

N = 10 3 rpm 

Therefore, for this numerical example 

e = 1.26 x 10~ 3 in (18) 


It would not be unusual to have a shaft-seal radial eccentricity of this 
magnitude. Even though the viscosity assumed in this calculation was quite high, 
the film thickness was much larger than normally would occur in a face seal. 
Further, if the equation for the seal leakage is written in dimensional terms, 
l . e . , 


2 rr (1 + 1.5 e 2 ) b 3 . (p. 

1 r X 

12 jZ’Q. R. 


P o ) 


+ 6 rr 0 ) r e 
o r 


(1 - R.) 

« 12/n R. h i C0S “ 


- Qo + Q 


pump 


then it will be noted that 


q oc h? 

O 1 


Q 0(5 h. 
pump i 


(19) 


( 20 ) 


Therefore, as the average seal film thickness, h^ decreases, the pumping component 
of the leakage becomes relatively more significant. 
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Thus it has been shown that for mechanical face seals the pumping leakage 
due to radial eccentricity can be a very important factor and for the small 
film thicknesses that normally occur only a small radial eccentricity is needed. 
In fact, it becomes a question whether or not this phenomenon can be completely 
eliminated from face seal experiments and applications, 

C* MISALIGNED-ECCENTRXC SEAL 

As will be explained later the experimental results were obtained for the 
misaligned-radially eccentric seal rather than the wavy-radially eccentric seal. 
Therefore, an analysis of this problem is of importance in verifying the seal 
hydrodynamic theory developed in this paper. 

The wavy and misaligned cases approach one another as the seal radius ratio, 
R^, approaches unity. The wavy solution is presented because it is much 
easier to obtain and the form of the final expression for leakage is simple and, 
thus, more readily understood. 

It was found to be more convenient to obtain solutions to the misaligned, 
full-film problem using numerical computer solutions. The FORTRAN IV program 
used on the GE 635 computer for this purpose is explained in more detail in 
Part IV. The approach used shows that the ,! short bearing 11 approximation is 
very good. 

For the misaligned seal the film thickness distribution is given by 
'H = 1 + Re cos 0 

and, therefore, is no longer a function of 9 

h b 3 g]+i le [h 3 H] = a 

- A 


“A 

The n short bearing 08 approximation could 
However, this approach requires using some series approximations. Since a 
numerical computer program was already available it was decided to use it 
instead. Equation (22) was solved in two parts. The first was for the hydrostatic 


( 21 ) 


only. Thus equation (1) becomes, 


:r ui ^ + U(? ^ 


r 3 r T u e se 

|(-Re r cos (3) (e t cos 0) 
+ (R + e r sin fO (- R © 


t sin 0)1 


Re e^. cos a - R e„ sin 0 
r t t 


[ 22 ) 


be used here to obtain a solution. 
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leakage where e r , the radial eccentricity ratio, was zero and the boundary con- 
ditions were P = P^ @ R = R^, P = P D @ R = 1. These solutions are plotted in 
Fig. 4 and take the following form, 


R i 

2 "<W 


Q q 12 Zn 

2 rt hJ(Pi-p o ) 


= 1 + a e 


where the slope, a, is a function of R^ and is shown in Fig. 5. 

The second part of the solution to equation (22) is for the pumping 
component and is obtained by using a finite value for and applying the 
boundary conditions P^ = P q = 0. As can be seen from equation (22) the result 
will be in terms of the constant factor A e r e t cos a, the same as for the wavy 
case. Values for 


A e e cos a ~ 2 

r t e e cos a h.u>r 

r t i o 


- f (\> e t > 


are plotted in Fig. 10, 

D. ECCENTRIC SEAL WITH CAVITY 

Once the equation (8) for the radial pressure gradient has been formulated, 
two equations can be determined for the cavity boundary shape from the cavity 
boundary conditions. Using the short bearing approximation, the upstream 
cavity boundary condition is 


<? P >c ■ III - 0 


and for the downstream cavity boundary condition 

®2 ■ V A [r ■ D J i c ] = H 2 3 [r2 


sp.&p de 

39 " 5R dR 


[ dp de 1 

dR dR 
c c J 


The circumferential pressure gradient, dP/S0 , has been neglected as in the 
original assumption. 
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The equation for the upstream cavity boundary shape is found by substituting 
equation (8) into equation (25), 


H" 


dH 

d0 


(P. 


m 


P ) 

n 


A 

4 


2R^/n R /R + (R 2 - R 2 ) + 4 e sin ph /n R /R + (R - 
c nm mn r Lc n m m 


The equation describing the downstream cavity boundary shape is obtained by 
substituting equation (8) into equation (27), 



|dR\ H 2 R c 


jdP) 

»d0' c (% - 

h 2 ) a 

1 dRl 



/ 

H 3 R 


P 

2 c 

J 

m 


R U’ 
c r 


U 0 


. n , A ^2 

' (H 1 - V U e A 1 R c/ n V R n 4H 3 2 d9 


2 2 

(R n - R J 
n m 

c ' R R /R 
c m n 


2R + 


t 4 e sin p I 1 + 


(R n " R J 

n m 

R Ja R /R 
c m n 


R U r 
c r 


u e 


(29) 


It should be noted that when the appropriate substitutions are made for a, 
and m, for the various regions discussed before, will be the inner cavity 
boundary, i,e», the boundary closest to R^, and R£ will be the outer cavity 
boundary, i.e. , closest to R o . This nomenclature is not to be confused with 
H^, which is the film thickness on the upstream boundary of the cavity and H 2 
the film thickness on the downstream boundary (Fig. 3). By definition, the 
upstream cavity boundary is where the fluid streamers leave the complete film 
region and enter the cavity region, and the downstream cavity boundary is where 
these same streamers re-enter the full film region. 

The point where the cavity boundary is tangent to the relative surface 
velocity vector distinguishes the difference between the upstream and the down- 
stream cavity boundaries. This point, therefore, is located where the following 
expression is satisfied, 



(30) 


Note that equations (29) and (30) are consistent since at this point of inter- 
section the upstream boundary condition, equation (25), also applies and equation 
(29) reduces to equation (30). 
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The remaining required expression is the one needed to solve for H^. The 
which corresponds to a given H 2 in equation (29), i.e., to a specific point 
on the downstream boundary, is found by tracing the streamer path through the 
cavity. That is, the proper for a given H 2 is obtained by determining the 
departure point on the upstream boundary for the fluid streamer entering the 
downstream boundary at H 2 (see Fig. 3). For a given point on the downstream 
boundary of the cavity (R d , 0<j), the required coordinates on the upstream 
boundary (R u , 9 U ) are obtained from the following equation 


R + 2 e R sin 0 
u r u u 


2 

R , + e R, sin B 

d r d ( 


(31) 


These equations are sufficient to calculate the dimensions of the gas 
cavity and thus yield the solution for pressure in the seal film. From this 
solution the load, flow and torque coefficients may be calculated, 



(32) 


(33) 

(34) 


Note that equation (34) represents the torque about the center of the seal. 


E. EXPERIMENTAL RESULTS 

The experimental apparatus used is described in Appendix C. (Figures 
referenced in this Part E are in Appendix C.) Figure 21 shows the comparison 
between theory and experiment for the full film outward pumping case. Fig. 22 
shows the comparison for the inward pumping case. In the latter set of data it 
was difficult to supply enough oil to the reservoir surrounding the outside of 
the seal because of the net inward pumping taking place. Therefore, only low 
speed data was obtained. 

When the outside of the seal is not flooded, gas cavities appear in the seal 
interface as shown in Figures 9, 10, 11 and 12. When this occurs the inward 
pumping effect is diminished as can be seen by the data plotted in Figures 19, 

20, 21 and 23. 
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F. DISCUSSION 

1 , Correlation of Theory 

For the full film case the correlation between experiment and theory was 
quite good. As yet numerical results have not been obtained from the theory 
developed for the radially eccentric seal containing gas cavities. The equations 
for the full film - two fluid problem have been derived but not set up for 
solution as yet. It is hoped that these two solutions can be pursued further. 


2 . Previous Investigations 

The work of Denny £l j demonstrated the existence of hydrodynamic forces 
mechanical face seals. He and Nau 31 both observed the "inward pumping" 

»> * r n v J 


m 


Nau 


M - M also observed gas cavities in transparent seal 


phenomenon, 

experiments., Because he observed these two phenomena simultaneously he was led 
to the conclusion that the cavitation produced the pumping leakage, Lyman and 


Saibel 




in their approximate cavitation seal analysis, tried to show 


theoretically that gas cavities did cause inward pumping in addition to hydro- 
dynamic lift forces. Orcutt and Cheng |l 2 ^ and Pape [13] also used approximate 
cavity solutions to demonstrate the relationship between this phenomenon op the 
hydrodynamic lifting force observed by Denny |\j . 

All these solutions J^ 5 j , [12] , [13J, neglected to use conservation of flow 
boundary conditions (Swif t-Floberg) at the edges of the cavities. Thus continuity 
of flow was not maintained. Pape JjL 3 ^ recognized this fact and Orcutt and Cheng 
[12] did not use their solution to calculate flow rates, for if they had, the 
same leakage rate would not have been obtained at both the inner and outer 
radii. This was the mistake that Lyman and Saibel £sj made. They integrated 
the local, radial flow rates at the inside radius of the seal and obtained an 
inward leakage when cavities were assumed. If, however, they had done the same 
at the outer radius they would have obtained an outward leakage for the same 
case. Thus a basic inconsistency existed in their solution and they were mislead 
to believe that cavitation produces inward pumping in face seals. 


III. CAVITY EXPERIMENTS 
A. WAVY SEAL 

Although the cavity experiments were carried out on a misaligned seal', . 
the analysis for a wavy seal interface that was carried out in Phase I 8 j 


can 
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be used as an approximation,, Thus, the film thickness variation for a mis- 
aligned seal is 

H = 1 + e R cos 0 (35) 

but can be approximated by 


H = 1 + e cos 9 


(36) 


for R 1, i.e. s small land width seal. For this case the upstream cavity 
boundary shape is given by the following expression 


Y 


H 


dH 

d0 


(P 


4 [ 2 ^ 


- V 

R. 


In ^ + R 2 
R c 
c 



(37) 


where 

1 dH _ e sin 9 
^3 d0 (1 Te cos 9)^ 


(38) 


Equation (38) is plotted in Figure 6 for various values of e, In addition, the 
curve which represents the limit of application of equation (37) is shown. This 
occurs at 



(39) 


which yields the following expression 

-cos 9 

g = - — — «— - ■ 

(3 sin^ 9 + cos^ 9) 


(40) 


Values of Y corresponding to equations (39) and (40) are designated ^crit* 

In order to find the seal conditions for which a cavity will first appear, 
R c should be set equal to R Q in equation (37). Therefore, from equation (37) we 
have 



13.57 

Y 

cnt 


(41) 
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where 


R. 

i 



2.250" 

3.50" 


0.643 


for the test seal. Using a 350 centistoke (4.95 x 1C)" 5 psi-sec) (silicone oil) 
the speed-pressure ratio as a function of film thickness and misalignment can 
be computed from equation (41) as follows, 


N/p. 


= 1.426 x 10^ 



Y 

cnt 


, rpm/psig 


(42) 


Taking the values of from Figure 6, the curves in Figure 7 were plotted. 

The parameters bu and e are related to the minimum and maximum film thick- 
ness by the following expressions 


h . = h. (1-e) 

mm i 


h = h . (1+e) 
max i 7 


h - h . 

- max min 

6 " h + h . 

max mm 


— h 4* h . 

h = max min 
i 2 


(43) 

(44) 

(45) 

(46) 


In the seal experiment it is the minimum and maximum film thicknesses that are 
set. 

Once the cavity has formed and is of a finite size, the position of the 
cavity nose can be found by calculating Y from the following expression 

Y = 1,426 x 10 5 tu Pi (47) 

~N 

and then finding 8 from Figure 6 using the correct e curve. 


B. MISALIGNED SEAL 

In the misaligned or tilted face seal the film thickness is a function of 
radius as well as 8. For this case the integration of the "short bearing" 
approximation to Reynold's equation is somewhat more complicated. 
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For the general case of seal misalignment and surface waviness the film 
thickness distribution is 

H - (1 + e cos n 8) + R fe t cos (0 + t|f)j ( 

= a + b R ( 


= -e N sin n0 - R e t sin (0 + ty) 


H = e fc cos (0 + ty) = b . (f 

The "short bearing" approximation to Reynold's equation for an eccentric face 
seal is 


1 S , ot 3 SPv_ '[ a dH x i r . ft l dH 

R 3R (EH 3R ~ A ( '*r cos P 3R + I 1 + R l “ S 30 


h o " 5 


* 

)= - A sin 0 e n sin n0j 


+ £e r cos 0 • b + e n sin n8 + 8^ sin P * sin (0 + \jf)J 


+ I e sin (9 + \|i ) 


]r 1 


= -A) A(0) + B (9) • R + C (0) • R 


Integrating this expression we have 


SR H 2H 3H RH 


P » -A ) A 




R 2 dR 


f dR_ 

J RH 3 


+ 5(0) ~ (0 ) 


F (R, 0) l+«(0) • G (R, 0) + t (9) • 
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Applying the boundary conditions 


P - P. @ R - R. 
1 1 


P = P c @ R = (at the cavity) 
the unknown functions §(0) and i|f(0) can be determined, i.e„, 


1(0) = 


_ (P. - P e ) -A [f(R c ,0) - R(R.,9)] 


G(R ,0) - G(R ,0) 
1 c 


(57) 


(58) 


♦ (0) = P c + A F(R c ,0) - §(0) • G(R ,0) • 


(59) 


Substituting equations (58) and (59) into (54) and (56) 
(P - P c ) = A [f(R c ,0) - F(R,0)] 


dP A_ 

dR - H 3 < 


A + BR + cr 

A + 2 + 3 


(60) 


+ j< P i - V ~ A [ F < R c’ 9 > “ F <V 9 >]j |^“ 


G(R,0) - G(R c ,0) 
0) - G(R c ,0) 


(P t - P c ) -A [f(R c> 0) - F(R.,0)] 
RH 3 [g(R., 0) - G(R c ,0)] 


(61) 


The integration indicated results in a singularity at the value of b = 0. 
Therefore, it was necessary to express the integral of equation (55) in the 
form of infinite series expansions. The functions of integration become: 

F(R,0) Z (-l) n+1 (|| n_1 R n + -4r 7 ( ' 1)n+1 f(l! R 

J n-1 2 la I 3 n^l 

a 2a 


n+1 


op n+1 , .-v , n-1 n+2 

3a 


G(R,0) = 


‘l_ / R . 1 

/-n — + — =■ 


La 


H 


a 2 H 


2aH 


(62) 
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where 


H - a + bR 
a = (1+e cos n9) 

b = j- e t cos (9 + ty)] 

A = [ e r s inpe n sin hb] 

= [b e r cos p + en sin n0 + ® r e t sin P sin (9 + ♦>] 
C = |e fc sin (9 + \|f)j 


B 


Since the upstream cavity boundary condition is 


|i£ ] . 

dR I 


(63) 


equation (61) is set equal to zero and is set equal to R q = 1 since this is 
where the cavity will first appear for P c = P . Thus, we have the following 
expression 



If the reference pressure in A is taken 
concentric seal (e = 0) then, equation 


(G. - G ) + (F - F.) (64) 

' i o o 1/ . 

as (p. - p ), then P. - P = 1. For the 
r i r c ’ l c 

(64) becomes 


1 

A 



(G. - G ) + (F 
l o o 


V 



' R i (l+b) ' 

+ 1 1 . 1 

1 

1 

1+bR. 

u i J 

1+bR. " 1+b 2 

l 

(1+bR.) 2 

(1+b) 2 


V 


+ E(-l ) n+1 

n=l 


& o "’ 1 * V* 


(65) 


/ 

The series in equation (65) does not converge as fast as desirable. Therefore, 
the following technique was used* First it was noted that 


lim <-l> n+l 


n(n+l) , n-1 / . „n+2l n i, i n-1 

20*2) b I 1 "* = 2 M 


( 66 ) 
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and 


(- 1 ) 


n+1 


n(n+l) ,n-l 
2(n+2) b 



for 0 < R. <1 and b <.0. 

l 


Therefore, 


(67) 


OO n+1 

S(-l) +1 

n=l 


n(n+l) ,n-l 
2(n+2) D 



2(1 b | ) z 


( 68 ) 


The series 



(69) 


converges faster than the series representing (F q - F^) . Also the expression 

1 


2(1" | bj) 


(70) 


gives an upper bound to the series. This information yields a series representa- 
tion of (F q - F^) which converges more quickly than the original one in equation 
(65). This new representation is 



1 

2(1- | b |) 2 


lim 

n ^ 


N 

L 

n=l 


n 

2 



N 

- 2(-D 


n+1 


n=l 


n(n+l) , n-1 /. „ n+2 

2(n+2) b 1 - R i , 


(71) 


The listing of the Time Shared Computer program which calculates equations 
(65) and (71) is given in Table I. The numerical results are shown in Figs. 8 
and 9. The line labeled 1/A cr ^ t represents the conditions for which the cavity 
will first appear. The other curves labeled e (= 0.1 for example) give the 
location of the nose of the cavity. The two figures represent conditions for 
two different test seals. The first one is for 


R = II - „ 2.25 

i ro do 3.50 


0.643 . 
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The second is for 

R. = m| = 0.83 * 

1 Z./5 

These are the conditions for the two transparent test seals used in the seal 
experiments. 

IV. SEAL PROGRAM (GE 635) IMPROVEMENTS 

Although it is nice to have closed form analytical solutions to the various 
face seal problems, it sometimes becomes necessary to resort to numerical 
computer solutions. Once the basic problem has been set up on the computer 
it will usually handle various seal film geometries equally well, whereas the 
analytical solutions hold for more limited conditions. Therefore, the GE 635 
SEAL PROGRAM has proven very useful. However, the computer solutions to the gas 
cavity problems require a large number of iterations and it is desirable to 
improve the numerical approach so as to reduce the computer time. 

A method has been developed that not only conserves computer time but also 
gives a numerical solution to the n short bearing" approximation as the first 
iteration in the process. This method is a variation of a method described in 
Reference 16. In Reference 16 the film pressures in a bearing are obtained by 
direct solution of a finite set of simultaneous equations. This approach still 
requires a large number of computations and in some cases may not reduce the 
computer time over that required by an iterative solution. 

It has been recognized that in the seal problem the influence of the cir- 
cumferential pressure gradient is much less than that of the radial pressure 
gradient. Thus the "short bearing" approximation, wherein SP/S0 is neglected 
altogether, has resulted in very good solutions to the face seal problem. 

The new method of solution used in the GE 635 SEAL PROGRAM makes use of the 
fact that a very good approximation can be obtained by neglecting dP/d0, and 
solves the "short bearing" approximation to Reynolds 1 equation numerically. 
However, this numerical solution is obtained by direct calculation of the 
simultaneous equations rather than by iteration. This is a very quick computa- 
tion since the number of simultaneous equations is equal only to the number of 
radial grid points taken in the numerical approximation. Once the "short 
bearing" solution is calculated and printed out, the SP/30 terms in Reynolds 1 
equation are included as correction terms. Several corrective calculations are 
then made in order to obtain the two dimensional solution to Reynolds 1 equation. 
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The development of this new method of numerical solution starts from the 
numerical formulation of equation (1) 


D.P. . = 
J 


1 

3 

3 H 1 


P P 

- ? 
2A0 

i+l,j + i-l,j J 4 A 0 Hij | 

1 30 | 

lij 

. i+l,j i-lJj 




2&R 


h >j+ i + p m-iJ 


3H_ +r 
Hij dR.. + j 


AR' 


U'H 


3 \ SR <U r H) + | R 
2H ij 


+ 57 


Dj 


[ —1 + - 4 ] 

L A0 J AR J 


ij 


. 2^m 


P -P 

i.J+l i 


H = (1 + a % + b| ) + ® cos n9 + e R cos (0 + ty) + e * 

t s 




These can be combined as follows 


P, , = A. . P. . + A.. P. _ . + A P. . .. + A. . P. - A /ri , 

ij ij l+l, j 2j l-l, J 3j i,j+l 4j i,j-l 5 j (72) 


where 


A ij 


A„ . = 


2j 


A. . 
3j 


A, . 
4j 


A . 
5j 


D. 

J 


[(1/2A9 2 ) + (3/4A9 H^) OH/30) ] / D 

[(1/2^0 2 ) - (3/4 A 0 H ) OH/30) ]/ 1 D 

|(R?/2AR 2 ) + (1/4 AR) [(3Rj/Hij) OH/3R).. + R^j / D. 
|(R?/2AR 2 ) - (1/4AR) [(3R^/H tj ) (3H/3R) ij + Rj]j / D. 

(AR 2 /2H.j) ||5 (U*H) + (U'H/R) + |g- (U^H)J / D. 


1/A 0 2 + R 2 /AR 2 


(73) 

(74) 

(75) 

(76) 

(77) 

(78) 
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If it is assumed that the pressures in the circumferential direction 
P P 

i+l,j i-l,j 

are known or assumed, then one may start at one of the boundaries of the seal and 
work across the seal obtaining expressions for P^j in terms of the starting point 
boundary pressure and P^ Thus proceeding across the seal to the other 

boundary pressure the last j + ^ is known, i.e., it is the second boundary pres- 

sure, and working back across the seal all P^j's are determined on that radial 
line. Going on to the next radial line the procedure is repeated. 

The equations used for this procedure are 

P U ' X J P i,j+1, + 

where 

x j ■ V U • % X J-1> (J - 2 ’ 

X x = 0 

Y j = (A lj P i+l,j + A 2j p i-l,j ” A 5 j + A 4j Y j-l )/(1 ' A 4j X 
= P^ (first boundary pressure) 

P i jmax = P n ( secon< * boundary pressure) (84) 


(79) 

(80) 

(81) 

j-l ); = 2 »*”» j max^ 

(82) 

(83) 


The calculational procedure is as follows, 
Step 1: Short Bearing Approximation, 


i.e 


• > 



(a) 

A lj = A 2j = 0 


(b) 

D. = R2/AR 2 


(c) 

Calculate Y^, (j = 2,3,. . . 

. , j -1) ; i fixed 

(d) 

Calculate P. . (j - j , , 

ij, VJ J max-1 

n, 

. , 3,2) ; i fixed 

(e) 

Repeat for all i’s 


2: 

First Correction 


(a) 

Calculate A^ , A^^ from equations 

(73) and (74) 

(b) 

Use P . . - . and P. , . from Step 1 

i+l, J 


(c) 

Calculate X ^ , Y ^ , (j = 2,3, . . . 

. , j , ) : i fixed 

’ J max-1 * 
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, 3 , 2) ; i fixed 


(d) Calculate P^ , (Eq„ 79), (J - j max _^ # • • • • • 

(e) Repeat for all i’s 

Steps 3 to N: Additional Corrections 

Repeat procedure under Step 2 always using the j and j values 

from the previous step (or iteration) . 

This method is now in use on the GE 635 SEAL PROGRAM and has proven to be 
very successful. The results obtained thus far indicate that for the full film 
seal the "short bearing" solution to Reynolds’ equation is very close to the 
two dimensional solution. 


V. TWO FLUID STUDY 

The analysis of this problem has not been completed. Some of the 
basic relationships have been derived and the problem set up but not solved. 
Therefore, only a summary of the basic equations and approach to the solution 
will be given here. 

The first step to this solution should be to take the single fluid, no 
cavity, misaligned seal solution and determine the regions in the seal inter- 
face occupied by the fluid coming from inside the seal and outside the seal. 
From this a general picture of the fluid regions will be obtained. This is- 
needed for a better insight into setting up this problem. 

From considering this problem qualitatively it appears that there may be 
multilayers of fluid in the seal. It looks quite likely that one fluid may be 
next to the two seal surfaces and the second fluid between these two strato. 

In the general case considered here only two layers of fluid have been con- 
sidered in the boundary conditions. 

The film cross-section geometry is shown in Fig. 12. For this case the 
problem is summarized as follows; 

Two Fluid Problem in Polar Coordinates 

Velocities 


u 


rl 


1 

2 3 r 


z (z-hp + U 





_ 1 dp 

u r2 2p» 2 3r 


z (z-h^) - h (z-h^) 


+ U, 


hr 


(z ' h l ) , (h-z) 

(h-h x ) U ri (h-t^) 
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91 2|x 1 r 59 


z (z-h!) + U o0 ( 1-^ ) + %i 


_ 1 5£ f, ■ . w , J . „ (z "V , (h-z) 

U 9 2 ~ 2(J> 2 r 59 [ (z_h ^ ( z ' h l)j + u h e (h-h^ + U 9i (h-h^ 


Flows 


h n * h n h. 

' ‘ 12^ I? + “or 2 1 + u ri T 

- (h ~ h l ) 3 3 d , „ (h - h l ) . 

1 2 n» 2 5r hr 2 1 


(h-hi> 


ri 2 


h l 1 5p , TT h l h. 

12^ r 59 + U o9 2 + ^Bi ^ 


< h - h l )3 lan <^1> ( h - h !> 

12^ r 59 U h9 2 + u 9i 2 


Continuity 


l_ , x ' 

5r ^ rq r^ + 99 


L. rh 3 3£ + 1S_ h 3 5£ =6 i5_ 

5r rh l 5r r 59 h l 59 5 ^1 \ 5r 


59 [ h l 


< U o9 + u 9i> 


h [ r(h_h i 


rh. (U + u .) 
1 or rx 


^3 ,15 

; 5r + r 59 


+ m <h ' h i )(u h0 +u ei ) 


<h - h i >3 m " 6i m! 7 r(h - h i )( V +u ri 


Shear Stress 


T . = U»- 
n r l 


3z I . ^2 \ 3 z 
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T 0i = ^1 Ji = ^2 \—] i 


3u . 
rl 

1 

dz 

2^1 

i auri i 

; 

\ dz 

i 

3 “r2 _ 

1 

dz 

2^2 

3u r2\ 

(h 

dz It 


6u ei _ 

1 


. U u . 

0 £ v N or n 

S? < 2z - h i> - -hi + -q 


(u .-U ) 


*1 ar + h i 


[22-Ch+hp] 


hr 

(h-hi> 


, -h) a (U, -u .) 
1 &£ hr rr 


(h-hl) 


^00 . U 0 i 


<2 - h l> - Ef + M 


Su 0i \ = i d£ 

3z 2^ r 30 h. 


8u 92 i lip r. ,, „ ,i , ^ D h9"'W 

ri£ L 2 z -< h+h l>J + (h-hj) 


(h.-h) 


(^UQ _U Q ,• ) 


6 2 1 = H2' 1 3p v h0 0r 

. . ’ oz ^ 2^2 r 30 (h-hj} 

Shear Stress Boundary Conditions 

h n . (u .-U ) (h.-h) .. (U, -u 

1 3p , ri or ; 1 3p , hr r: 

2 3r +p, l h 2 3r + ^ 2 (h-h^ 
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or 


hS£ , fil , ^2 

2 5r U ri (h-h^ 


. ^ = 0 

h t or (h-h x ) 


or 


u . 
n 


k 4. IT 4- ^ 2Uhr 

2 3r h 1 U ° r (h-tn ) 


JjL, , 2 

h l ^ h-h i^ 


Jhia£. (u 8i~ P o9 ) . (h l' h) 15£ ^"“e^ 

2 rS9 + ^l h x 2 r de + ^2 (h-h^ 


or 


— + u 

2r ae ei 


_I + _X 

h t h-h 1 


^l U o6 ^2 U h8 n 
h t " h-h 1 


or 


U 8i ” 


h_ S£ , ^1 P q 8 . ^ 2 P h8 
2r d8 h i (h-h i) 


*1 , **2 
+ (h-h^) 


Non-dimensional Form of Equations 


r = — p = -2— H = — U* = ^ 

R r ,P p * H h . ’ ru) 

o *o i o 


d_ 

SR 


RH 3 ai 


+ IS- 
r ae 


H 3 

H i ae 


‘MS 


RH. (U ' + u !) 
1 or n 
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H i (u o^ + V 


a_ 

aR 


e(h -v 3 H 


+ U- 
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(H-Hi> (U^+UqP 
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I 




1*1 o0. . ^2 _hr 

' R 39 H 1 + (H-H 1 > 



+ 


(Si) 

(h-h")' 


For the case where 


U = 0 and 
o 


U ' = - e cos P 
hr r 

« R + e r sin p , 


the seal equations become, 


L. UmI + IL L 3 

3R pi 3r] + R 39 [ 1 39 J 


= A 


1 3R 


. RH, 


- 3H 3P 


'|J- 2 e r cos p i 


A x bR V> 1 (H-H 1 ) J 


1 ^ 2 ^ 1 ^ 1 
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H T (H-H 1 ) j 


H, 


+ A 1 39 


- 3 H 13 P ^ ( R +e r smP) 

A x ‘ R 36 jJ- 1 


(H-H x ) 


1 _ + 


[H x T (H-H 1 ) 


3_ 

3R 
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+ I&- 
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^Kh ^ R(H -V 


+ V 


1 39 
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3H 3P ^ 2 e r cos j 
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First, it can be seen that a new parameter has been introduced into the set 
of simultaneous equations that must be solved. As might have been expected, 
this new parameter is the viscosity ratio, p2/pi« This parameter has been used 
in correlating the two fluid data in Appendix C. 

Second, it should be obvious that the solution to P is the same for both 
differential equations, That is, the last two equations given are to be solved 
simultaneously for P and H]^, the two unknowns. 

VI, CONCLUSIONS 

1. The short bearing cavity solution presented in has been confirmed 
experimentally, 

2. The pumping theory developed in this Appendix has been confirmed 
ex p er ime n t a 1 1 y , 

3. In concentric face seals of fixed geometry, cavitation produces a net 
hydrodynamic lifting force and moment, but does not alter the leakage 
rate, i 0 e«, it does not produce an "inward pumping" leakage. 

4. Pumping leakage is produced by a combination of radial eccentricity 
and seal misalignment (tilt) or one wave waviness. 

5. The above pumping effect can be either inward or outward depending, 
upon the relationship between the radial eccentricity and misalignment, 
and the direction of rotation, 

6. As long as a complete fluid film exists in the seal interface, the 
hydrostatic pressure differential does not affect the pumping leakage 
and produces a leakage which is independent of speed and, of course, 
always in the direction of decreasing pressure. 

7. The existence of gas cavities in the film reduces the pumping effect, 
in eccentric seals, rather than causing it as heretofore postuated. 

VII. FUTURE WORK 

An experimental program has been completed to verify the predicted eflects 
of fluid inertia jio], £l.l], and to study the influence of having fluids of 
different viscosities on either side of the seal. The solution to the two fluid 
problem and the single fluid pumping problem with cavities should be completed 
and included in the digital SEAL PROGRAM. 
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SECTION A-A 
SEAL INTERFACE 



Figure 1 

Mechanical Face Seal Geometry 


FIXED SEAL, ROTATING SHAFT 



SEAL COORDINATES ARE r,0 

FILM THICKNESS 

i u r o 

U r = — = - € r cos 

H = 1 + € cos 0 

r 0 w 

Uo = 6 = R + € r sin ft 

0< € < | 


r is measured from the center of the seal 
p is measured from the center of rotation 


Figure 2 
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UPSTREAM BOUNDARY CONDITION DOWNSTREAM BOUNDARY CONDITION 
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Cavity Boundary Conditions 




FLOW COEFFICIENT, -q InRj /2 tt 


FIGURE 4 

HYDROSTATIC FLOW COEFFICIENT 
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sin 


FIGURE 6 



0 , Degrees 
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SEAL CONDITIONS FOR INITIATION OF CAVITY 
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SEAL SPEED- PRESSURE RATIO 
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PUMPING LEAKAGE '-<lp ump /A < r « cos a 


Figure 10 

SEAL PUMPING LEAKAGE 
DUE TO RADIAL ECCENTRICITY 
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FIGURE 11 
SEAL LEAKAGE 



e 
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Figure !2 

FLUID FILM GEOMETRY 
FOR THE 

TWO FLUID PROBLEM 
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FIGURE 13 


EFFECT OF RADIAL MISALIGNMENT 
ON SEAL LEAKAGE 
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.6 
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LEAKAGE 
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ESTIMATED LEAKAGE 



-6 



-.8 


I 


LEAKAGE 



RADIAL 

MISALIGNMENT 


NOMENCLATURE 


A - seal area, in^ 

2 

- flow area, in 

h - film thickness, in 

- average film thickness at inside radius, in 

e - displacement between seal center and shaft center, in 
y - coordinate across the film, in 
p - coordinate measured from center of shaft, in 
p - density, lb/in^ 

r - coordinate measured from center of seal, in 
9 - angular coordinate, radian 

P - a + 0, radian 

a - phase angle between waviness and eccentricity, radian 

V - angle measured to p, radian 

i|r - tilt phase angle, radian 

t - time, seconds 

U - fluid velocity, in/sec 

(V - rotational speed, radian/sec 

Q - flow rate, in^/sec 

p - pressure, psi 

t - shear stress, psi 

W - seal load capacity, lbf 

SM - shaft friction moment, in-lb f 

Cp - specific heat, Btu/lb°F 

k - thermal conductivity, Btu/sec in°F 

T - temperature, °F 

4 - dissipation function, sec 
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dynamic viscosity, psi-sec 


Dimensionless Terms 
R - radius ratio = r/r 

o 

H - film thickness ratio = h/h 

P - pressure ratio = p/p 


, r 2 

seal number = (r~) 

P„ h ,* 


e - waviness parameter 

e - radial eccentricity 


tilt eccentricity 


load coefficient = 


W 


P r 
r o o 


friction coefficient = 


SM 


u 3 

h . p r /6 
i r o o 


flow coefficient = 


1 2 fu. Q 
h 3 P 


U' 


pumping velocity variable = — 

o 


radius variable = 


r - r. 

i 

r -r. 
o i 


film thickness variable = 


h - h, 
o 1 

h 




D. 


number of circumferential waves 
numerical grouping 

numerical grouping 


Subscripts 

0 - ; at the outside radius of the seal 

1 - at the inside radius of the seal 

C - within the cavitated film 
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TABLE I 


TIME SHARING COMPUTER PROGRAM 
FOR CALCULATING 1/A (EQUATIONS 65 AND 71) 


10 DIM T(50) ,1 (50) 

15 READ NO 
20 FOR 1=1 TO NO 
30 READ 1(1) 

35 LET T(I ) = Id) * 3.14139/180 
40 NEXT I 
45 LET P=0 
50 READ R 1 ,E 

53 PRINT "TILT ECCENTRICITY =",E 

59 PRINT 

60 PRINT "THETA", "1/L(FN)","1/L(FT)","N", "CONSTANT TERM" 

70 FOR J= 1 TO NO 

72 LET T0=0 

73 LET A =0 

74 LET L0=0 

76 LET B= E* CiS(T(J)) 

79 LET C= E*SINCTCJ)) 

81 LET A 1=0 
85 LET B1 =0 
90 LET Y= R 1 

92 LET L= 1/2*1 /<l-ABS(B))f2 
110 PRINT 

200 LET Xl= (L06 ( Y* ( HB) /< 1*B*Y) )>+< 1/(14B*Y)-1 /Cl+B) ) 

210 LET X2 = .5*((l/(l4e*Y)t2)-(l/<l+B)t2)) 

220 LET X3= X1 + X2 
225 FOR N= 1 TO 100 

230 LET T=(-l)t<N+l)*<(N+l)/<N4£))*N/2*Bt<N-l)*(l-Yt<N+2)) 

235 LET T0= TO+T 

236 LET Ll= l/2*<N*<ABS(B)t(N-l))) 

237 LET L0= L04L1 
240 LET D= L-CLO-TO) 

242 IF ABS ( (D-TO) /D) < .01 THEN 300 

260 NEXT N 

261 PRINT " CONDITIONS NOT MET IN 100 TERMS" 

263 60 TO 320 

300 LET X4 = (X3+T0)*C/3 
302 LET F = (X3+D)*C/3 
304 PRINT I(J) ,F,X4,N,X3 
320 NEXT J 
900 DATA 27 

910 DATA 180,182,184,186,188,190,192,194,196,198,200,202,204,206,208 

911 DATA 210,212,214,216,218,220,222,224,226,228,230,232 
920 DATA .643,. 9 

99999 END 
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APPENDIX E 
SEAL DESIGN II 
BY 

H.J. SNECK 

I. INTRODUCTION 

Included in this section are several analyses which are directed at furthering 
our understanding of mechanical face seals operating under hydrodynamic conditions 
with a Newtonian lubricant. 

The first, study examines the effects of heat generation in the face seal due 
to viscous shearing of the fluid. One of the purposes of this study is to evaluate 
the constant viscosity assumption made in our previous analyses. That is, how 
much does the fluid temperature (and thus the fluid viscosity) .vary across the 
seal face. As one might expect from the fact that leakage rates through the 
seals are low, the conclusion of this study was that the fluid in the seal inter- 
face film is almost isothermal. Thus the assumption of constant viscosity is a 
valid one for most seal problems. However, each seal application should be 
evaluated individually. 

The second study examines the interaction of the fluids on either side of a 
face seal due to inertia effects. That is, the fluids being separated by the 
seal may come in contact in the seal film and inertia forces may cause a reverse 
flow through the seal against the hydrostatic pressure drop. The analysis solves 
for the streamline separating the two fluids and determines the conditions for 
which this streamline will span the width of the seal and result in an interchange 
of fluid across the seal. 

II. THERMAL EFFECTS WITHIN AN ALIGNED FACE SEAL - LAMINAR FLOW 

A. INTRODUCTION 

In the previous analyses of face seal performance (Ref. 1) it was always 
assumed that the flow through the seal was isothermal. Since the density and 
especially the viscosity of the sealant fluid are reasonably sensitive to tempera- 
ture changes, it is important to establish the accuracy of this assumption. 

B. ANALYSIS OF LAMINAR FLOW 

Because of its simplicity, the aligned, laminar, flat seal will be 
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analyzed and the results of the analysis extended to seals in general. 

It was shown in Ref 1 that the radial and tangential velocity profiles 
in laminar flow are given by the equations 


u = (z 4 . ze 3 ) 

( * I2p, } 

V - r«»(f) 


(1) 

(2) 


where it is assumed that p and p. are independent of Z. 
The pressure distribution is given by the formula 

p * h + (r W> - “ « 


(3) 


Equations (1), (2) and (3) will be particularly useful in evaluating the 
various terms of the energy equation which is 


tz ( - k §> + e% a 0 " u + “ ' 


v 2 . 2* 


(4) 


Since the liquid film is quite thin it seems reasonable to assume that the 
temperature variation across the film in the axial (Z) direction is relatively 
small. This is consistant with the previous assumption that p. and p are not 
functions of Z. With this assumption the energy equation may be integrated 
across the fluid film in the axial direction. 

The first term in the energy equation (4) represents the axial conduction 
of thermal energy, which when integrated between the seal surfaces yields the 
total local heat transfer flux exchanged between the fluid and the surroundings. 
That is, c 

/ h ( - k S> dz - - (k §> - <- k i> < 5 > 

o CO 

II 

= q = total local heat flux out of fluid 

w 

The second and third terms (associated with the radial transport of energy) 
take the form when integrated in the axial direction 

(p c P S • f?> f* iz (6 > 


183 


where 


/ 


udz = -r^-" 

27 r 


( 7 ) 


Differentiating the pressure equation (3) and rearranging yields 


J 


udz 


27 r 


'C 

l^i 


3pr to' 
10 


. $2- 
S r 


( 8 ) 


for the integral coefficient in equation (6) . 

The last term in equation (4) represents the energy converted by viscous 
dissipation within the fluid film. From the radial velocity distribution given 
by equation (1) the dissipation associated with the radial flow is 

2 c£_ . 2 . Sp. 1 d 3 . or a) 2 

j oz lm &r 20u, (pr “ 5 v + ’ ( ' 

o 


& - ■ w <£> : 


7 


-) 


(9) 


The tangential flow contribution to the dissipation is 

n (f^ dz -feiiJil! 

dz c 




(10) 


With the aid of the foregoing evaluations, the energy equation integrated 
across the fluid film now takes the form 


( Q n ST = d dp. r 3 or _ Spl 

p ^27r Sr Sr [ 10 SvJ 

3 v 2 ,3 „ v , 3 ... 2 2 


K + pC P W S7 

+ £_ A 2 . £_ (or ^ 2 ) & + 1 + h, ( r (11 ) 

+ l^Jb ''Sr' 2Qi lp ^ Sr 7 |i ' 4 ' ^ c ' 


or, after some cancellation and rearranging 

.3 


a" + r ( Q 7 M. = 
% p( ^p ^2nr^ Sr 


3 3 2 2 

c . 2. Sp . 1 £ ,orii) v u(ru)) 

W (p m } t + 7 r ( 4 > « (1: 2> 


In order to integrate equation (12) in the radial direction , it will be 
assumed for the time being that p and p are independent of the radial coordinate. 
This is, of course, not known apriori. As a matter of fact, it is the possibility 
of just such a radial variation which is under investigation here. The results 
to follow will fortunately justify this seemingly bold assumption. 
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When equation (12) Is integrated radially it becomes 

R„ 


f & II 
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2rr u w , 9 

Qc 100 


+ pC p (T 2 -T x ) = fp P ou 2 (R 2 2 - R x 2 ) + 


2 

2 . 4 „3 
w ) ttC 


4 4 

V- R 1 

< - 4 ) 


The numerator of the first term, is obviously the total heat transfer exchanged 

with the surroundings „ f „ 

$ » | 2nr $ dr 

w R X J w 

With this definition the energy equation may be rewritten as 


o R, u) 2 f R? 2 " 

f + p VV T i> - 5 p <-f-> 1 - <r> • 

' pR i U)2 l L i 2-1 

1 + i- l*2_ fx + 200 <_£_ ) 2 1 1 + (5 1 ) 

20 -3 0 9 V /V R 2 J 

TIC J l p(JCfe J t 


where the heat transfer, enthalpic, and kinetic energy terms are easily identified. 


C, ANALYSIS OF TURBULENT FLOW 

The basic difference between the laminar and turbulent flow analyses 
is the way in which the dissipation (viscous) term in the energy equation is 
formulated. The general form of this term in dyadic potation is 

t: V V (15) 

where t is the stress dyadic and V is the velocity vector. For face seal flow 
this quantity may be approximated by 

T:VV - V <H> +T rz H (17) 

When the flow is laminar then 


T : W * 


, tat, 2 + (S£) 2 ' 


which appears in equation (4), For isotropic turbulent flow it has been shown that 
(Ref, 1) 

rZ 8Z dZ / dv (19) 

\ / 3Z/ 


which means that 


T: W - T 


@z dZ 


H+®: 
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where 


= T, 


8 z 0 z 


= 0.0225 p (V) 


1/A 
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r 7 

(^y-) 

c/2 


T 1/4 


( 21 ) 


Following the same procedure used in the laminar analysis, the integrated 


form of the energy equation for turbulent flow is 
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" 3T 

c, + pC 
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C p Sr J 
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udZ = r yj t 


ez + 7 J' 1 : 
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(Pq) dZ 


( 22 ) 


Noting that Q - I 2rrrudZ, the working form of the energy equation can be written 
as cr 

c 


h p ^ q q j r 2 

+ “= r WT. + P„|u(v) dz (23) 
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2rrr dr 


e z 


w 




The velocity distributions to be used when evaluating the last term of this 
equation will be the same as used in Reference 1 


u T 


v = 
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v 2 = r uj 


r 17-71 

[ x - 1/2 ( tt! J’ 


c/2 < z < c 


The pressure gradient from the same reference is 

1/4 


dp _ 835 / Wv 2 0.0195 ,cUJ. U) 3 x 

d? " 816 P r ( T> “ ~ Q 


(24) 


(25) 


(26) 


The evaluation of the last term in equation (23) leads to 

T z , - I u( ’-' )2 ‘ ,z ■ - 1-2090 i ( L ~) 3 fy 2 IttIpc^) 2 
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( 28 ) 


It should be observed that this rather complicated expression is really a simple 
function of the radius of the form 


yj u<v> 


2 d z - Ar 5/4 + 
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(29) 


Returning to equation (23) and integrating with respect to r yields 
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Carrying out the integration assuming that A and B are. independent of the radius 
yields 

q + pQC (T -T ) = 0.0225 p( v ) 1/4 (2n ou ) . 

W p C 1 
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71 
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4) (R 2 19 ' 4 - R x 19 ' 4 ) + 2rr A (f) (R 2 13 / 4 - Rl 13 / 4 ) 


+ 27 B (i) (R 2 2 - R t 2 ) 

Rearranging and replacing A and B with their earlier representations gives the 
final form 


(31) 
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(32) 


This equation plays the same role for turbulent flow which equation (14) plays 
"in -laminar flow. 


D. CONCLUSIONS 


The implications of equation (15) may be most conveniently assessed by 
first considering the case when Q=0. The equilibrium temperature of the fluid 
trapped in the clearance space is determined by the equation 
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(33) 


The term q , which contains the ambient temperature and the fluid film temperature, 
w 

is determined by heat transfer characteristics of the seal support and surroundings. 

Sirre the leakage rate through the seal is by design small, the temperature 
level within the seal even with non-zero leakage is expected to be very nearly equal 
to the equilibrium temperature determined by equation (33). This means that the 
temperature change across the seal is approximated by 

„ RJ0 2 

(T 2 -T x ) - f p (-4") |l -(^) | (34) 


P C . 


R 1 2 


So it would seem that the level of the temperature within the seal fluid is mainly 
established by equation (33), with the external heat transfer characteristics of 
the seal support system as the major controlling influence. This temperature level 
is then perturbed by through- flow leakage which does not substantially change the 
level, only the distribution of the temperature within the seal around that level. 

Calculations of the temperature rise predicted by equation (34) for typical 
operating conditions indicates that this rise is ordinarily small, the order of 
10°F or less. The most important temperature for the evaluation of the fluid 
properties is, therefore, the temperature level, calculated from equation (33). 
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It is now clear that the sealant leaks through the seal nearly isothermally 
relative to the no-leakage equilibrium temperature (the temperature level). This 
supports the isothermal assumption which was made during the analysis. The analysis 
itself has provided the means of estimating the isotherm to be used when selecting 


p, and p. 

Furthermore, it appears unnecessary to investigate the temperature distributions 
in misaligned, eccentric, wavy, or other seal configurations, since the temperature 
variations caused by these aberrations are undoubtedly small and inconsequential 
relative to the operating temperature level. 

The arguments presented above apply with equal validity to laminar or turbulent 


flow. The heat transferred in turblent flow with Q = 0 is 

r r i9/4 ir 

<^ =0 = 19 3TE 2 2 (R i”> 1 -<rJ> I 0 - 0225 » ^ IM 


0.0225 p v 




Estimates of the temperature rise associated with this term using typical seal data 
indicates values of the order of 1.0 F or less, again validating the isothermal 
assumption. 
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NOMENCLATURE 


c 


C 


P 


k 

P 


Q 


H 



r 


R. 


T 

T. 


u 

V 

z 


M* 

V 


p 


U> 


T 

zr 


T z0 

w 


seal clearance 
fluid specific heat 

fluid thermal conductivity 
fluid pressure 
supply pressure 

volumetric leakage rate 

total local heat transfer flux from seal 

total heat transfer from seal 

radial coordinate 
inner seal radius 

outer seal radius 

temperature 

temperature at inner radius 

temperature at outer radius 

radial velocity component 
tangential velocity component 
axial coordinate 
viscosity 

P7p 

density 

angular velocity of seal surface 
radial shear stress 

tangential shear stress at the wall 
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III. A STUDY OF REVERSED FLOW REGIMES IN PARALLEL FACE SEALS - LAMINAR AND 

TURBULENT FLOW | 

A. INTRODUCTION 

In an earlier investigation (Ref. 1) it was noted that when the net 
leakage rate was zero or nearly zero the radial velocity was directed outward I 

i 

along the rotating disk and inward along the stationary disk. The hydrodynamic 1 
features of this multi-directional flow were not explored in the earlier study, 
other than to suggest that a zero net leakage rate did not necessarily mean that 
there was no exchange of sealant fluid between the inner and outer boundaries of 
the seal. 

The analysis which follows is a more detailed investigation of the nature of 
these reversed flow regimes. 

B. STREAM FUNCTION - LAMINAR FLOW 

The analysis of the reversed flow regimes is most conveniently carried 
out with the aid of the stream- function for the stream surface in this axisymmetric 
flow. In Reference 1 it was shown that the continuity and momentum equations were 
satisfied by a stream- function of the form 

ty(r,|) * r u) Re c W(£) + U(fj) (1) 

where, to a first order approximation, 

U = U 1 = 2 r - 3 5' (2) 


W'S' w 


1 







(3) 


In the above equations £ ■ ^ where z is measured from the rotating surface 
toward the stationary one. The values of the stream-surface at the bounding surfaces 
are 


y (r,o) - 0 

* < r >» - - &• 


w 
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C. CASE OF LAMINAR NET RADIAL OUTFLOW (Q jfr 0) 


In this case we seek a stream-surface separating internal and external 
fluids which branches from the stream-surface for the stationary surface, \[/(r,l). 
For this purpose equation (1) will be written in the form 

W(l) - W(!) + p[u(l) - U(i)] = 0 (5) 

where 


B = — ^ 

p 2 3 

2p(r U) ) c tt 


( 6 ) 


From the initial investigation it is known that reversed laminar flows occur 

when 



3»Q 

p(r o)) 2 c 3 n 



(7) 


or 


- 0.01666 < B < 0.01111 (8) 

The shape of the stream-surface separating the fluid coming from the inside and 
the fluid coming from the outside is determined by obtaining the roots of equation 
(5) for P in the range given relation (#) . Figure 1 shows the result of such 
an evaluation for 0 < £ < 0.01111 (Q > 0) • 

In a practical sense, after having designed a seal and determined its leakage 
rate using the analysis given in Reference 1, one would next calculate P with 
r = R^. If P < 0.01111, then for R^ < r < R^, P will be in the range given by 
expression (7) and reversed flow will be present. Figure 1 allows one to calculate 
the shape of the separating stream-surface, and in particularly at what radius it 
branches from the stationary surface. 

A typical plot will look like that shown in Figure 2. If the separating stream 
surface branches within the seal itself as shown in Figure 2, no external fluid 
will reach the inside. If, on the other hand, this stream-surface branches at a 
radius less than R^ there will be an exchange of external and internal fluids 
even though the net flow is radially outward. 

D. CASE OF LAMINAR NET RADIAL INFLOW (Q < 0) 

In this case the separating stream-surface branches from the rotating 
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surface, so we seek zeros of the equation 


W(6) +'p U(6) - 0 m I 

with 

- 0.01666 < 0 < 0. | 

Whether there is an exchange of inner and outer fluids is determined, as 
above, using in this case Figure 3. If it is determined that the separating j 

stream-surface branches at a radius greater than as shown in Figure 4, there 
will be an exchange of inner and outer fluid. j 

E. CASE OF ZERO NET RADIAL FLOW (Q = 0) 

The case of no-net flow is special situation which cannot be treated 
the same way as the two previous cases. In this case the streamf unction reduces 
to 

%lr (r,£) = r cd ReC W(|) (IQ) j 

From equations (4) it is clear that the z axis, as well as the stationary and 
rotating disks, lie on the same stream-surface, From Reference 1 it is I 

obvious that W(£) is roughly parabolic in the range 0 < £ < 1.0 so that i|r is ' 

multivalued at any given radius. It is not difficult then to show that the , 

radial velocity vanishes on a surface | = 0.475 for all radii. With this I 

information the general streamline picture can be sketched (see Figure 5). 

The solution given by equation (IQ) is actually valid for disks with Rj*G ] 

as shown in Figure 5. When R^ > 0 the stream- surfaces shown in Figure 5 for the 
region r < R- are, of course, meaningless. The flows shown in Figure 2 and 

1 i 

5 being the limiting case of Figure 2, where the separating stream-surface has 
merged into the rotating surface. 

Once again there is an exchange of inner and outer fluids, the separating 1 

stream-surface being the one which just reaches the inner radius before turning 
around and returning to the outside. The number value of this particular stream 
line is obtained by inserting the coordinates of the turning point into equation 
(10), i. e., 

ilr . « i|r(JL, 0.475) = R Re C W(0.475) = -0.00262 Rju Re c 

Separating Y 1 1 1 

an 

The inflow and outflow leakage rates are equal in magnitude and given by the equation 
q = 0.00262 (2TR ) R <D Re d XU) 
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It should be noted at this point that when the separating stream-surface 
branches from the seal surfaces, as in the earlier cases studied(see Figure 4), 
the fluid exchanges are determined in exactly the same way as just described. 
These exchanges will not be equal in magnitude, of course, but their algebraic 
difference will be equal to the net flow rate Q. 

To find the turning point coordinates of the stream-surface just reaching 
a radial boundary in the cases where Q 0, one would calculate g from the 
formula 




1 + 4 <| - A) 


(13) 


setting r equal to or in 


A * lLkL-2 

" 3 2 

fJT c (ru) 


(14) 


Inserting the proper r and £' into equations 1, 2 and 3 yields the number value 
of the streamline bounding the exchange flows. Using the relations given in 
equation (4) the exchange rates of the inner and outer flows can be calculated. 


F. TURBULENT FLOW 

r 


Since the equations of motion for the turbulent case are more con- 
veniently solved in integral rather than differential form, the streamline technique 
of laminar flow will not be used with turbulent flow. An alternative approach 
is to calculate the flow rate through a cylindrical area whose height is measured 
from the seal surface and determine at what distance from the surface the net 
flow-rate equals zero. In this way it is possible to construct a functional re- 
lationship between the separating stream-surface and the radial coordinate. 

From reference 1 we have the following relationships for the radial velocity 
component and the pressure gradient 


and 
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These equations may be combined to give 
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It should be noted that z in this case is measured from the stationary surface. 
Defining the net flow-rate through the cylindrical surface of height z as 




<J “ / udz 
o' 


( 20 ) 


the position of the stream- surface is determined by the zeros of the expression 
z 
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For this purpose the following integrals are required: 
0 < z < c/2 
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c /2 < z < c 
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G, CASE OF TURBULENT NET RADIAL OUTFLOW (Q < 0) 

It was determined analytically in reference 1 that reversed flow, would 
occur when 


-27.7 < (- 


36 HQ ( pr me 


p ttc (r (D) 
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2 M. 


3/4 

) < 23.6 


(24) 


The zeros of equation (21) for values in this range are the only ones which need 
be considered* They are most conveniently determined with the aid of a computer. 
Figure 6 shows the computed results for th = positive portion of the range given 
in (24). In order to be consistant with the laminar results, relation 24 is re- 
written as 2/4 

~ . « ,or oo Cs . ^ 

(25) 


0.385 < 8 ( P r — ■ ) < 0.328 

~ p ~ 


nr 00 c 3 / 4 

and the product 3 (£— — a is plotted as a function of | = 1 - z/c. The stream- 

V> 

surface for these conditions looks like that sketched in Figure 2. 

Ho CASE OF TURBULENT NET RADIAL INFLOW (Q < 0) 

Since the separating stream- surface for this situation will look like the 

one sketch in Figure 4, the zeros of the following equation must be found 
z z 

0 = /udz - udz 

The results of this computation are plotted in Figure 7* 


(26) 
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I. CASE OF ZERO NET RADIAL FLOW 


The case of zero net flow once again requires special treatment. Referring 
to Figure 5, the exchange flow is determined by first locating the separating stream- 
surface, On this surface u * o so that £ is found from equation (17) by setting 
u = Q «= 0 , i. e. , 


G 


835 

816 ( c/2 ) 


Using the function F and G from reference 1 yields 


(27) 


z . - 0.55 

t/c 


( 28 ) 


Since the inflow and outflow exchange rates are equal, the exchange rate can be 


calculated from 


z = 0.55c 

t 


q t 


2 TT R 


1 / U(R 1 


where 


u(R 15 z) 


R„ R, u; 

) 


5 z) dz 
2 


2 C 2 } 2 ' 


T 0z 

UJ 


[ G( ^2> - Mf F <^2> 


(29) 


(30) 


T 0z (R l ) = °* 0225 P (V) 
UQ 


1/4 


R w 

<-r> 


1/4 


L c/2 

Carrying out the indicated integration yields 


(31) 


q = - 0.01696 (2 ttR^ ujc 2 

which is the equation for turbulent analogous to equation (12) for turbulent 
flow. 

J. CONCLUSION 

The means for determining whether reversed flow conditions result in an 
exchange of inner and outer fluids have been developed. When there is an exchange, 
a technique for dividing the net flow-rate into individual exchange rates is 
outlined. The results are restricted to the case where both inner and outer 
boundaries are continually immersed in the same fluid. 


. K (JO 

£ (— ) 
2 2 ; 


”1 


(32) 
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Because the special case of no net leakage (Q~0) results in a simple stream- 
surface geometry, the equations for the exchange flows can be deduced without 
reference to other the specifics of the seal operation. These equations (12) 
and (32) and Figures (1,3, 5, 6) provide a complete description of the reversed 
flow phenomena. 
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NOMENCLATURE 


A 

= 

36 nQ 

3 2 

p rrc (r <») 

c 

3E 

Axial clearance 

q t 


Exchange leakage rate for Q«0 

Q 

= 

Net leakage rate 

r 

= 

Radial Coordinate 

R 1 

= 

Inner Radius 

R„ 



Outer radius 

2 


c^> 

Re 


V 

u 

= 

Radial velocity 

U 1 


First approximation to U 

W 1 

= 

First approximation to W 

z 

= 

Axial coordinate measured from non-rotating disk 


P - 23 

2p(r («) c rr 

|i = viscosity 

V = |i/p 

I = 1 - z/c 

| = Coordinate of streamline turning point at seal radial boundary 

p = Density 

T z q = Tangential wall shear stress 
w 


\|r = Stream function 

U) = Angular velocity of rotating disk 


1 
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APPENDIX F 
SEAL DESIGN III 
BY 

H. J. Sneck 

ECCENTRIC FACE SEAL WITH A TANGENTIALLY VARYING FILM THICKNESS 
I. INTRODUCTION 

Previous investigations of the effects of geometry on the performance of 
face seals have considered such aberrations as misalignment waviness in the 
tangential direction, and film thickness variations in the radial direction (Ref.l). 
In each case, however, it was assumed that the seal faces themselves were centered 
on a common centerline. It has been suggested by Findlay (Ref, 2) that the 
an omolous pumping effects observed in face seals might be attributed to the ec- 
centricity of the seal faces which occurs when their centers are not on a common 
centerline. Preliminary analytical and experimental studies of the effects of such 
eccentricities indicate that the proposed pumping effect is indeed initiated with 
the introduction of an eccentricity (Ref. 3). In the analysis to follow the effects 
of tangential fluid film thickness variations and eccentricity will be investigated 
in order to determine how, and to what degree, they influence face seal performance. 



II. ALIGNED FLAT SEAL WITH ECCENTRICITY 

Before launching into the analysis of the seal with a tangentially varying 
fluid film thickness it should be demonstrated that the geometric effect of 
eccentricity alone cannot cause a face seal to perform in any unusual way. 

Figure 1 shows the geometry of the problem. The velocity of a point in the 
lower surface relative to the origin is given approximately by the formula 


VqSsC 


e r + re Q 


The tangential velocity of the fluid in the clearance space (Ref.l) is given by 


V = r 0 (^) 


( 2 ) 


The equation of motion governing the fluid flow within the clearance space 
has been shown to be (Ref.l) 


|E- - P si 

or r 



( 3 ) 
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Integrating equation (3) with the aid of equation (2) yields for the radial velocity 
distribution 


-v 2 2 4 

u “ h + f (r,9)z + g(r ,9) 


dr v 2^' 2 ' 12 


(4) 


Because of the seal eccentricity, the radial velocity must satisfy the boundary 
conditions 


(5) 


u = efi S in (§ - 9) , z * c 
d = 0 , z ■ 0 

This results in a radial velocity distribution of the form 


2 . 2 4 3 

, 2 . - cz v dp ,r a v ,z - c z 

u = <■— — ■> ,'S- ( — ITU 


) + ea (”) Sin($ - 9) 


(6) 


The integral form of the continuity for this case is 
c 


If 


[■ d_ 
Sr 


(ur) dz =0 


(7) 


This results in a differential equation for the pressure distribution 

d / rc 3 &\ m 3 , p« 2 c 3 r . . eqc Sin ft - 9) 

d? ( l*i dr J 5 ( fzTJ ■ + . 2 

The net leakage rate is computed from equation (6) using the formula 
.'j2nr e 


( 8 ) 


Q = r T f udz | d 9 (9) 

6T Lo J 

The separating force is obtained from the integration of equation (8) for p 


which is then inserted in 

R_ 


■II 


p rdrd9 


( 10 ) 


o R. 


Since the determination of Q and F both require 0 integrations from 0 tp 2D , it 
is obvious that the calculated flow-rate and separating force will be the same as 
the concentric seal due to the periodicity of Sin (§ - 9) . The mathematical reason 
for this is the assumption of a uniform clearance e as will be demonstrated in the 
analysis to follow* 


III. ALIGNED WAVY SEAL WITH ECCENTRICITY - LAMINAR FLOW 
A. Velocity Distribution 
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If the upper surface is allowed to be wavy, as shown in Figure 2, while the 
lower surface is eccentric and rotates, the tangential velocity distribution cah 
be represented by 

v ■ tq ^1 * J Y (11) 

where h is a function of © only. When this equation is combined with equation 
(3) and the result when integrated is 


u 


2 2 
z__ dR prq 

^ dr 


> 3 . .4 1 

z , 1 z 

’ 3h + 12 .2 

h J 


+ F(r ,9) z + G(r ,0) 


The boundary conditions to be applied to fhis equation are 
u = oe Sin ( $ - 0) , z ■ 0 


u = 0 


z = h 


so that the radial velocity distribution finally takes the form 


u * * (z2 ' hz > Sr 


dp. p r ft 
dr 

+ -ft e Sin ($ - ©) 


zh 

L' ^ 


2 3 - 4 

, £ z_ . .L. 

2 3h 12 ,2 , 

h J 


M] 


( 12 ) 


(13) 


(14) 


B. Pressure Distribution 

The integrated form of the continuity equation for this case is 
h h h 


Ji: 


(ur) dz + 


+ /rfl dZ + f 

qJ o J 


dw 


dz = 0 


(15) 


Sr "" T | r d© ' | dz 
Using the boundary conditions that w (0) = w(h) =0, this equation reduces to 

h h 


\ d_ 

dr 


. I , f dV 

r J u<,z + J Se 

L QJ J O 


dz = 0 


(16) 


The second term in equation (16) can be simplified as follows 
h h 


d_ 

d© 


/ vd2 ' J 


dv , , /uS dh 

ae dz + v(h) ae 


(17) 


with the aid of Leibnitz Rule. Since v(h) = 0 
h h 




vdz 


(18) 


This simplification, along with the velocity distributions given by equations 
(11) and (12), result in a partial differential equation for the pressure 
distribution 
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(19) 


d /I*L. ^R\ - 3 / Pfl 2 h 3 r N , e.Qh Sin (§ - 6) rn dh 
"Sr '■l^x Sr* 5 ^ 12> ; 2 + 2 39 

This equation, of course, reduces to equation (8) when h = c and also consistent 

with the earlier analysis of Reference 1 when e =0. Recalling that h is a function 

of 6 only, equation (19) may be easily integrated to yield 

P - P rV + Sin 9 - (0 + (3i4»» § - *nr +■ 


2h' 


( 20 ) 


( 21 ) 


The arbitrary functions of integration are determined using the pressure boundary 
conditions 

P » , r = R x 

P “ P 2 » r " r 2 

The pressure distribution satisfying these conditions is given by 
P - P 2 - “ pn 2 (r 2 -R 2 2 ) +(^4) Sin (I - 0) (r - Rj 

Vi / 


+f?<* 2 -V>f+ “* 2 


. jfB 


*1 


R_ 


P - P 
1 r 2 


3 2 2 t, 2v 

20 PO ( R 1 " R 2 ^ “ 2h 3 


2 2. Sh 

R 2 ' S9 


Sin (ft - 0) (R t - R 2 )V 


C. Leakage Flow 

The net leakage-rate is computed from the equation 
_2rr h 


•JJ 


udzrd0 


( 22 ) 


(23) 


Performing the first integration with the aid of equation (14) gives 

; . 2rr 


= f f-ia a + rd0 

o J L > & r J 
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a » 

J f(6) <36 

2T 

-I f 

L^n R 2 / 


[ P 1 - P 2 + lo <* 2 (S 2 2 - E l 2 >] ^ 


+ 6 v n eh Sin (# - 6)(R 2 - d0 (24) 

If it is now assumed that the surface waviness is approximated by the film thick 
ness equation 

h = c (1 + a Cos k 0) (25) 

the formula for the leakage rate is 


6n$n R„ 


3 _ „ . 3 2,_ 2 _ 2v 

C P 1' P 2 + 20 pfl (R 2 - R 1 > 


(1 + | a 2 ) 


, 6(J. n eca(R 2 -R^)Sin $ 

2 TT 


Cosk9 Cos 0 d0 


The total leakage rate consists of two parts, one due to the waviness alone and 
one due to eccentricity and wavihess together. The waviness contribution is same 
as found earlier (Ref. 1). The eccentricity contribution is seen to exist only 
if there is a "once-per-rev" (k=l) distortion of the type observed by Nau(Ref. 4). 
Then equation (26) becomes 


R = 1 


*Mn R 0 


■P + — 
r 2 ^ 20 


m 

2 /tJ 2 _ 2. /1x 3 2 N 

pn (R 2 - r ! ) ( 1+ 2 a ) 


+ 3 |i q e(R 2 “R 1 ) ca Sin $ j (27) 

Whether the eccentricity cases inward or outward pumping depends on 
i) the direction seal rotation 

ii) the phase angle of the waviness relative to the line of centers,#. 
Since it is quite difficult to determine the phase angle # , it is not surprising 
to find reports of both inward pumping effects and outward pumping effects exhibited 
by the same seal. A very logical explanation is a shift in# between disassembly 
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and reassembly, or a change in rotational direction. 


D. Torque 


Since the tangential velocity distribution is the same as the one used earlier 
(Ref. 1) the torque is unaffected by eccentricity, being given by the formula 


( 1 + f 2 ) 


neglecting terms of order a . 


E. Separating Force 


The separating force is computed from equation (10) using the pressure dis- 

tirbution given by equation (22). It will be noted that integration with respect 

X 

to 0 will eliminate terms containing — ^ ^ if the film thickness is described by 
equation (25) . The remaining terms fa^l into two categories, those which contain 
S and those which do not. The terms which do not contain e do not contain h 
either, so their contribution to the separating force is the same as a concentric, 
aligned, flat-surface seal. 

The eccentricity effects the separating force through the terms containing 

Sind - 8) . This expression can be expanded using equation (25) as follows 

■ ,2 : ' 


Sin($ - 9) _ .Sin $ Cos 9 - Sin 9 Cos? 


(1 - 2a Cos k9 + 3a 2 Cos 2 k0 + *V ) 


If k - l 


a i 

/ Sin 
h 


(4 - 9) ,r, ~ -2a Sin & 

2 dQ = 2 

c 


„ 2 „ -2a Tt Sin 9 

Cos 0 d9 = — “ 2 — 

c 


The total separating force using this approximation for k = 1 is 

* *i 2 > +1 2 (1 §r> ( r 2 2 ' h 2) + 


P 2' P 1 " < ^" )(R 2 2 " R 1 2) 




12|j.a&TTeR- 3 Sin 0 I , f R 2 
+( — ^ ■> * 


R 2 tn R 2^ R 2 -R 1 

2 “ ( 4 } 
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1 1 
1 + 


l -<r 1 /-r 2 )‘ 


1 - R /R R R 

In R 2 /R 1 2 ( R^ *R^ 


k - 1 


The geometric factor contained in the last term is shown in Figure 3. 


F. Numerical Example 

Equation (27) for the net leakage is comparable to that reported earlier 
(Ref 2) with the exception that the inertia contribution is included here. Using 
the same data as assumed in Ref. 2 for a sample calculation, we have 

R x - 1.125" 

R 2 = 1.75" 

» = 4.95(10" 5 ) psi/sec 

-3 

c * 10 inches 
a * 1.0' 

e - 0.003" 

n- (i° 3 ) £> 

$ * 90° 

The contribution to the total leakage rate due to the eccentricity alone is given 

ky f\ r. Z' D _» C-Jr. & "1 


** " 2 


fle ca (R 2 ~R 1 ) Sin $ 
^n(R 2 /R 1 ) 


tt(10 3 ) (0.003) (10 -3 ) (1,75 - 1.125) 


( 1 * 25 . ) 

n.i25 ; 


= 5.86 (10 ^) in 3 /sec. 


= 15.86 (10" H ) (3600) (24) 

= 50.6 in 3 /day 

The load carrying capacity contribution due to eccentricity done is given by 


12 ^ Oe R 2 3 Sin* I , f R 2 ] , 


— I i - - (>—) 
2 V R J 

z ! 2 


1 " > 
R 2 


215 


“1 - Rj/R 2 “ 


Jin R 2 / R i 

_ 


1 ^ 

2 <5? *” <l£> - 


2. 1 - <W 


12(4.95) (10 -5 ) 10 3 (~) (0.003) (1.75) 


10 


-6 


0.00382 


3.82 lb. 


IV . ALIGNED WAVY SEAL WITH ECCENTRICITY - TURBULENT FLOW 


A. Velocity Distribution 

For turbulent flow the tangential velocity profile will be taken to be the 
1/7 power law measured from the surfaces. Using the velocity matching condition 
that the mid- channel velocity is one-half the rotating surface velocity, the 
tangential velocity profile is then given by 


= r a 


1 t*—\ 

2 Mi/2 ; 


1/7“ 


0 < z < h/2 


(32 ) 


v 


2 


2 


(—) 

K h/r 


1/7 


h/2 < z < h 


(33) 


It will be recalled that these profiles are consistent with the conce'pt' th'at the 
tangential shear stress is onstant across the clearance space assuming that the 
shear flow follows the "law of the wall", Ref. 1. It is easily verified that 
equations (32) and (33) also have the same slope at mid-channel with the result 
that the velocity profile is "smooch" all across the clearance space. 

The radial velocity profile can now be computed from the equation of motion 


u T 


Z0 


w 


- v £f (r,0) + * - £ J v 2 dzj 

+ J £ v 3 dz - Jv ^ dz + g(r,0) 


(34) 


which was derived in Ref. 1 with the aid of the concept of isotropic viscosity. 

This radial velocity profile must satisfy the boundary conditions given in 
equation (13). Evaluating the arbitrary functions f(r,0) and g(r,0) yields 


ur 


Z0 




Sin (4 - 0) + 
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( 35 ) 



When the velocities given by equations (32) and (33) are inserted into 


equation (35) the radial velocity distribution takes the form 


u T 


z6 


w 


F <g72> - f <?> <!> G <i 72 > 


where 


^h72> 




= (—) 

( h/2> 


1/7 


1/7 


1 " 8 ( h/2^' 


and 


1/7 


G ^h/2^ = ( h/2^ | + 12 + 2 ( h/2^ 


ri+i 

I 4 8 

[• 


<i^2> 


, 0 < z < h/2 

, h/2 < z < h 
8/7 


1 _ ( ^ 9/7 . i _ ( ^. ) 10/7 1 

18 'h/2^ 90 'h/2' J , 0 < z < 

. ■ 9/7*1 

Ht!> J 


<E72> 


i/ 7 r 


35 ]__ 

36 " 90 


h/2 
h/2 < z < h 


(36) 


(37) 


(38) 


With the exception of a sign change due to a change in direction of rotation, 
equation (36) is identical to one reported earlier for turbulent flow (Ref. 1). 

It should also be noted that z is measured from the rotating surface in this analysis j 
whereas it was previously measured from thV stationary surface. 


B, Pressure Distribution 

The radial pressure gradient is determined from the continuity equation (16) 
with the aid of equations (32), (33), (36), (37) and (38) . The result is 
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(39) 


^ - m » & * + *vk - °-° 225 f 

r 




j£ 


<f> 


L (h/2) 9/4 


0.022S P (v) 1/4 e fl (||) 


ra 

2 


1/4 


h/2 


5/4 


Sin (§ - 0) 


Integrating this equation with respect to r yields the general solution 

P - !rt § (?> 2 - 0-0225 f (v) 1 ' 4 • 


816 £ v 2 

r- 3/4 


11/4 


ii_ # 


i ( — ) 

7 ■49'' 


(h/2) 9/4 | 

■ Iff ' 4 ' 


(h/2) 


5/4 


~ 4 a.0225 p ;(y)^ /4 ea sin; ($ - e) ; 


r 7/4 + f f (9) r 3/4 + g (0) 


5h 

& 


(40) 


The arbitrary functions f(0) and may now be evaluated using the pressure boundary 
conditions 


p = P^ @ r = R^ 


p = P 2 @ r = R 2 


(41) 


The pressure distribution which results is given by the formula 

835, ox ,.(k 2, 2 „ 2, 


p-? 2 -fg<e$vv> 


+ 0.0225 (f§>- <fj-) f- 


1/4 


3/4 _ , 


fj (f) 


.9/4 


r r 3 / 4 k2 3 / 4 1 

v> 3 /4 r 3/4 
R 1 R 2 J 


(h/2) 

p r p 2 + lit f (f>V - *L* 


(rll/ 4 . R2 ll/4 ) I 


0.0225 (f£) (fj) f (v) 1M 


FI 3 / V) 

2 1 

* A I / 


(h/2) 


9/4 


11/4 11/4. Sh 

(R 1 R 2 7 ae 


+ 0.0225 p (v) 1/4 (^)(||) ea Sin (i-0) 


, 3/4 

II 


(h/2) 


575 


(r 7/4 - R 2 7/4 ) 
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• 3/4 _ 3/4 
r - R 2 

R 3/4 3/4 

h r 2 


0.0225 p (v) 1/4 (y)(||) en Sin (t - 9) 


(h/2) 


(R L 7/4 - R 2 7/4 > 


When 6=0 this equation reduces to the earlier result given in Ref. 1 with 
the sign of the angular velocity reversed. 


C. Leakage Flow 

The net leakage-rate is calculated from equation (23) using for the velocity 
profile the formulations given by equation (36) and (39). The resulting leakage 
rate is given approximately by 


p.p + M5 (£) |v.) r A)1[i + 

*1 z + 816 VI 2 ' 1 V/ J L 

mm ’ ' J- Z 


45 2 
64 a 


0.0195 


1/4 <f> (R, 374 - Rp 3/4 ) 




R 7/4 R 7/4 *1 
R 1 ~ R 2 

R 3/4 R 3/4 

R 1 R 2 J O. 


h Sin (i - 9) d0 


neglecting terms of the order a . 

The last term in this equation Is the contribution due to eccentricity. As 
in the laminar case the integral in this term is non-zero only when k = 1 . In 


this case 


h Sin ( i - 9) d9 = n ca Sin $ ; k = 1 


so that the total leakage rate for k = 1 is 


p _p + 

1 *2 + 816 


, R 2 n , 2 f. A. 1 , 45 2 

( 2 )( 2 ^ 1 “ ( R > J 1 + 64 a _ 

•» y •» 


0.Q195 Aa 3 

TT 9 

C 


<t> <* 2 3/4 - -h 3/4 > 
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+ 7 (f) e r 2 ca rr Sin $ 


v /4 

i -(if) 

R 2 


1 - ( R> 
R 2 


Once again it should be noted that eccentricity alone is not sufficient to cause 
additional leakage, it must be coupled with a specific kind of surface waviness 
(k = 1). Even then the additional leakage can be outward or inWbrd depending 
on the phase angle $ or the direction of rotation. 

D. Torque 

The total torque exerted by the fluid on the rotating foace of the seal is 
given by 2 t R 2 

T -/ | T . r 2 rdrde (46) 

oJ RjJ 


T . r rdrde 

Z0 

w 


The wall shear stress used to obtain equation. (39) is the same as used in Ref. 1, 
namely r . »3/4-t 1/4 


T zQw - 0.0225 p (v) 1/4 (f) — J (47) 

If the film thickness relation given by equation (25) is expanded binomially 
before integration with respect to 0, it is clear that the torque is the same 
as for the concentric case, i.e., 

r o ^ i a -I 1/4 to I/. 


T = 0.00333T 


pVr 2 19 


*1 

1 - irk 
R 2 




neglecting terms of the order a . 


E. Separating Force 

The last two terms in the expression for the pressure distribution present 
the contribution of the eccentricity. The separating force is determined by the 
integral 2r R 


j/" 


p rdrd0 


When this integration is carried out the separating force for the concentric case 
will be retrieved with the addition of a term proportional to s. Since both 


eccentricity terms of equation (42) contain Sin($ 

(h/ 2) ' 


the integral of this factor 
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will be determined first. Using the film thickness relation given by equation (25) 


7 Sind - 9) d9 


Sin 4 Cos 9 - Cos 4 Sin 9 

j (h/2) 5/4 

<«/*> 5/4 J 

(1 + a Cos k&) 5/ ^ 4 


as 


2TT 


(c/2) 


5/4 


(Sin 4 Cos e: - Sin 9 Cos 4) (1 • ^ a Cos kS + ••••) d0 


5 

4 


.• a tt 


(c/2) 


5/4 


Sin 4 , 


k - 1 


(49) 


It should again be noted that only if k « 1.0 is there any contribution of ec- 

2 

centricity to the order of a . 

The total separating force for k = 1 and small a computed from equations 
(10) and (42) is approximately 


» * <v - y > +5 <it> § - y> 


2 ,„ 2 


+ 2 TT 


P 2 " P 1 " 816 ^^ (R 2 " R 1 * 


3/4 _ 3/4 

R 2 R 1 


4_ 11/4 11/4 

11 (R 2 R 1 ) 


R 3/4 

K 1 2 2 

2 < R 2 - Y 5 


a e 


R 2 fi 


3/4 


10 w 30 n 


1/4 


+ 0.0225 (y-) (g) p (v) 


<T-> 


Sin 4 


L (c/2) 


5/4 


< u 


r 1 -(^) 

R 2 


7/4 


1 - R 

( R > 

2 


3/4 


11/43 


1 - <R > 
R 2 


4_ 

15 


^1 

1 -<rT> 


15/4 


1 - 


w 7/4 

1 

1 ' 


1 " *1- 
( rT> 


3/4 


1- R 1 ' 

' 


} ; * - i 


( 50 ) 


The geometric factor contained in the last terra is shown in Figure 4. 


V. CONCLUSIONS 


From the analysis of both laminar and turbulent flow it is clear that ec- 
centricity contributes to both the flow-rate and separating force if the surface 
is wavy. 

For small amplitudes of waviness relative to the nominal clearance (c) this 
contribution is proportional to the product of the eccentricity and waviness, which 
means that the eccentricity will have no effect unless the surface is also wavy. 

The analysis also indicates that a "once per rev" waviness (k ='l) is the primary 
contributor to the pumping effect. 

An unusual feature of the pumping due to eccentricity is the effect of rota- 
tion and phase angle ($). Depending on the direction of rotation and the magni- 
tude of $ , the puling contribution due to eccentricity can be positive or nega- 
tive augmenting diminishing the leakage due to pressure difference, centrifugal 
forces and waviness. Since in a practical sense the phase angle and eccentricity 
are uncontrollable factors, one would expect to see some random variations in 
leakage and separating force between apparently identical seals which were actually 
assembled differently. 
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ROTATING, FLAT, ECCENTRIC SURFACE 
FIGURE 2 

SEAL GEOMETRY FOR A WAVY SURFACE AND 
RADIAL ECCENTRICITY 
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FIGURE 4 

GEOMETRY FACTOR FOR TURBULENT FLOW 


EQUATION (50) 






NOMENCLATURE 



u 

v 


V 

ft 

w 

z 

e 


e 

v- 


ratio of the waviness amplitude to the nominal clearance 
nominal clearance 
unit vectors 

total force in clearance space tending to separate the seal surfaces 
local film thickness 
an integer 

pressure in clearance space 

pressure at the inner radius of the seal 

pressure at the outer radius of the seal 

volumetric leakage - rate 

radial coordinate 

inside radius of seal 

dutside radius of seal 

total torque on seal face 

temporal mean radial velocity component 

temporal mean tangential velocity component 

seal surface velocity vector 

temporal mean axial velocity component 

axial coordinate 

seal eccentricity 

tangential coordinate 

fluid viscosity 


v 

P 

T 

§ 

ft 


= Wp 

= fluid density 

= shear stress 

= phase angle between eccentricity line of centers and waviness reference 

line 

= seal surface angular velocity 
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APPENDIX G 


LITERATURE REVIEW 
A # J # ORSINO 

J. A. REILLY & H. T. ROBISON 

I. INTRODUCTION 

An important part of the technical work under this contract is the review 
of available literature on seals and sealing mechanisms. Such a review permits a 
continuously expanding background upon which to evaluate and appraise the new 
fundamental knowledge as well as the general broad technology of seals and sealing 
mechanisms . 

The original work in this area, completed in February, 1963, produced a 
rather extensive bibliography of all the available data and information on seals 
and sealing mechanisms* Over 2200 documents were considered for abstract and 
entry into the bibliographies which were published as part of the final report for 
contract NAS 7-102. Subsequent phases of NAS 7-102 continued the review of the 
available literature and the results have been issued as a part of the subsequent 
final reports. It is expected that a separate addendum or supplement will even- 
tually be required to formally update the original bibliographies. 

A considerable portion of the material which was reviewed was of a purely 
theoretical nature and was essentially an engineering data gathering to obtain 
information for use in conjunction with the various main task efforts. This ma- 
terial appears referenced in the appropriate sections where it is treated more 
extensively rather than being listed here in abstract form as part of this litera- 
ture review* The remaining references obtained are included in this appendix and 
are arranged alphabetically according to author; or, as in the case of NASA docu- 
ments, were arranged numerically by accession number. Some minor departure from 
this pattern may be present but this is not considered sufficiently critical to 
detract from the usefulness of the material listed* Where possible the listed 
reference contains an abstract prepared from a reading of the text* For some 
references the author* s prepared abstract was used. In several instances no text 
was available, however, the document was included since it was considered to con- 
tain information of interest in some area of seals and sealing mechanisms. 

II. DESCRIPTION OF REVIEW 

In the original work more than twenty-two separate sources from the English 
and American literature were completely reviewed. These sources are fully 
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identified in NASA report N63-19597, which is the bibliography on the open litera- 
ture developed under the original study. The bibliography of over 700 abstracts 
from the Defense Documentation Center (formerly ASTIA) is contained in NASA report 
N63-19596. 

The review conducted for contract NAS 7-434 was concerned with a further 
updating of the bibliographies by searching the more recent issues of the various 
primary sources. The references chosen, along with suitable summaries or ai 
stracts, are contained in this Appendix as supplements to material presented in 
previous reports. 

The following sources were reviewed, to the dates indicated, as a continu- 
ation of previous reviews conducted on this study contract and contained in NASA 
reports N63-19596, N63-19597, X64-16399 and X65-105B2. 

Source 

Applied Science and Technology Index 
ASM Review of Metal Literature 
Chemical Abstracts 
Engineering Index 
International Aerospace Reports 
Nuclear Science Abstracts 
Science Abstracts (A) (Physics) 

Science Abstracts (B) (Elec. Eng'g.) 

Scientific and Technical Aerospace Reports 
Technical Translations 
U.S. Government Research Reports 

III. BIBLIOGRAPHY SUPPLEMENTAL MATERIAL 

The supplemental material developed in the course of updating of the bibli- 
ographies is present here. The "Supp" numbers used follow sequentially from pre- 
vious reports. Additional references, primarily of a theoretical nature, are 
contained in the text presented as a report of the work effort for each major task. 


Date 

January 1968 
January 1968 
January 1968 
January 1968 
January 1968 
January 1968 
January 1968 
January 1968 
January 1968 
January 1967 
January 1967 
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SUPPLEMENTAL BIBLIOGRAPHY OF 
OPEN LITERATURE ON SEALS* 


*(See introductory text of this appendix for 
a full description. of this supplement and 
other bibliographic material issued under 
NASA Contract NAS 7-102 jand NAS 7-434, 
Phase I) 
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Supp. 


200 

Anonymous NSA21-67 3225 

Loads to Deform Elastomer 0 Seals in Grooves and Between Plain Flanges 
Atomic Weapons Research Establishment, Aldermaston, (England), October 
1966, 27 pp. (AWRE-0-82/66) 

The loads necessary to deform 0 seals in various types of grooves and 
between plain flanges have not been extensively investigated. As there 
is a recurrent demand for this information, a series of practical tests 
was carried out. The results of 210 tests are presented on a nomogram. 


Supp. 201 

Anspach, W. F. Eng. Ind., July 1967 

Hydrofluorocarbon High Temperature Integral Fuel Tank Sealants 
Rubber Chem. & Technology, Vol. 39, No. 4, Pt. 2, September 1966, 
pp. 1200-1214 

Interim solution to critical fuel tank sealant problem in March 3 air- 
craft; high temperature sealant based on mixture of high and low 
molecular weight hydrofluorocarbon polymers has solids contents of 
85$ and cures at moderate temperatures; it is one-component system 
stable for several weeks at ambient temperatures; sealant exhibits good 
room temperature and high temperature. Mechanical properties both 
before and after aging in hydrocarbon fuel and fuel vapor at temperatures 
up to 500F adhesion to metal substrates is good. 


Supp. 202 

Appleton, B. Eng. Ind., January 1967 

Design Considerations for Polyurethane Sealing Devices 

Nat'l Conference on Fluid Power, Proc., Vol. 19, 1965, pp. 79-83 

Designer must examine carefully what is ultimately required of seal 
and what fringe benefits seal material can yield to product to which it 
is to be fitted; chemical and chemical properties required in seal, cost 
of material and cost of processing must be considered. Types of seals 
made from polyunethane resins include: 0-rings, U-cups, V-rings, flange 

packings, cup packings, rod wipers, scrapers, bellows, diaphragms, and 
special balljoint and bearing seals. 


Supp. 203 

Bazzarre, D. F. , Spain, R. G. , Withers, J. C. Eng. Ind., July 1967 

Metal Clad Elastomeric Cone Seals for Liquid Rocket Propellants 
Products Finishing, Vol. 31, No. 1, October 1966, pp. 66-77 

Development of seals resistant to storable liquid rocket propellants was 
motivated by two facts, i.e., (a) that current storable propellants 
differ widely in composition and necessitate elastomeric vulcanizater to 
secure reasonable compatibility, and (b) that propellants are new fluids 
and existing elastomers are eliminated rather quickly as elergetics of 
these fluids advance; fabrication and preliminary evaluation of limited 
quantities of metal clad 0 -rings, using variety of metals in thickness of 
1 to 10 mils; determination of techniques for preparing conductive surface 
on 0 -ring and application of metals by electroplating; evaluation of metal 
clad 0-rings in N 2 O 4 and in acetone at 75F. 
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Supp. 204 

Bellanca, C. F . , Salyer, J. 0. Eng, Ind. , July 1967 

Effect of Liquid Rocket Fuels and Oxidizers on Elastomeric 0-Ring Seals 
Rubber Chem. and Technology V. 39, No. 4, Part 2, September 1966 
pp. 1215-1221 

Tests of selected elastomeric and compliant materials were made at 73F 
with N 2 O 4 , CIF 3 , and 90 $ H 2 O 2 (oxidizers) and with mixed hydrazines, 
pentaborane and hybaline A-5 (fuels); at 73F only poly tetraf luroethylene 
was resistant to all fluids; butyl and ethylene propylene rubber appear 
satisfactory for long-term storage in presence of mixed hydrazines and 
hybaline A-5 ; fluorocarbons appear best suited for use with pentaborane 
and hydrogen peroxide. 


Supp. 205 

Billington, J. J, , Toronchuk, J. P. Eng. Ind., December 1967 

Characteristics of Self Energized Hydrostatic Seals for Rotating Shafts 
Eng. Inst. Canada Trans. Paper EJC- 66 -MECH 7, May 1966, 10 pp. 

Theoretical and experimental study of controlled leakage shaft se#ls 
with fixed gap maintained between rotating and stationary members is 
described; fixed gap is by self-energized hydrostatic actuators; effect 
of various geometrical parameters on seal performance was investigated 
theoretically; good correlatibn was obtained between experiment and 
theory when operating in low viscosity fluid (water); effects of 
pressure-induced seal face deformation and effects of thermal transients 
were studied. 


Supp. 206 

Billington, J. J. , Fitzsimmons, T. E. NSA67-21 17432 

Final Report on Shaft Seal Development 


A coolant pump shaft seal development program is summarized; particular 
reference is made to the final endurance test series on hydrostatic 
self-energized, finite clearance seals. The endurance program resolved 
the seal geometry and materials. 


Supp*, 207 

Bupara, S. Q., Walowit, J. A., Allen, C. M. IAA67-22705 10-15 1658 

Gas Lubrication and Distortion of High Pressure Mainshaft Seals for 
Compressors 

B.H.R.A., Int. Conf. on Fluid Sealing, 3rd; Cambridge, England, April 
1967, Paper B3, 24 pp., 7 ref. 

Theoretical analyses and static experiments have been performed for 
high-pressure, gas lubricated face seals of the type used for mainshaf t 
seals for jet-engine compressors. The analyses consider the combined 
effects of seal-element distortion apd gas flow through the seal interface. 
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Eng. Ind., June 1967 


Supp. 208 

Davies, A. R. , 0* Donoghue, J. P. 

Lubrication of High-Pressure Face Seals 
ASME-Paper 66-WA-LUB-7 for Meeting November-December 1, 1966, 8 pp., 

12 Refs. 

Theoretical analysis of lubrication; conditions of hp face seal is 
presented which uses very simple treatment for estimating effects of 
elastic and thermal distortion, c. Agreement jhetween tkepryuand experimental 
results obtained for carbon seal is good; from results obtained, it is 
concluded that there is inherent stiffness in seal and this enables 
leakage to be controlled by load clamping seal If aoesrtGge 
of speed on temperature and problems of hydrodynamic pressure generation 
under seal face are discussed; it seems that from order of pressure used 
(1000 lb. ft./sq. in.) hydrodynamic effects are not of major importance* 


Supp. 209 

Decker, 0. Eng. Ind., September 1967 

Advances in Dynamic Seal Technology 

ASME-Paper 67-DE-50 for Meeting May 15-18, 1967, 16 pp. 

Need for advanced seal concepts has been brought about by increasingly 
severe requirements imposed by applications such as supersonic trans- 
port, dynamic space power machinery, nuclear applications cryogenic 
machinery, and liquid-metal technology; conventional sealing techniques 
prove inadequate to meet challenges improved by such applications; 
recent trends and technological developments &r,e 

sealing problems; some dynamic seals are described in applications 
employing liquid metals, gas, and move standard organic fluids for both 
ground, aerospace and space turbomachinery./ 

Supp. 210 

Decker, 0. IAA67-28881 14 2328 

Advances in Dynamic Seal Technology 

ASME, Design Engineering Gonf. & Show N.Y. , May 1967, Paper, 

67-DE-50, 17 pp., 21 Refs. 

Description of recent trends and technological developments in solving 
advanced sealing problems. A variety of dynamic seals is described in 
applications employing liquid metals, gas and more standard organic 
fluids for ground, aerospace and space turbomachinery. 


Supp. 211 

Dickenson, R. W. NSA21-67* 14356 

Bearing and Seal Requirement for Large Liquid Metal Cooled Reactor 
Systems 

Atomics International, Canoga Park, Calif., MTI-66TR66, pp. 8, 1, 1-6 

The requirements of bearings and seals in pumps for high power output 
liquid-metal cooled reactors are outlined. Present and potential uses of 
sodium lubricated bearings are described. The need for a reliable, 
economical seal for mechanisms in sodium reactors is emphasized. 
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Supp. 212 

Eibert, P. NSA21-67 8002 

Nuclear Reactor Coolant Circulator with a Liquid Shaft Sealing 
To Siemeas-Schuckerverke Aktiengesellschaft German Patent 1, 218, 076, 
June 2, 1966, Filed September 29, 1962 

The sealing comprises an annular gap portion between the shaft and 
the housing filled with a sealing liquid, and a second gap portion 
filled with a sealing gas, the latter portion being located adjacent to 
the interior of the reactor vessel. In the middle of the annular portion 
which is filled with gas there is a chamber connected to a storage tank 
containing clean reactor coolant gas. The liquid sealed portion and the 
gas sealed portion are separated by another chamber connected to the 
interior of the housing, into which the sealing liquid and the sealing 
gas escape and which is provided with an exit through which the gas- 
liquid mixture is drawn off. 


Supp. 213 

Gerlach, T. Eng. Ind. , January 1967 

Self-Adjusting Seal with Fluid Friction 

Polska Akadenia, Prace Instytutu'Maszyn Prezeplywouych, No. 22, 1965, 
pp. 15-33 

Self-adjusting seal, with constant clearance between seal elements 
independent of pressure and fluid viscosity, enabling frictionless 
transfer of dynamic loads in beariig test rigs; equation of equilibrium 
of sealing element is analyzed; dynamic properties of seal are examined 
for constant load and for periodically changing load occuring in test 
rig operation; sufficient conditions for attenuation of vibrations are 
determined. Results of experimental tests verify theoretical 
assumptions. In Polish. 


Supp. 214 

Green, J., Levine, N. B. , Sheehan, W. Eng. Ind., July 1967 

Elastomers Resistant to Rocket Propellants 

Rubber Chem. & Technology, Vol. 39, No. 4, Pt. 2, September 1966, 
pp. 1222-1232 


Resistance of polymers to hydrazine - type fuels, hybaline A-5, nitrogen 
tetroxide, and fluorine - containing oxidizers; butyl and ethylene 
propylene rubber are found compatible with hydrazine ; 1 carboxy-nitrose 
rubber is recommended for long-term exposure to N 2 O 4 at 165F; fluoro- 
carbon plastics and silicone rubber are most suitable for use with 
oxygen difluoride, chloryl f luoride-tetrafluoro-hydrazine mixture; 
teflon, Kel-F81, and" uncured nitrose rubber gum are suitable for limited 
applications in chlorine trifluoride. 
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Supp. 215 

Halliger, L. Eng. Ind. , June 196 r 7 

Sealing of Anti-Friction Bearings 

T2 Fuer Praktische Metallbearbeitung, Vol. 60, No. 4, April 1966, 
pp. 207-218 

Important features and properties of common shapes of protective seals 
are described and factors to consider in their selection indicated; 
development of grease-charged bearings with built-in seals for maintenance 
free bearings; examples of seal combinations which proved satisfactory 
under plant conditions. In German. 


Supp. 216 

Hamaker, J. Eng. Ind., January 1967 

New Materials in Mechanical Face Seals 

Nat'l Conference on Fluid Power, Proc., Vol. 19, 1965, pp. 58-65 

Fundamentals of seal operation, basic seal types, available and 
advantages and limitations inherent in seal's design and materials of 
construction are discussed; special emphasis is given to materials 
introduced in recent years; materials for secondary seals and seal rings. 


Supp. 217 

Hawkins, R. J. NSA21-67 8007 

Improvements in or Relating to Sealing Means 

To Babcock and Wilcox, Ltd., British Patent 1,045,910, October 19, 1966, 
Filed September 11, 1964 

A method is described for making nuclear reactor seals by means of which 
a flow path for coolants, from the reactor core to heat exchangers, is 
formed. Ducts which provide for circulation of coolants are given a 
sealing contact by sliding blacks of graphite. The blocks are held in 
place,; against movement caused by thermal expansion', by springs, pins, 
or toggle arms, but may be moved vertically or horizontally. 


Supp. 218 

Hazard, R. E. Eng. Ind., September 1967 

Testing Seals for Aerospace 

Fluid Power Int. , Vol. 32, No. 373, April 1967, pp. 38-40 

New sealing technique developed in view of increasingly higher 
temperatures and pressures required to work aerospace hydraulic equip- 
ment; design of testing machitie at British Aircraft Company at Stevenage, 
England, for evaluating dynamic performance of seals and to control 
temperature and pressure of hydraulic fluid, fluid leakage part seals 
and piston movement; types of seals investigated were metal, plastic 
dhd combined metal and plastic. 
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Supp* 219 

Hooke, C. T. , Leais, D„ J t , O^Donaghue, J. P. 

Elasto Hydrodynamic Lubrication of O-Ring Seals 

IME Proc, , Vol. 181, Pt. 1, No. 9, 1966-1967, pp. 205-210 

The paper presents a theoretical solution to the ideal problem of an 
0-ring seal constructed of an isotropic elastic material moving at a 
constant axial velocity along a cylinder. The net leakage past the seal, 
which is the difference between the fluid flow across the seal when it 
is "pumping” against the pressure and that when the seal is "monitoring" 
or moving with the pressure, is shown to increase with sealed pressure. 


Supp. 220 

Hummel, B. L. - Editor IAA67-22610 09-15 1509 

Seals: Reference Issue (3rd Edition) 

Machine Design, Vol. 39, March 9, 1967, 219 pp. 

The Design Data section contains completely revised and updated informa- 
tion which clearly and concisely describes how to select and apply seals 
and sealants. The all-new Product Directory is set up so that the: 
reader in search of information about manufacturers of seals and sealanti 
can easily find exactly which company can fill his needs, 

Supp. 221 

Johnson, D. F. Eng. Ind . , January 1967 

Controlled Gap and Segmental Seals 

Nat f l Conference on Fluid Power, Proc. Vol. 19, 1965, pp. 39-43 

Performance materials and construction of controlled gap and segmental 
carbon shaft riding type seals are discussed; seals can operate to 
1400F for short periods of time and to 1200F for long periods of time; 
they are easy to install over shaft and have capability to handle axial 
tolerance or shaft movements. 


Supp. 222 

Kinsbury, J. E. IAA67-15788 04-18 0642 

New and Proposed Cryogenic Temperature Resistant Plastics 
Society of Automotive Engineers, Aeronautics and Space Engineering and 
Manufacturing Meeting, Los Angeles, Calif., Oct. 3-*7, 1966, Paper 660638, 

5 pp. 

Discussion of the need for cryogenic- temperature-resistant plastic 
material in the space program; the selection of a plastic material for 
cryogenic application in dependent on the cryogenic fluid with which it 
will come in contact. Problem areas discussed are insulation of 
cryogenic materials, adhesives used on plastics, use of seals and 
gaskets and the expulsion of cryogenic propellent in a zero gravity 
environment. It is believed that although the research being done in this 
field at present appears to have limited applications in the aerospace 
industry, it may be of great significance for the future. 
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Supp 


. 223 

Koehler, H. Etig. Ind., January 1967 

Apparatus for Sealing in Plastic Forming Machines 

U. S. Patent 3,238,566, March 8, 1966 

Seal for sealing space between two members having sheet of plastic 
material interposed there between when one of the members is moved into 
engagement with sheet; continuous, seal receiving groove around outer 
edge of one of members and yieldable sealing member mounted in groove, 
sealing member being formed upwardly angled flat surface with projects 
above upper edge of said groove, and centrally disposed, continuous 
passage provided in one side wall of seal member. 

Supp. 224 

Lawrie, J. M. , Summersgill, J. Eng. Ind., March 1967 

Design Specifications for Hydraulic Seals 

Fluid Power Int., Vol. 31, No. 368, November 1966, pp. 330-334 

Information required for seal design specifications; categories con- 
sidered include applications, fluid to be sealed, temperature, pressure, 
relevant dimensions, additional technical : i nf or^ , special, life, 
function limitations, petmissible leakage and form of lubrication. 


Supp. 225 

Lindsey, H. 

Mechanical Seals: Carbons Key Role 

Chemical Engineering, Vol. 74, No. 5, February 27, 1967, p.,160 

Carbon f s unique qualities have made carbon sealing ring an indispensable 
part of the mechanical seal. Now, a denser, more impervious carbon is 
widening mechanical seal applications, and a new type of components is 
reducing costs. Carbon - bedause it is self-lubricating, resists 
chemical attack excellently, withstands large and rapid temperature 
change, and can be machined readily by conventional methods - finds wide 
application in the petroleum and petrochemical industries which handle 
many gases and liquids that cannot provide the lubrication that is 
necessary between moving surfaces. 


Supp. 226 

McManus, S. P., Piken, S. Eng. Ind., July 1967 

Elastomeric Seals for Army 1 ^ Lance Missile 

Rubber Chem. &: Technology, Vol. 39, No. 4, Pt. 2, September 1966, 
pp. 1233-1246 

Missile propulsion system requires dynamic and- static seals which must 
operate under extremely severe conditions for short periods; resin-cured 
chlorobutyl rubber meets requirements for low and high temperature 
operation and retains properties on exposure to propellants; comparisons 
with vinylidene f luoride-hexafluoropropylene copolymer (Viton). 19 Refs 
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Moore, B. C. , Bergquist, L. E., Clements, P. R. Eng. Ind., July 1967 
Cryogenic Low-Pressure Seal 

J. Vacuum Science & Technology, Vol. 4, No. 1, January- February 1967, 
pp. 45-46 

Communication to Editor describes development of simple, reliable metallic 
seal aluminum or stainless steel flanges at low pressures and cryogenic 
temperatures; self metal wire gasket is retained in rectangular groove 
in aluminum spaces. 
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Mowers, R. E. IAA67-31990 17 2864 

Static and Dynamic Seals Used on Rocketdyne Engines 

So c. Auto. Eng. Aerospace Systems Conf., Calif. , June 1967 

The paper reviews typical sealing concepts and materials used through- 
out the propellent feed systems and pneumatic and hydraulic control 
systems of liquid propellent rocket engines. Emphasis is placed on the 
sealing concepts utilized to compensate for extreme conditions. 
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Murray, F. S. NSA21-67 14345 

Bearing and Seal Materials for Liquid Metal Fabrications 
Mechanical Technology Inc., Latham, N.Y. , MTI-66TR66, pp. 6.3.1-26 

Flame and plasma-sprayed hard coatings for bearing and seal components 
of process-fluid lubricated bearing systems are evaluated. The 
effectiveness of some of these coatings in protecting the surfaces 
against sliding damage in inert environments has been demonstrated in a 
number of practical applications. All of the test results show that 
these coatings can be used effectively for liquid metal-lubricated 
systems. The coatings simplify the design and fabrication of the bear- 
ing and seal components, but have the disadvantage of lack of adequate 
quality control. Necessary quality control procedures for the use of 
coatings must be developed to ensure a reliable product. 
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Nau, B. S. IAA67-27889 14-15 2324 

Hydrodynamic Lubrication in Face Seals 

B.H.R.A. Int. Conf. on Fluid Sealing, 3rd, Cambridge, England, April 1967, 
Paper H3, 48 pp. , 54 Refs. 

The literature of fundamental studies of face seals is reviewed in 
detail with the emphasis on the conditions in the interface region. 
Attention is concentrated on hydrodynamic face seals. Among the 
subjects covered are the effects of surface roughness; wear-rates; 
leakage and inward pumping; and film thickness measurements. 
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Supp. 231 

Paxton, R* R* Eng* Ind* , December 1967 

Carbon and Graphite Materials for Seals, Bearings, and Brushes 
Electrochem Technology, Vol* 5, No 0 5-6, May-June 1967, pp* 174-82 

Developments of past ten years in carbon and graphite materials for 
seals, bearings, and electric motor brushes are briefly reviewed; 
properties of wide range of currently available materials are compared, 
and variety of applications are described indicating type of carbon or 
graphite which has been found suitable for each; bibliography of recent 
literature is included* 242 Refs 0 


Supp* 232 

Ruhl, F, F., Wendt, A* B* , Dalenberg, P, H, 

A New Approach to Face Type Sealing 

Lubrication Eng* 23(6), pp* 241-244, 1967, England 

The following 17 major factors affect face type seal wear: design, 

materials, speeds, temperatures, pressure, alignment, lubrication or 
lack of it, cooling, mechanical forces, chemical potential, abrasive 
foreign material, care in assembly, the medium being sealed, secondary 
mechanical factors, secondary chemical factors, possible electrical 
effects, and the interactions of all factors* These factors are 
appraised in relation to each other and in relation to the overall 
problem* 


Supp* 233 

Schnetzer, E* NSA21-67 14353 

Alkali Metal Bearing and Seal Development at Space Power and Propulsion 
Section 

General Electric Co*, Cincinnati, Ohio, Space Power and Propulsion 
Section, MTI-66TR66, pp* 6*2*1-46, Includes Appendices A through E 

Three major programs associated with the development of Rankine cycle 
power conversion technology for space applications were undertaken; 
low viscosity bearing stability investigation, development of dynamic 
seals, and preparations for a liquid metal bearing test, simulating the 
conditions of a space generator* All programs were supported by a 
corresponding analytical effort, presently geared towards predicting 
bearing-rotor response* Experimental results are expected in 1966 and 
1967* Screw seals, rotating channel and slinger squeeze seals were 
investigated in the dynamic seal program* Interface instability was 
identified as one of the major problems* It was resolved for the 
case of the slinger squeeze seal, which was successfully tested in 
potassium* 


Supp. 234 

Shevckenka, R* P* IAA67-22537 09-28 1560 

Shaft, Bearing and Seal Systems for a Small Engine 
Society of Automotive Engineers, Automotive Engineering Congress, 
Detroit, Mich*, January 9-13, 1967, Paper 670064, 10 pp. 
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Discussion of the key mechanical difficulties inherent in designing 
shaft bearing and seal systems for three high-speed turbine engines. 

A review is made of the bearing shaft and seal system of a family of 
large turbine engines which has demonstrated its high reliability, high 
efficiency, and low cost per unit of power as the thermodynamic, 
material and structural component technology has been improved and as 
the design power requirement has been reduced. Three different solutions 
to the problems discussed are presented; two for sophisticated, efficient 
engines for weapon system consideration, and the third for' a potential 
industrial, low cost, high-ef f iciency , moderately heavy engine. 


Supp. 235 

Smaley, E. M. Eng. Ind., March 1967 

Sealing with Gaskets 

Machine Design, Vol. 38, No. 25, October 27, 1966, pp. 171-187 

Design guide presented shown best practices for creating seal in 
gasketed joints and insuring seals longevity; two sets of parameters 
are surveyed, those that create seal and those that affect life in 
service, e„g., flange movement; gaskets considered are cork composition, 
cork and rubber asbestos fiber and rubber cellulose fiber and rubber, 
or rubber above; flat and grooved - face types of flanges are discussed. 


Supp* 236 

Stair, W. K* , Hale, R. H. IAA67- 22706 10-15 1659 

The Turbulent Visco Seal - Theory and Experiment 

B.H.R.A. Int. Confo on Fluid Scaling, 3rd, Cambridge, England, April 
1967, Paper $2, 32 pp«, 12 Refs. 

Experimental evaluation of the laminar viscosial for a number of seal 
geometries has shown that the Boon and Tal f s analysis is in better 
agreement with experiment than any other analytical approach available. 
It has been observed further that the seal performance improved and 
became dependent upon Reynold's number as turbulent conditions were 
encountered. 


Supp. 237 

Stair, W. R. IAA67-15356 04-15 0630 

Effect of Groove Geometry on Viscoseal Performance 

ASME Winter Annual Meeting and Energy Systems Exposition, New York, 

N.Y., November 27-December 1, 1966, Paper 66-WA/FE 28, 8 pp., 19 Refs. 

An experimental investigation of the effect of groove geometry on the 
performance of the grooved-shaft viscoseal is being conducted. Using 
distilled water as the sealed fluid, the concentric and eccentric 
performance of ten seals has been determined under laminar and 
turbulent conditions. An approximate method is devised to predict the 
empirical factors in the turbulent sealing equation. The findings reveal 
that the optimum seal geometry for laminar operation is not optimum for 
turbulent operations* The seal eccentricity has slight effect on 
performance, and the phenomenon of air injection may present a major 
problem in the application of the viscoseal. 


241 


Supp. . 238 

Ward, J. R. Eng. Ind. , June 1967 

Piston Seal Development for Oil-Free Compressors 

Naval Engrs. J.,Vol, 78, No. 6, December 1966, pp. 995-1002 

Report of experimental work at U. S. Navy Marine Engineering Laboratory 
leading to design of nonlubricated seal for air compressors intended 
for shipboard use, which successfully operated 1500 hrs. at discharge 
pressure of 4500 psi; seal consisted of sleeve made of 25 % glass-filled 
TFE and follower made of polyimide material; long term operation 
occured in chromium-plated cylinder bore; in later evaluations, follower 
material was charged to metal impregnated carbon-graphite. 


Supp. 239 

Wasil, T. J., McCleary, G. P. 

Sealing Corrosive Materials 

Lubrication Engineering, Vol. 23, No. 6, June 1967, pp. 234-240 

This paper discussed tungsten carbide, tungsten titanium carbide, 
and titanium carbide as applied to axial mechanical sealing materials 
for sealing corrosive fluids. It states that the corrosive resistance 
of the carbides is dependent mainly upon the compatibility of the 
corrosive media with the carbide binder material. The tungsten carbide 
grains are inert to most corrosive fluids. The general mechanism of 
the carbide corrosion is explained. Galvanic corrosion, concentration 
cell corrosion, erosion corrosion, direct chemical attack and oxidation 
are discussed in detail. 


Supp. 240 

Zuck, J., Strom, T. N., Ludwig, L. P., Johnson, R. L. , IAA67- 28797 14-15 
2327 

Convective Inertia and Gas Ingestion Effects on Flow Regimes of the 
Viscoseal - Theory and Experiment 

ASLE, Annual Meeting, 22nd, Toronto, May 1967, Paper, 61 pp., 39 Refs. 

An experimental and visual investigation of the viscoseal operating 
regions was conducted. Importance of convective inertia forces and 
changes in the functional relations between variables in the super 
laminar flow regime are discussed. Correlation between modified 
Reynolds number and sealing capacity improvement is found. 
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^SUPPLEMENTAL REFERENCES FROM 
SCIENTIFIC AND TECHNICAL AEROSPACE REPORTS (STAR) 

AND 

TECHNICAL ABSTRACT BULLETINS (TAB) 


#The major portion' of the references included in 
this listing are the result of government sponsored 
research and development. 
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N67-11171 AD 638972 

Air Force Systems Command, Wright-Patterson AFB, Ohio, Foreigh Tech. 

Div. 

Operation of Cascade Contact Seals 

A. M. Yevtushenko, 31 August 1966, 25 pp., Refs., (FTD-HT-66-281 ; 77-66- 
62243) 

Translation of Russian research: operation of cascade contact seals. 


N67-13998 

Atomic Energy of Canada, Ltd., Ottawa (Ontario) 

Performance Testing of Conical-Face High Pressure Rotary Shaft Seals 
R. D # Watson, November 1965, 36 pp., Refs. (AECL-2551) 

The performance characteristics of a special type rotary shaft seal 
have been determined at pressures up to 1000 lb ft/in 2 gage and speeds- 
at 1900 rev/min. The seal is made up of two rings, with one flat face 
ring in contact with the very slightly cone-shaped face of the other 
ring. Leakage rates for the seal were studied. 


N67-17506 

NASA Lewis Research Center, Cleveland, Ohio 
Screw-Seal Performance in Vacuum Environment 

J. C. Hudelson and L. P. Ludwig, Washington, NASA, February 1967, 

16 pp., Refs. (NASA-TN-D-3842) 

Non-contact seals are being proposed for use in space power systems. 
The non-contact seals offer reliable operational lifetimes of over 
10,000 hours. The leakage rate of an ether lubricant through a non- 
contact slinger pump and screw-seal assembly to a simulated space 
environment was measured. 


N67-19126 

Pratt and Whitney Aircraft, East Hartford, Conn. 

Development of Main Shaft Seals for Advanced Air Breathing Propulsion 
S ystems 

Semiannual Report No. 2, 1 January-30 June, 1966 

R. L. Thomas, A. J. Parks, and R. M. Slayton, 9 September 1966, 130 pp. 
(Contract NAS3-7609) (NASA- CR- 8 245 3 ; PWA-2879) 

Detailed analytical studies of all seal configurations to be tested are 
continued. Approval was received for the following seal concepts: 
orifice compensating hydrostatic face seal and face contact seal with 
bellows secondary. 
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N67-22711 

Royal Aircraft Establishment, Farnborough (England) 

Temperature Limitations of Polymeric Materials. Part 6: Elastomers 

G. Wood, April 1964, 40pp., Refs. (RAE-TN-CPM-52) 

Elastomers find many applications in aircraft, e.g., seals (static and 
dynamic), sealants, gaskets, 0-rings, and different properties are 
significant for particular cases. This report reviews the properties 
of the more important commercial elastomers and sealants. 


N67-23037 

British Hydro Mechanics Research Association, Cranfield (England), 
Information Dept. 

References on or Related to Screw-Type Rotary Shaft Seals 
March 1966, 14 pp„. Refs. (BIB-22) 

A compilation of abstracts related to screw-type rotary shaft seals of 
interest to hydromechanics is presented. Areas covered include screw 
, viscosity pumps, the sealing action of a viscous fluid, helical groove 
seals, high speed lubricating oil pumps, and screw seals in laminar 
and turbulent operation. 


N67-25487 

Pratt and Whitney Aircraft, East Hartford, Conn. 

Development of Compressor End Seals, Stator Interstage Seals, and 
Stator PivOt Seals in Advanced Air Breathing Propulsion Systems 
Semi-annual Report, 1 January-30 January 1966 

R. M. Hawkins, 20 July 1966, 249 pp„. Refs. (Contract NAS3-7605) 

(NASA- CR- 83 786 ; PWA- 2875; SAR- 2) 

The design of compressor end seals, stator interstage seals and stator 
pivot seals is discussed in detail. One-side floated shoe seals, two- 
side floated shoe seals, and thin strip seal designs are considered. 


N67-25842 

NASA Lewis Research Center, Cleveland, Ohio 

Improving Performance of Face Contact Seal in Liquid Sodium (400° to 
1000°F) by Incorporation of Spiral-Groove Geometry 

L. P. Ludwig, T. N. Strom, G. P. Allen and R. L. Johnson, Washington, 
NASA, May 1967, 39 pp.. Refs. (NASA-TN-D-3942) 

Conventional face contact seal performance was improved by incorporation 
of spiral-groove geometry. In comparison with conventional face contact 
seals, seals with spiral grooves had negligible leakage. The wear and 
contact patterns indicated that the spiral-groove seal operated with 
separation of the sealing surfaces, which is necessary for long life. 
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N67-26213 

Atomic Weapons Research Establishment, Aldermaston (England) 

Loads to Deform Elastomer 0-Seals (Toroidal Sealing Rings) in Grooves 
and Between Plain Flanges 

The loads necessary to deform O-seals in various types of grooves and 
between plain flanges have not been extensively investigated. This 
report describes and presents the results of 210 tests in the form of 
a nomogram. 


N67-31128 

Pratt and Whitney Aircraft, East Hartford, Conn. 

Development of Main-Shaft Seals for Advanced Air-Breathing Propulsion 
Systems 

Semiannual Report, 1 July-31 December 1967, 85 pp, (Contract NAS3-7609) 
(NASA- CR- 85908; PWA-2996; SAR-3) 

Four main-shaft seals for advanced gas turbine applications have been 
designed. The seals studied are an< of ifice- compensated hydrostatic face 
seal, an extremely pressurized hydrostatic face seal, a carbon-face- 
contact seal with a bellows secondary seal, and a carbon- face-contact 
seal with a piston-ring secondary seal. 


N67-31151 

Pratt and Whitney Aircraft, East Hartford, Conn. 

Development of Compressor End Seals. Stator Interstage Seals, and 
Stator Pivot Seals in Advanced Air Breathing Propulsion Systems 
Semiannual Report, 1 July-31 December 1966 

R. M. Hawkins and A. H. McKibbin, 20 January 1967, 94 pp. , Refs. 
(Contract NAS3-7605) (NASA-CR-85891; PWA-2995; SAR-3) 

The basic design concept for the '’OC" diaphragm thin strip seal is 
discussed. The final design of the one-side floated shoe seal is dis- 
cussed in some detail. The final designs of the vane pivot seals and 
test rig are also discussed. 


N67-31390 

Tennessee Univ. , Knoxville, Deptw of Mechanics and Aerospace Engineering 
Theoretical and Experimental Studies of Visco Type and Buffered Shaft 
Seals 

Semiannual Progress Report, 15 October 1966-15 April 1967 
William K. Stair and C. F. Fisher, Jr., May 1967, 25 pp., Refs. 

(Grant NsG-587) (NASA- CR- 85 893; ME-67-587-9) 

This report outlines progress made on the investigation of the viscoseal 
and the buffered seal. The experimental performances of fine grooved 
shaft viscoseal are presented. 
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N67-35568 

NASA Marshall Space Flight Center, Huntsville, Alabama 

Conference on Design of Leak-Tight Separable Fluid Connectors, Vol. 1 

25 March 1964, 340 pp. , Refs. Proc. hpld 24-25, March 1964 (NASA-TM- 

X-60234) 

Contains twelve papers given at the conference. 
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AD-379 053L Fid. 21/5 

Pratt and Whitney Aircraft, East Hartford, Conn, 

Component Propulsion Program for Future High-Performance Strategic 
Aircraft. Vol. XVI Seals 

Semiannual Report No. 3, 1 August 66-31 January 67 

R. C. Barningham, Contract AF33 (657) -14903, Proj. AF-61139A 

Confidential Report 


248 


